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Preface 




This book has one purpose: to be a time saver. In all the time spent 
writing and debugging programs, one of the most frustrating oc- 
curences is to need a formula for the particular program at hand. 
More times than not, precious time is spent looking for that elusive 
formula or its derivation. In this book I've tried to put in one single 
place many different types of programs with the hope that they will 
provide the necessary information you seek. Needless to say, 
many of the programs were direct results of my own interests. 
However, I've tried to add others from different fields to provide a 
suitable cross section. 

The biggest problem when writing this book was in hitting a 
happy medium between simple, mundane programs and more com- 
plicated, obscure ones. Ohms law, for instance, does not require a 
separate program. Its usage is wide enough that the formula should 
be at one's fingertips. In the same way that many of us have become 
dependent on the hand calculator (Quick: what's 15 x 12?), myself 
included, we should not let the computer rule our existence down 
to the smallest detail. Remember, the calculator is an excellent 
time saver, not a crutch. Therefore, in this book you will see some 
basic indentity-type programs as well as some complicated 
single-purpose programs. Each, though, has been selected to pro- 
vide a source of information. With this information at your finger 
tips, you, too, will hopefully, have more time for using your 
computer. 



Dedication 

For My Very Organized Father 
From His Very Disor^nized Son: 

"Time has value, and we each spend 
that value for our own satisfactions/' 
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Chapter 1 
Numbers 



This chapter contains only seven different programs, but each 
program is useful in many different ways. The first three programs 
form the basis for many others, since complex operations and 
polar/rectangular conversions are used in many different facets of 
electronics. And, needless to say, where else do you find different 
bases but in computers? 

Fve tried to present the programs in such an order that the 
basic programs are first, followed by the more complex programs 
which may, or may not, use parts of the basic programs. In this way 
the later programs will be easier to understand, and you should be 
able to pull parts out as your own needs dictate. 
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POLAR/RECTANGULAR CONVERSION 

The polar/rectangular conversion is one of the most funda- 
mental programs around. However, as with most of the other 
programs in this book, it's not one that you can recite without some 
extra time for thought. Figure 1-1 shows a diagram that relates the 
various components of the conversions. 

To convert from polar to rectangular use the following 
equations: 

X = McosT 
Y = McosT 

where M is the magnitude and T is the angle in the polar notation. 
For rectangular to polar, use: 

M = VX 2 +Y 2 
T = arctan(Y/X) 

Polar/Rectangular Conversion Program 

5 'P/R CONVERSION 

10 CLS 

20 PR I NT "SELECT APPROPRIATE FUNCTI 

ON P"CHR$<94)" R (1) OR R"CHR$ 

(94) "P (2)" 
30 INPUT N 
40 IF N=l PRINT"ENTER MAGNITUDE, P 

HASE ANGLE"! INPUT M,T ELSE GOTO 

80 
50 X=M*COS(T*. 01745329) 
60 Y=M*SIN(T*. 01745329) 
70 PRINT" X="X" AND Y="Y:END 
80 PRINT" ENTER X AND Y" : INPUT X»Y 
90 M=SQR<Xt2+Yt2> 
100 T=ATN(Y/X>*57. 29578 
110 PRINT" MAGNITUDE^' M "AND PHASE 
ANGLE="T 
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Fig. 1-1. How the various components of polar/rectangular conversions 
relate to each other. 



Sample Program 

SELECT APPROPRIATE FUNCTION P ♦ R(l) OR R ♦ P(2) 

?1 

ENTER MAGNITUDE, PHASE ANGLE 

?3,60 

X = 1.5 ANDY = 2. 59808 

SELECT APPROPRIATE FUNCTION P+R (1) OR RtP (2) 

?2 

ENTERXANDY 

?5, 12 

MAGNITUDE = 13 AND PHASE ANGLE = 67.3802 
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COMPLEX OPERATIONS 

The complex operations comprise the four basic arithmetic 
functions: addition, subtraction, multiplication and division of 
complex numbers. The equations for the functions are: 
addition - (x t + ty,) + (x 2 + iy 2 ) = (jc, + \) + (y t + y 2 ) i 
subtraction - (x t + iy : + (x 2 + iy 2 ) = (x, + x 2 ) + (y t + y 2 )i 

-(x 2 -(y, + ty,) 

multiplication- (x, + iy 1 )(x 2 + ty 2 ) = M/M 2 T ; 

division- (Xj + ly x )/(x 2 + ty 2 ) = M/N^ 
where M and T are the respective magnitudes and phase angles 
from the rectangular to polar conversion. 

Complex Operation Program 

5 'COMPLEX OPERATIONS 

10 CLS 

20 PRINT" SELECT THE APPROPRIATE CO 
MPLEX OPERATION" CHR$ ( 13 ) TAB ( 20 ) 
"ADDITION 1"CHR*(13)TA 

B< 20) "SUBTRACTION 2"CHR$< 
13) TAB (20)" MULTIPLICATION 3" 
CHR$< 13) TAB( 20)" DIVISION 
4" 

30 INPUT N 

40 PRINT" ENTER FIRST COMPLEX NUMBE 
R (REAL* IMAGINARY)" : INPUT FR»F 
I 

50 PRINT" ENTER SECOND COMPLEX NUMB 
ER (REAL? IMAGINARY)" : INPUT SR* 
SI 

60 IF N=l GOTO 190 

65 IF N=2 GOTO 170 

70 FM=SQR(FR+2+FI*2) 

80 FT=ATN(FI/FR)*57.2957B 

90 SM=SQR<SR*2+SIt2) 

100 ST=ATN(SI/SR)*57.2957BsIF N=3 
GOTO 150 

110 AR=FM/SM 

120 AI=FT-ST 

130 BR=AR*C03(AI*. 01745329) 

140 BI=AR*SIN(AI*. 01745329) iGOTO 2 
10 
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150 AR=FM*SM 
160 AI=FT+ST:GOTO 130 
170 BR=FR-SR 
180 BI=FI-SI:GOTO 210 
190 BR=FR+SR 
200 BI=FI+SI:GOTO 210 
210 PRINT" ANSWER <REAL, IMAGINARY) 
— "BR?BI 
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Complex Operation Sample Program 

SELECT THE APPROPRIATE COMPLEX OPERATION 

ADDITION 1 

SUBTRACTION 2 

MULTIPLICATION 3 

DIVISION 4 

?3 

ENTER FIRST COMPLEX NUMBER (REAL, IMAGINARY) 

?3,5 

ENTER SECOND COMPLEX NUMBER (REAL, IMAGINARY) 

?6,-9 

ANSWER (REAL, IMAGINARY) 

--63 3.0 

SELECT THE APPROPRIATE COMPLEX OPERATION 

ADDITION 1 

SUBTRACTION 2 

MULTIPLICATION 3 

DIVISION 4 

?4 

ENTER FIRST COMPLEX NUMBER (REAL, IMAGINARY) 

?-12,3 

ENTER SECOND COMPLEX NUMBER (REAL, IMAGINARY) 

?4, — 5 

ANSWER (REAL, IMAGINARY) 

-1.53658 1.17073 

These two samples, multiplication and division, are the more 
difficult. They use the previous polar/rectangular conversion. 
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COMPLEX FUNCTIONS 

The four complex functions— absolute value, square, recip- 
rocal, and square root— also use the polar/rectangular conversion. 
The respective formulas are: 

absolute value - |x| = M 

square -x 2 = W /2T 

reciprocal — = -rr- /-T 

square root=\/x"= ± (n/M /T/2) 

with M and T the magnitude and phase angle from the 
rectangular/polar conversion. 

Complex Function Program 

5 'COMPLEX FUNCTIONS 

10 CLS 

20 PRINT" SELECT THE APPROPRIATE CO 
MPLEX OPERATION" CHR$< 13) TAB (20) 
"ABSOLUTE VALUE 1"CHR*(13)TA 
B( 20)" SQUARE 2"CHR*( 

13)TAB(20)" RECIPROCAL 3" 
CHR$ < 1 3 ) TAB ( 20 ) " SQUARE ROOT 
4" 

30 INPUT N 

40 PRINT"ENTER COMPLEX NUMBER ( REA 
L» IMAGINARY)" ! INPUT R, I 

50 M=SQR<R*2+I* 2) 

60 T=ATN<I/R)*57. 29578 

70 IF N=l PRINT"ANSWER — "M:END E 
LSE IF N=2 GOTO 120 ELSE IF N=3 
GOTO 100 

80 AR=SQR(M)*C0S(T/2*. 01745329) 

90 AI=SQR<M)*SIN(T/2*. 01745329): 
GOTO 140 

100 AR=l/M*COS(-T*. 01745329) 

110 AI=1/M*SIN<-T*. 01745329) :GOTO 
140 

120 AR=R«2*C0S<2*T*. 01745329) 

130 AI=R*2*SIN<2*T**01745329) sGOTO 
140 

140 PRINT" ANSWER <REAL» IMAGINARY) 
— "AR»AI 
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Complex Function Samples 

SELECT THE APPROPRIATE COMPLEX FUNCTION 
ABSOLUTE VALUE 1 

SQUARE 2 

RECIPROCAL 3 

SQUARE ROOT 4 

?2 

ENTER COMPLEX NUMBER (REAL, IMAGINARY) 

2 —3 

ANSWER (REAL, IMAGINARY) — 5.0 - 12.0 

SELECT THE APPROPRIATE COMPLEX FUNCTION 
ABSOLUTE VALUE 1 

SQUARE 2 

RECIPROCAL 3 

SQUARE ROOT 4 

?4 

ENTER COMPLEX NUMBER (REAL, IMAGINARY) 

?-3.5,7 

ANSWER (REAL, IMAGINARY) -- 2.37973 - 1.47076 
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CONVERSION TO BASE 10 

In this program, which converts a number to base 10, I've 
actually treated the number as a string and then operated on each 
integer individually. For example, the first of four digits is 
multiplied by the base to the third power (1234 -- 1 x 3 3 , 2 x 3 2 , 3 
x 3\ and 4 x 3° = 58). The basic equation is: N, n = X B"" 1 + 
X^B"- 2 + X_ 2 B"- 3 ...X 2 B+X ° " 

Conversion to Base 10 Program 

5 'BASE CONVERSIONS 

10 CLS 

20 INPUT "ENTER NUMBER TO BE CONVE 

RTED TO BASE TEN" 5 X# 
25 INPUT "ENTER BASE OF NUMBER TO 

BE CONVERTED" ;B 
30 N=LEN<STR$<X#))-i 
35 M=N-1 
40 FOR 0=2 TO N 

50 S=S+VAL < M I D* < STR* < X# > » Q » 1 > > *B* M 
55 M=M-1 
60 NEXT Q 

70 S=S+VAL(MID$<STR$<X#>,0>1)> 
80 PRINTXtt" BASE "B" EQUALS" S" BASE TE 

N" 
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Base 10 Conversion Samples 

ENTER NUMBER TO BE CONVERTED TO BASE TEN 

?3456 

ENTER BASE OF NUMBER TO BE CONVERTED 

?7 

3456 BASE 7 EQUALS 1266 BASE TEN 

ENTER NUMBER TO BE CONVERTED TO BASE TEN 

?8776 

ENTER BASE OF NUMBER TO BE CONVERTED 

?9 

8776 BASE 9 EQUALS 6468 BASE TEN 
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CONVERSION FROM BASE 10 

As with the previous program, the conversion from base 10 to 
a new base is handled by treating the individual digits as part of a 
string. Each digit is separated from the string and operated on with 
the new base and its appropriate exponent. A FOR NEXT loop is 
used to evaluate the entire number. During the loop, the individual 
digits are again stored in a string until the completion of the 
process. 

A number is changed to a new base by the following equation: 

N = X-I fir 1 -I R n ~ 2 ...IB n 

x n-1 n-2 n 

Conversion From Base 10 Program 

5 'BASE CONVERSION 

10 CLS 

12 CLEAR 1000 

15 P=l 

25 DIM Z*<300> 

30 INPUT "ENTER THE NUMBER TO BE C 

ONVERTED (BASE 10) "?N 
40 M=N 

50 INPUT "ENTER NEW BASE" 5 B 
60 FOR 1=0 TO 50 
70 X=INT(B*I) 
80 IF X>N GOTO 100 
90 NEXT I 
100 J=I 
110 1=1-1 
120 R=l 

130 M=INT(M-INT<B*I)> 
140 IF M=0 Z*(P)=STR$(R): GOTO 220 
145 IF M=INT(B*I> GOTO 160 
150 IF M<INT(BM) GOTO 180 
160 R=R+l:GOTO 130 
180 Z$(P)=STR*(R) :PRINTZ$(P) 
190 P=P+1 
200 1=1-1 
210 IF M>=INT(BU) GOTO 120 ELSE R 

=0:GOTO 180 
220 IF 1=0 GOTO 250 
230 FOR X=P+1 TO J 
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240 Z*<X)=STR*(A> 

245 NEXT X 

250 FOR Y=l TO J 

260 W*=W$+Z$(Y) 

270 NEXT Y 

280 PRINTN?"BASE 10 = "W*" BASE" 5 B 
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New Base Conversion Samples 

ENTER THE NUMBER TO BE CONVERTED (BASE 10) 

? 3112 

ENTER NEW BASE 

?7 

1 

2 



3 

3112 BASE 10 = 12034 BASE 7 

ENTER THE NUMBER TO BE CONVERTED (BASE 10) 

? 1234 

ENTER NEW BASE 

?2 

1 





1 

1 



1 





1234 BASE 10 = 10011010010 BASE 2 
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SIMULTANEOUS EQUATIONS 

Use of this program will allow you to solve for two unknowns- 
Each of the unknowns is independently solved by using the coeffi- 
cients of the terms in a matrix. Should the coefficients A°E - B*D 
= 0, the program will indicate a / error. This indicates that 
either no solution or a unique solution exists. To solve for X and Y, 
the following equations are used: 



X = 



1 c 


B | 


1 E 


D ' 


A 


B 1 


D 


E J 



CE- BF 
AE- BD 



1 A 


C | 


1 D 


E 1 


A 


B 


D 


E 



AF- DC 
AE-BD 



Simutoeous Equations Programs 

5 'SIMULTANEOUS EQUATIONS 

10 CLS 

20 PR I NT "TO SOLVE SIMULTANEOUS EQU 
ATIONS FOR TWO UNKNOWNS* " CHR*< 1 
3)" REARRANGE THE EQUATIONS INTO 
THE FOLLOWING FORMAT. " CHR*< 1 
3)TAB(20)"AX+BY=C"CHR*(13)TAB<2 

0)"DX+Ev=F" 
30 INPUT "ENTER THE COEFFICIENTS -- 

A»B»C»D»E,F" ;A»BjC*D>E»F 
40 X=(C#E-B*F)/<A*E--B*D> 
50 Y=<A*F--D*C)/<A*E--B*D) 
60 PRINT" X="X" AND Y="Y 
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Simultaneous Equations Examples 

In the following examples, the program will be used to solve 
for X and Y in the two groups of equations: 

2X + 3Y = -5 . 9.5X - 8.5Y = 3 
-5X-2Y = 4 ma 2X + 4Y = -5 

TO SOLVE SIMULTANEOUS EQUATIONS FOR TWO UN- 
KNOWNS, 

REARRANGE THE EQUATIONS INTO THE FOLLOWING 
FORMAT. 

AX + BY = C 
DX + EY = F 

ENTER THE COEFFICIENTS - A,B,C,D,E,F, 

? 2,3,-5,-5,-2,4 

X = -.181818 AND Y = -1.54545 

As proof: 

2.-. 181818 + 3.-1.54545 = -5 

-5--. 181818 - 2«- 1.54545 = 4 

For the second equation: 

TO SOLVE SIMULTANEOUS EQUATIONS FOR TWO UN- 
KNOWNS, 

REARRANGE THE EQUATIONS INTO THE FOLLOWING 
FORMAT. 

AX + BY = C 
DX + EY = F 

ENTER THE COEFFICIENTS - A,B,C,D,E,F, 
? 9.5, -8.5,3,2,3,-5 
X = -.554546 AND Y = -.972727 
As proof: 9.5'-. 554546 - 8* - .972727 = 3 
2»-. 554546 + 4«-. 972727 = -5 
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NUMERICAL INTEGRATION 

This program for numerical integration provides an easy and 
convenient way to solve the common integral calculus problem. 
One especially useful case is to solve for the area under a curve. In 
fact, by solving for the area bounded by the curve which describes 
an antenna directivity pattern, the gain of that antenna can be 
calculated. 

In this program you must use an odd number of points, and the 
points must be equally spaced. If data is not available for equally 
spaced points, try using one of the curve fitting programs to 
determine the equation of the line bounding the area. From this 
equation, equally spaced points can be selected and solved to 
provide data for the numerical integration. 

The actual integration is solved by use of Simpson's Rule 
which is: 

J x " S 

x o f(x) dx = y 

[f(x ) + 4f(x,) + 2f(x 2 ) . . . 2f(x n _ 2 ) + 4f(x n _ 1 ) + f(x„)] 

Numerical Integration Program 

5 'NUMERICAL INTEGRATION 

10 CLS 

15 DIMV*<100> 

20 PR I NT "ENTER THE SPACING BETWEEN 
THE SELECTED POINTS. AN ODD"CH 
R$( 13) "NUMBER OF POINTS MUST BE 
USED. " 

25 INPUT S".S=S/3 

30 PR I NT "ENTER THE VALUE FOR A SEL 
ECTED POINT. USE A / TO DENOTE 
"CHR*(13)"THE END OF DATA." 

40 FOR 1=1 TO 100 

50 INPUT V$(I> 

55 IF V$(I)="/"GOTO80 

60 V(I)=VAL(V*(I>) 

65 0=0+1 

70 NEXT I 

80 NI=S*V<1)+S*V(Q> 

90 FOR 1=2 TO Q-1STEP2 

100 A=A+VU)*S*4 
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110 NEXT I 

120 FOR 1=3 TO 0-2 STEP 

130 B=B+V(I)*S*2 

140 NEXT I 

150 NI=NI+A+B 

160 PRINT" ANSWER — "NI 
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Fig. 1-2. The area (under the curve) which can be solved by using the 1st 
sample program for the numerical integration program. 

Numerical Integration Sample Problems 

This first problem (see Fig. 1-2) will solve for the area under 
the indicated curve. The greater the number of points used, the 
more accurate the answer. 

ENTER THE SPACING BETWEEN THE SELECTED POINTS. 
AN ODD NUMBER OF POINTS MUST BE USED. 
? .314159 

ENTER THE VALUE FOR A SELECTED POINT. USE A / TO 
DENOTE THE END OF DATA. 
?0 

? .0955 
? .3455 
? .6545 
? .9045 
? 1.0000 
? .9045 
P.6545 
?.3455 
? .0955 
?0 
?/ 
ANSWER -1.5708 

If the program were set up for double precision, the actual 
answer would be 1.570794939994812. To see the difference in 
answer accuracy, the next example will be run using seven instead 
of eleven points. 

ENTER THE SPACING BETWEEN THE SELECTED POINTS. 
AN ODD NUMBER OF POINTS MUST BE USED. 
? .523598 
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ENTER THE VALUE FOR A SELECTED POINT. USE A / TO 

DENOTE THE END OF DATA. 

?0 

?.25 

P.6545 

PI 

P.6545 

?.25 

?0 

ANSWER - 1.5704 

With double precision, the answer is 1.570412255525589. 
The actual difference in answers is very small, but it is used to 
illustrate that the number of points will influence the final answer. 

And finally, since the area under the curve is really a semicir- 
cle, the actual answer is 1.570796327. 
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Chapter 2 
Finance 



The following chapter focuses on financial programs and deals with 
three main areas: loans, savings, and compounded interest. In each 
case, we're generally concerned with what any payment might be, 
how much time or number of payments required, and total value. 

The most basic programs deal with compounded amounts/ 
interest. The main variables in the different equations are: initial 
deposit (A), interest rate (I), future value (V), number of years (T) 
and number of times the principal is compounded per year (TY). 
Equations for the following programs are: 

V= A*(l + I/TY)t(TY«T) 

T= log(V)/log(l + I/TY)"TY 
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COMPOUND INTEREST 

This program will compute the final amount after some initial 
principal has been compounded periodically, with no further 
deposits. Additional deposit programs will be shown later in this 
chapter. 

Compound Interest Program 

5 'COMPOUND INTEREST 

10 CLS 

20 INPUT "ENTER AMOUNT TO BE COMPO 

UNDED" 5 A 
30 INPUT "ENTER INTEREST RATE" 5 1 : 1 

=1/100 
40 INPUT "ENTER NUMBER OF YEARS MO 

NEY IS COMPOUNDED" ;T 
50 INPUT "ENTER NUMBER OF COMPOUND 

PERIODS PER YEAR" ;TY 
62 V=A*<1 + I/TY)MTY*T> 
70 PRINT"*" INT<V*100+.5)/100 
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Sample Compound Interest Problems 

In these problems, the interest rate is entered as a percentage 
number rather than a decimal, i. e. , 6 percent instead of 0. 06. What 
will be the value of $1500 in 9 years if it is deposited at 7. 5 percent 
(year interest rate) and compounded quarterly? 

ENTER AMOUNT TO BE COMPOUNDED 

?1500 

ENTERINTERESTRATE 

?7.5 

ENTERNUMBER OF YEARS MONEYIS COMPOUNDED 

?9 

ENTER NUMBER OF COMPOUND PERIODS PER YEAR 

?4 

$2927.69 

After $15,000 is deposited for 20 years at 5.5 percent yearly 
interest (compounded semi-annually), what will be the total value 
of the principal? 

ENTER AMOUNT TO BE COMPOUNDED 

? $15,000 

ENTERINTERESTRATE 

?5.5 

ENTER NUMBER OF YEARS MONEYIS COMPOUNDED 

?20 

ENTER NUMBER OF COMPOUND PERIODS PER YEAR 

?2 

$44,298 
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YEARS OF INTEREST 

This program deals with the number of years some principal 
will have to remain on deposit to reach a specific value. 

Years of Interest Program 

5 'NO. YEARS OF COMPOUND INTEREST 

10 CLS 

20 INPUT "ENTER AMOUNT TO BE COMPO 

UNDED" ;A 
30 INPUT "ENTER INTEREST RATE";I:I 

=1/100 
40 INPUT "ENTER FINAL DESIRED AMOU 

NT" ;V 
50 INPUT "ENTER NUMBER OF COMPOUND 

PERIODS PER YEAR" ;TY 
60 T=LOG(V/A)/< (LOG(l+I/TY) )*TY) 
70 PRINT" COMPOUNDED" TY" TIMES PER Y 

EAR WILL REOUIRE THE MONEY TO B 

E DEPOSITED" INT<T*10+. 5 ) /10" YEA 

RS" 
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Sample Problems 

How many years will $1000 dollars have to remain on deposit 
if the interest rate is 9.5 percent, the money is compounded 
quarterly, and the final desired amount is $2000? 

ENTER AMOUNTTO BE COMPOUNDED 

? 1000 

ENTER INTEREST RATE 

?9.5 

ENTER FINAL DESIRED AMOUNT 

?2000 

ENTERNUMBER OF COMPOUND PERIODSPER YEAR 

?4 

COMPOUNDED 4 TIMES PER YEAR WILL REQUIRE 

THE MONEY TO BE DEPOSITED 7.4 YEARS. 

If the annual inflation rate averages 10 percent, in how many 
years will the price of some item double? For this problem assume 
that the price is currently $100 and will double to $200. 
ENTER THE AMOUNTTOBE COMPOUNDED 
?100 
ENTERINTERESTRATE 

?10 

ENTER FINAL DEPOSITED AMOUNT 

?200 

ENTER NUMBER OF COMPOUND PERIODS PER YEAR 

?1 

COMPOUNDED 1 TIME PER YEAR WILL REQUIRE THE 

MONEY TO BE DEPOSITED 7.3 YEARS. 

Or, in other words, with an annual inflation rate of 10 percent, 
prices will double every 7. 3 years! 
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INTEREST RATE 

With the many different types of savings instruments avail- 
able, a person often has to evaluate each instrument versus his 
particular needs. This program is designed to calculate the amount 
of interest required to attain some final amount of money. This is 
especially useful with the varying amounts of time required by 
different savings plans. Note that the interest rate can also be 
considered as a rate of return for some initial deposit or 
investment. 

Interest Rate Program 

5 'RATE OF RETURN 

10 CLS 

20 INPUT "ENTER AMOUNT TO BE COMPO 

UNDED" ?A 
30 INPUT "ENTER FINAL AMOUNT" 5 V 
40 INPUT "ENTER NUMBER OF YEARS MO 

NEY IS COMPOUNDED" ?T 
50 INPUT "ENTER NUMBER OF COMPOUND 

PERIODS PER YEAR" ?TY 
60 I=( <V/AH <1/(TY*T))-1)*TY 
70 PRINT" RATE OF RETURN IS" 1*100" P 

ERCENT PER YEAR" 
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Sample Interest Rate Problems 

What interest rate is required to double your money in 2 
years? 

ENTER AMOUNT TO BE COMPOUNDED 
?1 

ENTER FINAL AMOUNT 
?2 

ENTER NUMBER OF YEARS MONEY IS COMPOUNDED 
?2 

ENTER NUMBER OF COMPOUND PERIODS PER YEAR 
?1 
RATE OF RETURN IS 41.4214 PERCENT PER YEAR 

Obviously, doubling your money requires a high rate of 
return, a rate that is not available through normal savings interest. 
Therefore, to double your money in a short period of time requires 
a more speculative type of investment, one that causes the chances 
of losing your money to increase in proportion to the rate of return. 

You wish to make an initial deposit of $5000 into a college 
education fund which will yield $25000 in 15 years. If the money is 
compounded quarterly, what rate of return must the fund yield per 
year? 

ENTER AMOUNT TO BE COMPOUNDED 

?5000 

ENTER FINAL AMOUNT 

? 25000 

ENTER NUMBER OF YEARS MONEYIS COMPOUNDED 

?15 

ENTER NUMBER OF COMPOUND PERIODS PER YEAR 

?4 

RATE OF RETURN IS 10.8748 PERCENT PER YEAR 
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FUTUREVALUE 

This group of programs deals with the same basic savings 
theme, though with periodic payments made into some savings 
instrument. First, the most basic: what is the future value of X 
dollars deposited Y times at Z interest? The equation for this 
computation is: 

Future Value Program 

FV = -j-((l+I) T+1 -(l+D) 

5 'FUTURE VALUE 

10 CLS 

20 INPUT "ENTER AMOUNT OF MONTHLY 

DEPOSIT" ;A 
30 INPUT "ENTER INTEREST RATE";I:I 

=1/1200 
40 INPUT "ENTER NUMBER OF YEARS" 5T 

:T=T*12 
50 FV=A/I*< (l+I)t(T+l)-(l+I) ) 
60 PRINT"$"INT(FV*100+.5>/100 
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Sample Future Value Problems 

How much money will a person have if he deposits $200 each 
month for 10 years? The interest rate is 6 percent per year. 
ENTER AMOUNT OF MONTHLY DEPOSIT 
?200 

ENTER INTEREST RATE 
?6 

ENTER NUMBER OF YEARS 
?10 
$32939.60 

Note that the actual money deposited was $24000 ($2400/ 
year for 10 years). 
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REQUIRED PAYMENT 

The next program deals with the required deposit each month 
to attain some specific goal. This is extremely useful when trying 
to determine how much money is available from a budget for the 
purchase of some article, i.e., a car payment. To determine the 
monthly payment, use the following equation: 

A" (1+I)T + 1_ (1+I) 

Required Payment Program 

5 'PAYMENT FOR AMOUNT 

10 CLS 

20 INPUT "ENTER FINAL DESIRED AMOU 

NT" ;FV 
30 INPUT "ENTER INTEREST RATE" Si: I 

=1/1200 
40 INPUT "ENTER NUMBER OF YEARS" 5 T 

sT=T*12 
50 A=(FV*I)/<<1STEPRETURNS0NLEFT$ 
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Sample Payment Problems 

You will need $15000 in 7 years. How large a monthly 
payment do you need if the annual interest rate is 8. 5 percent? 

ENTER FINAL DESIRED AMOUNT 

? 15000 

ENTER INTEREST RATE 

?8.5 

ENTER NUMBER OF YEARS 

?7 

$ 130.38 MONTHLY PAYMENTS 

A new car costs $7000. What is the monthly payment over 4 
years if the interest rate is 12.5 percent? Disregard any interest 
payments. 

ENTER FINAL DESIRED AMOUNT 

?7000 

ENTER INTEREST RATE 

112.5 

ENTER NUMBER OF YEARS 

?4 

$ 111.98 MONTHLYPAYMENTS 

This amount is necessary to pay off the principal in 4 years. As 
you'll see later in this chapter, the amount required for interest 
payments raises the monthly payment considerably. 
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NUMBER OF PAYMENTS I: SAVINGS 

The final program in this group of three (Future Value and 
Required Payment are the other two) will determine the number of 
equal payments necessary to achieve some final amount. As 
before, you will enter the monthly deposit, the total desired 
amount, and the interest rate. 

To solve for the number of payments, use the following 
equation: 





K 


Io 8(V^i +1 + I) 










log(l + D 




Number of Payments Program 








5 ; 

10 

20 

30 

40 


' NUMBER 

INPUT « 
ED" 5P 
INPUT " 
ERCENT" 
INPUT " 


OF PAYMENTS 

'ENTER AMOUNT TO 

'ENTER INTEREST 1 
■*I 11=1/1200 
'ENTER AMOUNT OF 


BE FINANC 
RATE IN P 
AFFORDAB 



LE MONTHLY PAYMENT" ; A 
50 N=LOG(A/(A-P*I))/LOG(l+I) 
60 PR INTN" MONTHLY PAYMENTS" 
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Sample Payments Program 

How long will it take to save $1200 if you are making monthly 
payments of $15 at 10 percent interest? 

ENTER AMOUNT OF MONTHLY PAYMENT 

?15 

ENTER INTEREST RATE 

?10 

ENTER FINAL AMOUNT 

? 1200 

$ 1200 WILL REQUIRE DEPOSITS OF $15 FOR 61.1569 

MONTHS 

You make quarterly deposits of $450 into a savings account. If 
the interest rate is 7.5 percent, how many months will it take to 
accumulate $15000 in the account? 

ENTER AMOUNT OF MONTHLY DEPOSIT 

? 150 (450/3) 

ENTER INTEREST RATE 

?7.5 

ENTER FINAL AMOUNT 

? 15000 

$15000 WILL REQUIRE DEPOSITS OF $150 FOR 77.5398 

MONTHS 

(Or, quarterly deposits for 6. 5 years) 
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MONTHLY LOAN PAYMENT 



This program will calculate the total monthly loan payment 
given the amount, interest rate, and number of payments. The 
program is based on: 



A = P[ 



1/12 



1- (1-1/12) " 



■] 



Loan Payment Program 

5 'MONTHLY LOAN PAYMENT 

10 CLS 

20 INPUT "ENTER AMOUNT TO BE FINAN 

CED" ;P 
30 INPUT "ENTER INTEREST RATE IN P 

ERCENT";I:I=I/100 
40 INPUT "ENTER NUMBER OF YEARS OF 

FINANCING" 5N:N=N*12 
50 A=P*<I/12)/<l-<i + <I/12)M-N> 
60 PRINT" MONTHLY PAYMENT IS $" INT( 

A*100+.5)/100 
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Loan Payment Program Examples 

A new car costs $7000. What is the monthly payment over 4 
years if the interest rate is 12. 5 percent? 
ENTER AMOUNTTO BE FINANCED 
?7000 

ENTER INTEREST RATE IN PERCENT 
? 12.5 

ENTER NUMBER OF YEARS OF FINANCING 
?4 
MONTHLY PAYMENT IS $186.06 

You'll note that this is considerably above the $112 payment. 
Payment for the interest amounts to about two-thirds of the 
principal payment. 
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LOAN BALANCE 

This section deals with the balance on a loan financed over 
some period of time. As demonstrated in the previous example 
(Monthly Loan Payment), payments on the interest are a signifi- 
cant portion of a monthly payment. In fact, in the case of a home 
mortgage, the interest payment is virtually the entire payment for 
the first several years. This Loan Payment program is designed to 
show the principal balance and accumulated interest after a desired 
number of payments. With this information, you can generate a 
schedule of interest paid, a number that is used when computing 
your Federal Income Tax return. 

To compute the balance and interest, use: 

E _ i fA , (i+ir-i rl 



AI = B-B ,+(y-X + l)A 

y x-1 w ' 

Loan Balance Program 

5 'LOAN BALANCE 

10 CLS 

20 INPUT "ENTER AMOUNT TO BE FINAN 

CED" ;P 
30 INPUT "ENTER INTEREST RATE IN P 

ERCENT";I:I=I/1200 
40 INPUT "ENTER MONTHLY PAYMENT" ;A 
50 INPUT "ENTER BEGINNING PAYMENT 

NUMBER" ;X 
60 INPUT "ENTER ENDING PAYMENT NUM 

BER" ;Y 
70 B=i/((I+l)t-Y)#((A*(((l+Ilt-Y)- 

1))/I+P) 
80 PR I NT "BALANCE AT THE END OF PAY 

MENT NUMBER" Y"— $" INT<B*100+. 5 

)/100 
90 AI=B-(l/( (I+l)t-(X-l))*((A*(((l 

+ IH-<X-i))-m/I+P)) + <Y-X+l)*A 
100 PRINT "ACCUMULATED INTEREST AT 
THE END OF PAYMENT NUMBER "Y" — 
*"INT(AI*100+.5>/100 
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Sample Loan Balance Problems 

A $40000 mortgage is arranged such that the first payment is 
in January. The mortgage is 9 percent over 30 years and the 
payment is $321.85 (use the previous program for Monthly Loan 
Payment to compute the payment). Determine the interest paid 
after the first, fifth, twentieth, and thirtieth year. 

ENTER AMOUNT TO BE FINANCED 

? 40000 

ENTER INTEREST RATE IN PERCENT 

?9 

ENTER MONTHLY PAYMENT 

? 321.85 

ENTER BEGINNING PAYMENT NUMBER 

?1 

ENTER ENDING PAYMENT NUMBER 

? 12 (12* month) 

BALANCE AT THE END OF PAYMENT NUMBER 12 - 

$39726.70 

ACCUMULATED INTEREST AT THE END OF PAYMENT 

NUMBER 12 -$3588.92 

ENTER AMOUNT TO BE FINANCED 

? 40000 

ENTERINTERESTRATEIN PERCENT 

?9 

ENTER MONTHLYPAYMENT 

? 321.85 

ENTER BEGINNING PAYMENT NUMBER 

?49 

ENTER ENDING PAYMENT NUMBER 

?60 

BALANCE AT THE END OF PAYMENT NUMBER 60 - $38352 

ACCUMULATED INTEREST AT THE END OF PAYMENT 

NUMBER 60- $3471.02 

The total amount paid per year on the mortgage is $ 3862.20. 
So even at the end of the fifth year, almost 90 percent of the 
monthly payment is still going for interest payments. However, 
look at the situation in the fifteenth and twentieth years: 
ENTER AMOUNT TO BE FINANCED 
? 40000 

ENTER INTEREST RATE IN PERCENT 
?9 
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ENTER MONTHLY PAYMENT 

? 321.85 

ENTER BEGINNING PAYMENT NUMBER 

?169 

ENTER ENDING PAYMENT NUMBER 

?180 

BALANCE AT THE END OF PAYMENT NUMBER 180 - 

$31732.10 

ACCUMULATED INTEREST AT THE END OF PAYMENT 
NUMBER 180 -$2903.31 

ENTER AMOUNT TO BE FINANCED 

? 40000 

ENTER INTEREST RATE IN PERCENT 

?9 

ENTER MONTHLY PAYMENT 

? 321.85 

ENTER BEGINNING PAYMENT NUMBER 

?229 

ENTER ENDING PAYMENT NUMBER 

?240 

BALANCE AT THE END OF PAYMENT NUMBER 240 - 

$25407.30 

ACCUMULATED INTEREST AT THE END OF PAYMENT 

NUMBER 240- $2360.87 

Even at the end of the twentieth year, more than 60 percent of 
the total payment is still going for interest. 
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NUMBER OF PAYMENTS II: LOAN 

How long will it take to pay something off with a known 
monthly payment? This is not an unusual question in today's 
finances. To calculate the number of payments, use the following 
equation: 



log (1 —*££-) 



N = - 



log(l + I) 
Number of Payments Program 



5 ' 


■NO. OF PAYMENTS 


10 


CLS 


20 


INPUT "ENTER AMOUNT OF MONTHLY 




DEPOSIT" ;A 


30 


INPUT "ENTER INTEREST RATE";Is 



I 

=1/1200 
40 INPUT "ENTER FINAL AMOUNT" 5 V 
50 N=<LOG(V*I/A+l+I)/<LOG(l+I))-l) 
60 PRINT" $"V" WILL REQUIRE DEPOSITS 
OF *" A" FOR" N" MONTHS" 
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Number of Payments Examples 

If I can afford a monthly payment of $150 and the prevailing 
interest rate is about 12 percent, how many payments will it take to 
pay off a loan of $5000? 

ENTER AMOUNT TO BE FINANCED 

?5000 

ENTER INTEREST RATE IN PERCENT 

?12 

ENTER AMOUNT OF AFFORDABLE MONTHLY PAYMENT 

?150 

40.7 MONTHLY PAYMENTS 
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LOAN AMOUNT 



Conversely, instead of how many months it will take to pay 
something off, let's try: "How much can I afford to spend if I commit 
so much per month to another payment?" The equation for this 
program is: 



, =A . (i-q+p-"\ 



Loan Amount Programs 

5 'LOAN AMOUNT 

10 CLS 

20 INPUT "ENTER NUMBER OF PAYMENTS 

DESIRED" ?N 
30 INPUT "ENTER INTEREST RATE IN P 

ERCENT" Si: 1=1/1200 
40 INPUT "ENTER AMOUNT MONTHLY PAY 

MEIMT" ;A 
50 P=A*U-(1 + I)t-N)/I 
60 PRINT" YOU CAN AFFORD TO BORROW 

$"P 
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Sample Loan Amount Programs 

If I can afford to pay $100 per month for the next 15 months, 
how much new furniture can I purchase if the interest rate is 12 
percent? 

ENTER NUMBER OF PAYMENTS DESIRED 

?15 

ENTER INTEREST RATE IN PERCENT 

?12 

ENTER AMOUNT OF MONTHLY PAYMENT 

?100 

YOU CAN AFFORD TO BORROWS 1386.49 

Current car loans have an interest rate of 10.75 percent. If a 
new car loan runs for 4 years, how much can I afford if my monthly 
payment is $ 150? 

ENTER NUMBER OF PAYMENTS DESIRED 

?48 

ENTER INTEREST RATE IN PERCENT 

? 10.75 

ENTER AMOUNT OF MONTHLY PAYMENT 

?150 

YOU CAN AFFORD TO BORROW $5831.05 
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SINKING FUND 

The last program in this chapter is for a sinking fund type of 
payment. The sinking fund method of payment is based on ac- 
cumulating money through periodic deposits which also earn some 
interest. To calculate the monthly deposit use the formula: 



A = 



P»I 



(1+ 1) N - 1 



Sinking Fund Program 



5 ' 

10 

20 

30 

40 

50 
60 



SINKING 

CLS 

INPUT 

ULATED 

INPUT 

ERCENT 

INPUT 



FUND 



"ENTER AMOUNT TO BE ACCUM 

" ?P 

"ENTER INTEREST RATE IN P 

" 51:1=1/100 

"ENTER NUMBER OF YEARS OF 
ACCUMULATION" ;N 
A=P*I/<(<i + I)+N>-l> 
PRINT"YEARLY DEPOSIT IS *"INT< 
A*100+.5>/100 
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Sinking Fund Examples 

How much money must be deposited at the end of each year in 
a 5-year sinking fund to accumulate $5000 if the annual interest rate 
is 7. 5 percent? 

ENTER AMOUNT TO BE ACCUMULATED 

?5000 

ENTER INTEREST RATE IN PERCENT 

?7.5 

ENTER NUMBER OF YEARS OF ACCUMULATION 

?5 

YEARLYDEPOSITIS$860.82 
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Chapter 3 
Statistics 



This chapter on statistics will deal with several different areas, not 
just the common mean standard deviation. Included are programs 
on curve fitting, probability, and number evaluation. Some of the 
problems will find use within other even larger problems. And, as 
mentioned in Chapter 2, some of the programs can be used to 
supply data for other areas. 
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PERMUTATIONS 

The first program in this chapter deals with the operation 
known as permutations. Or, how many different ordered groups can 
be obtained from some set of objects without repetition. As an 
example, how many 3-letter permutation groups are there in the 
letters W.X.Y, and Z? 

To solve the permutation problem requires the use of the 
factorial operation. To explain, 3 factorial, represented as 3!, is the 
product of the form 3-2-1 = 6. The 3 denotes the highest integer in 
the series. On the other hand, 5! would equal 120 (5-4-3«2-l = 
120). 

The general equation for a permutation is: 

m p n= _JL 



(M-N)! 



Therefore, the previous question of how many 3-letter 
permutation groups are there in the letters W, X, Y, and Z could be 
expressed as 4 P 3. 
Evaluating this expression produces 4 P 3 = 



(4-3)! 1! 



4'3»2'1 _24 i 



Permutation Program 

5 'PERMUTATIONS 
10 CLS 

20 PRINT" INPUT M" : INPUT M 
30 PRINT" INPUT N": INPUT N 
40 IF N=0 P=l:GOTO 110 
50 IF N=l P=N:GOTO 110 
55 IF M-N<2 THEN N=M 
60 FOR 1=2 TO N-l 
65 RU)=M-I + 1 

70 IF $=1 GOTO 80 ELSE P=M*R<I> 
75 Q=l sGOTO 90 
80 P=P*R(I) 
90 NEXT I 
100 P=P*(M-N+1) 
110 PRINTM5"P";N"="P 
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Sample Permutation Problems 

Evaluate 7 P 5 
INPUT M 
?7 

INPUT N 
?5 
7P5 = 2520 

Evaluate 3 P 2 
INPUT M 
?3 

INPUT N 
?2 
3P3 = 6 

In this example you should notice that the program has 
changed the problem from 3 P 2 to 3 P 3. This is a special case and the 
program changes N from 2 to 3. To check the answer 

-3! = 3«2'1 _ 6 
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COMBINATIONS 

The difference between combinations and permutations is 
that a combination is without regard to order. For example, in the 
group of letters X, Y, and Z, there are 6 permutations (XY, YX, XZ, 
ZX, YZ, ZY). With XY and YX there is only 1 combination, XY. So, 
for the letters X, Y, and Z, there are 3 combinations, XY, XZ, and 
YZ. 

The equation for combinations is: 



M C N = — - M! 



(M-N)!«N! 
Combination Program 



5 ' 


COMBINATIONS 




10 


CLS 




20 


PRINT" INPUT M" : INPUT 


M 


30 


PRINT" INPUT N" : INPUT 


N 


35 


IF M-N<N THEN N=M-N 




40 


IF N=0 C-l: GOTO 130 




50 


IF N=l OM: GOTO 130 




60 


R<1)=M-N+1 




70 


FOR 1=2 TO N-l 




80 


R(I)=(M-N+I)/I 




90 


IF 0-1 GOTO 100ELSE C=! 


95 


0=1: GOTO 110 




10(Z 


) C=C*R(I) 




11C 


) NEXT I 




12E 


) C=C*<M/N> 




132 


) PRINTM; 5" C" ;N"="C 





R<i)*RU> 
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Sample Combination Problems 

In how many ways can a group of 3 people be chosen from a 
larger group of 10 people? (10 c 3) 
INPUT M 
?10 

INPUT N 
?3 

10C3 = 120 

There are 120 different ways to choose the 3 people for the smaller 
group. 

A box contains 4 cards numbered 1 through 4. If 3 cards are 
drawn from the box, how many different combinations of numbers 
are possible? 

INPUT M 
?4 

INPUT N 
?3 
4C3 = 4 
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ARITHMETIC MEAN 

The arithmetic mean or common average of a set of numbers is 
the sum of those numbers divided by the quantity of numbers. For 
example, in the arithmetic mean, you would be asked to determine 
the mean of 2, 4, 3.1, 5, 6.7, 9, a total of 6 different numbers. The 
mean is the sum of these numbers divided by 6. Or, more 
specifically, 

II. i- ±u 

1 = 1 

where the numbers are represented as aj, a 2 , a 3> . . . a n . 
Arithmetic Mean Program 

5 'ARITHMETIC MEAN 

10 CLS 

20 PRINT" ENTER NUMBERS FOR THE ARI 

THMETIC MEAN" 
30 INPUT N* 
40 IF N$=" /" GOTO 90 
50 N=VAL(N$) 
60 A=A+N 
70 X=X+1 
80 GOTO 30 
90 PRINT"THE ARITHMETIC MEAN IS"A/ 

X 
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Sample Arithmetic Mean Problems 

Determine the mean of 2, 4,3. 3.1, 5, 6.7, 9 (note that in the 
use of this program a / is used to denote the end of the information). 

ENTER NUMBERS FOR THE ARITHMETIC MEAN 

?2 

?4 

?3.1 

?5 

?6.7 

?9 

?/ 

THE ARITHMETIC MEAN IS 4.9667 

Population figures list the following as the number of children 
from the different families in a small rural town, 2, 4, 1, 5, 3, 10, 0, 
6. What is the arithmetic mean for these families? 
ENTER NUMBERS FOR THE ARITHMETIC MEAN 
?2 
?4 
?1 
?5 
?3 
?10 
?0 
?6 
?/ 
THE ARITHMETIC MEAN IS 3.875 

You'll note that even through the seventh family did not have 
any children, they were counted in the total number of families (8). 
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GEOMETRIC MEAN 

The geometric mean is another type of average based on 
determining a specific root of a set of numbers. Using the same set 
of numbers 2, 4, 3. 1, 5, 6.7, 9, the geometric mean is determined 
from the sixth root of the product of the numbers. More 
specifically, 

A = »V V VV..n B 

Geometric Mean Program 

5 'GEOMETRIC MEAN 

10 CLS 

20 PRINT" eNTER NUMBERS FOR THE GEO 

METRIC MEAN" 
25 A=l 
30 INPUT N$ 

40 IF N$=OPEN/" GOTO 90 
50 N=VAL(N$) 
60 A=A*N 
70 X=X+1 
80 GOTO 30 
90 PRINT" THE GEOMETRIC MEAN 

1/X) 
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Geometric Mean Sample Problems 

Determine the geometric mean of the examples given in the 
Arithmetic Mean section — 2, 4, 3. 1, 5, 6.7, 9 and 2, 4, 1, 5, 3, 10, 
0,6. 

ENTER NUMBERS FOR THE GEOMETRIC MEAN 
?2 
?4 
?3.1 
?5 
?6.7 
?9 
?/ 
THE GEOMETRIC MEAN IS 4.42205 

Note the difference between the arithmetic and geometric 
means — 4.9667 vs. 4.42205. 
ENTER NUMBERS FOR THE GEOMETRIC MEAN 
?2 
?4 
?1 
?5 
?3 
?10 
?0 
?6 
?/ 
THE GEOMETRIC MEAN IS . 

This example was selected to illustrate that when trying to 
determine a geometric mean you must have some amount for each 
data sample. In other words, zero is not valid data. Referring back 
to the equation that solves for the geometric means, we see that a 
zero forces a zero product. Deleting the zero produces: 
ENTER NUMBERS FOR THE GEOMETRIC MEAN 
?2 
?4 
?1 
?5 
?3 
?10 
?6 
?/ 

THE GEOMETRIC MEAN IS 3. 5567 
The arithmetic mean would be 4.41857. 
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HARMONIC MEAN 

The last type of mean to be illustrated is the harmonic mean. 
Instead of being based on either a sum or product as in the 
arithmetic and geometric mean, the harmonic mean uses the 
reciprocals of the numbers in the data set. The formula for the 
harmonic mean is: 



M = 



(where numbers are a Jt a 2 , a 3 , * . a n ) 
Harmonic Mean Program 



5 ' 


'HARMONIC MEAN 


10 


CLS 


20 


PRINT" ENTER NUMBERS 




MOIMIC MEAN" 


30 


INPUT N$ 


40 


IF N*="/" GOTO 90 


50 


N=VAL(n$) 


60 


A=A+(1/N) 


70 


X=X + 1 


80 


GOTO 30 


90 


PRINT "THE HARMONIC 1 




/X) 



FOR THE HAR 



MEAN IS"1/(A 



63 



Sample Harmonic Mean Programs 

Using the numbers from the examples in the geometric mean 
section, determine the harmonic mean. 
ENTER NUMBERS FOR THE HARMONIC MEAN 
?2 
?4 
?3.1 
?5 
?6.7 
?9 
?/ 

THE HARMONIC MEAN IS 3. 91403 
ENTER NUMBERS FOR THE HARMONIC MEAN 
?2 
?4 
?1 
?5 
?3 
? 10 
?6 
?/ 
THE HARMONIC MEAN IS 2. 7451 

You'll again note that since reciprocals are being used, zero is 
invalid data since it would force a division by zero. 

For the original set of data, 2, 4, 3.1, 5, 6.7, and 9, the means 
are: 

arithmetic -4. 96670 

geometric — 4.42205 

harmonic — 3.91403 
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GROUPED DATA MEAN 

There is one other type of mean that will find use in other 
programs, the mean of grouped data. Instead of finding the mean of 
fifty individual numbers, the numbers are grouped into specific 
categories. For example: 

category number 

10-19 14 

20-29 14 

30-39 22 

50 
Mechanically, the mean of this grouped data is found by 
determining the midpoint of each group, multiplying this midpoint 
by the number in that group to produce a category product, 
summing the category products, and finally dividing this sum by the 
total number: 

category number midpoint category product 

10-19 14 14.5 203 

20-29 14 24.5 343 

30-39 22_ 34.5 759 
50 1305 

The mean of this grouped data is 1305/50 = 26.1. The actual 
formula is: 

a* 2mc 

Mean = 

n 

where n is the number of pieces of data and mc is the category 
product. 

Mean Grouped Data Program 

5 'MEAN GROUPED 

10 CLS 

20 PRINT"ENTER THE RANGE OF THE PA 

RTICULAR CLASS (HIGHEST. LOWEST) 
it 

25 INPUT N*(2)»N*(D- 

30 IF N$(2>="/" GOTO 100 

40 PRINT-ENTER THE NUMBER WITHIH T 

HIS CLASS" 
45 INPUT N*<3) 
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50 N1=VAL(N$(1)) :N2=VAL<N*<2>> sN3= 

VAL(N$<3>) 
60 M=< (N2-Nl)/2)+Nl 
70 A=<N3*M)+A 
80 NT=NT+N3 
90 GOTO 20 
100 PRINT" THE MEAN IS" A/NT 
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Grouped Data Examples 

The time between repairs of various city buses is: 

time number of buses 
0-7 34 

8-14 12 

15-21 56 

22-28 7 

Find the mean time between repairs (to indicate the end of data, 
enter /, / for the range) . 

ENTER THE RANGE OF THE PARTICULAR CLASS (HIGH- 
EST, LOWEST) 
?7,0 

ENTER THE NUMBER WITHIN THIS CLASS 
?34 

ENTER THE RANGE OF THE PARTICULAR CLASS (HIGH- 
EST, LOWEST) 
?14,8 

ENTER THE NUMBER WITHIN THIS CLASS 
?12 

ENTER THE RANGE OF THE PARTICULAR CLASS (HIGH- 
EST, LOWEST) 
? 21,15 

ENTER THE NUMBER WITHIN THIS CLASS 
?56 

ENTER THE RANGE OF THE PARTICULAR CLASS (HIGH- 
EST, LOWEST) 
? 28,22 

ENTER THE NUMBER WITHIN THIS CLASS 
?7 

ENTER THE RANGE OF THE PARTICULAR CLASS (HIGH- 
EST, LOWEST) 

THE MEAN IS 13. 156 

Therefore, for the 109 buses in the city, each bus averaged 13. 156 

days between needed repairs. 

Another example is one that is often seen in evaluations. 
Students in an Introduction to Drama course were asked to grade 
their instructor on a scale of 1 to 5 as to the value they received 
from that course, with 1 indicating little or no value and 5 indicating 
tremendous value. The evaluations produced the following infor- 
mation: 



67 



Score - 12 3 4 5 
Number -15 8 9 3 21 

What is the mean rating on this course? 

Since the scores represent the actual midpoints of five 
infinitely small categories, the range information can be arbitrarily 
entered as plus or minus 1 from the midpoint. In other words, for 
the score 1, enter 2.0. 

ENTER THE RANGE OF THE PARTICULAR CLASS (HIGH- 
EST, LOWEST) 
?2,0 

ENTER THE NUMBER WITHIN THIS CLASS 
?15 

ENTER THE RANGE OF THE PARTICULAR CLASS (HIGH- 
EST, LOWEST) 
?3,1 

ENTER THE NUMBER WITHIN THIS CLASS 
?8 

ENTER THE RANGE OF THE PARTICULAR CLASS (HIGH- 
EST, LOWEST) 
?4,2 

ENTER THE NUMBER WITHIN THIS CLASS 
?9 

ENTER THE RANGE OF THE PARTICULAR CLASS (HIGH- 
EST, LOWEST) 
?5,3 
ENTER THE NUMBER WITHIN THIS CLASS 

?3 

ENTER THE RANGE OF THE PARTICULAR CLASS (HIGH- 
EST, LOWEST) 
?6,4 

ENTER THE NUMBER WITHIN THIS CLASS 
?21 

ENTER THE RANGE OF THE PARTICULAR CLASS (HIGH- 
EST, LOWEST) 

?/,/ 

THE MEAN IS 3. 125 
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STANDARD DEVIATION (UNGROUPED DATA) 

The various types of means have been used to determine the 
center of specific groups of data. The next most commonly 
requested piece of information is the amount of variability or the 
spread of the data. This standard deviation is based on the amount 
of variance from the mean of the data. Standard deviation is 
measured in the same units as the scores, whereas variance, the 
basis for standard deviation, is measured in the scores squared, 
i.e., s 2 is the symbol for variance and s is the symbol for standard 
deviation. The recognized formula for standard deviation is: 



s _ / 2 (x - x) 2 



n-1 

Standard Deviation Program 

5 'STANDARD DEV 

10 CLS 

20 PRINT" ENTER NUMBERS FOR THE STA 

NDARD DEVIATION (UNGROUPED DATA 

)" 
30 INPUT N* 
40 IF N$="/" GOTO 100 
50 N=VAL(N$) 
60 A=A+N 
70 X=NC2+X 
80 Y=Y+1 
90 GOTO 30 

100 PRINT" THE STANDARD DEVIATION I 
S"SOR( <X-A*2/Y)/<Y-1>) 
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Sample Standard Deviation Problems 

Find the standard deviation for the data 2, 4,3.1, 5, 6. 7, 9. 
ENTER NUMBERS FOR THE STANDARD DEVIATION (UN- 
GROUPED DATA) 
?2 
?4 
?3.1 
?5 
?6.7 
?9 
?/ 
THE STANDARD DEVIATION IS 2. 54925 

Find the standard deviation for the data 2, 4, 1, 5, 3, 10, 0, 6. 
ENTER NUMBERS FOR THE STANDARD DEVIATION (UN- 
GROUPED DATA) 
?2 
?4 
?1 
?5 
?3 
?10 
?0 
?6 
?/ 
THE STANDARD DEVIATION IS 3. 18198 
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STANDARD DEVIATION (GROUPED DATA) 

Standard deviation for grouped data is determined in much the 
same way as the mean of grouped data and standard deviation. To 
explain, as in the mean of grouped data, a midpoint of each class is 
determined and then used in the standard deviation formula. This 
formula is: 



S 



/sf(x-x) a 
V n-1 



Standard Deviation (Grouped Data) Program 

5 'STANDARD DEV. GROUPED 
10 CLS 

20 PRINT" ENTER THE RANGE OF THE PA 
RTICULAR CLASS (HIGHEST, LOWEST) 



25 


INPUT N*<2)>N*U) 


30 


IF N$<2>="/# IF110 


40 


PRINT"ENTGR THE NUMBER WITHIN T 




HIS CLASS" 


45 


INPUT N*<3) 


50 


N 1 =VAL ( N* ( 1 ) ) : N2=VAL < N$ < 2 )) 




N3=VAL(N*<3) ) 


60 


M=< (N2-Nl)/2)+Nl 


70 


A=(N3*M)+A 


80 


B= ( N3*M*M ) +B 


90 


NT=NT+N3 


100 GOTO 20 


110 PR I NT "THE STANDARD DEVIATION I 




S"LOG( <B-A»2/NT)/<NT-1) ) 
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Sample Problemsf or Standard Deviation (Grouped Data) 

Again, the time between repairs of various city buses is: 

time number of buses 
0-7 34 

8-14 12 

15-21 56 

22-28 7 

Determine the standard deviation. 

ENTER THE RANGE OF THE PARTICULAR CLASS (HIGH- 
EST, LOWEST) 
?7,0 
ENTER THE NUMBER WITHIN THIS CLASS 

?34 

ENTER THE RANGE OF THE PARTICULAR CLASS (HIGH- 
EST, LOWEST) 
?14,8 

ENTER THE NUMBER WITHIN THIS CLASS 
?12 

ENTER THE RANGE OF THE PARTICULAR CLASS (HIGH- 
EST, LOWEST) 
?21,15 

ENTER THE NUMBER WITHIN THIS CLASS 
?56 

ENTER THE RANGE OF THE PARTICULAR CLASS (HIGH- 
EST, LOWEST) 
? 28,22 

ENTER THE NUMBER WITHIN THIS CLASS 
?7 

ENTER THE RANGE OF THE PARTICULAR CLASS (HIGH- 
EST, LOWEST) 

?/,/ 

THE STANDARD DEVIATION IS 7. 1504 
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LINEAR REGRESSION 

The next three programs deal with the investigation of the 
relationship between two variables. Considering the two variables 
X and Y, what mathematical formula best describes the relation- 
ship between the given points? This first program will deal with 
linear regression, or a set of data points or a curve that 
approximates a straight line. The general equation for a straight 
line is Y = MX + B. From the data points it will be necessary to 
solve for the constants M and B which give the closest agreement. 
In addition to providing the constants, the program also solves for 
the coefficient of determination, a measure that indicates how 
closely the equation fits the initial data. The coefficient of 
determination is between and 1, and the closer it is to 1 the better 
the equation fits. 

Solving for the constants uses a method of least squares. The 
equations are: 



2 _ 


X V 

n 


2x 2 - 


(2) 2 

. x 7 

n 


T = y-SX 




where y = — - 




2 
x ~ n , 





R = 





Linear Regression Program 

5 'LINEAR REGRESSION 
10 CLS 

20 PRINT "ENTER THE GROUPS OF DATA 
POINTS <X»Y> 
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30 INPUT N*(1)jN*(2) 

40 IF N$<1>="/" GOTO 110 

50 N1=VAL(N$<1>) :N2=VAL<N$<2) ) 

60 A=N1*N2+A 

70 XT=N1+XT 

80 YT=N2+YT 

90 XS=N1*N1+XS 

95 YS=N2*N2+YS 

100 X=X+1 

105 GOTO 30 

110 S=<A-<XT*YT/X))/<XS-< (XT*2)/X) 

) 
120 T=YT/X-(S*<XT/X> ) 
130 R=(X*A-<XT*YT) )/<SQR(X*XS-(XT* 

2>>*<SQR<X#YS-(YT*2> ) )> 
140 PRINT" FOR THE GENERAL EQUATION 
Y=MX+B» M="S"AND B="T" . " CHR*< 

13) "THE COEFFICIENT OF DETERMI 

NATION IS" RE 2 
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Sample Linear Regression Problems 

You want to use a thermometer to measure the outside 
temperature at the top of a mountain. Unfortunately, the only 
thermometer that you can buy locally is measured in Centigrade. 
However, you do have an inside thermometer which is marked in 
Fahrenheit, but it only goes down to freezing, 32 degrees. How do 
you measure the outside temperature? 

One method might be to immerse both thermometers in a 
bowl of very hot water and record the temperatures of both 
thermometers as the water cools. From this information, you could 
determine the exact relationship between the two thermometers 
and calculate the outside temperature in Fahrenheit. 

The readings as the water cools are: 

Fahrenheit Centigrade 

180 82 

142 61 

133 56 

105 40 

90 32 

73 23 

64 18 

50 10 

ENTER THE GROUPS OF DATA POINTS (X, Y) 

? 82, 180 

? 61, 142 

? 56, 133 

? 40, 105 

?32,90 

? 23,73 

? 18,64 

? 10,50 

?/,/ 

FOR THE GENERAL EQUATION Y=MX+B, M=l 1.80864 

AND B= 31.8274. 

THE COEFFICIENT OF DETERMINATION IS . 999922 

This example was selected since the conversion equation 
between Centigrade and Fahrenheit is well known (F = 1.8C + 32) 
and readily fits the general equation for a straight line. Linear 
regression can be used in many ways to help correlate/extrapolate 
data. Consider Fig. 3-1 and find the equation that best describes 
the data. 
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ENTER THE GROUPS OF DATAPOINTS (X,Y) 
? 1.25,1.1 
? 3,3.4 
? 4.5,5.2 
? 6,7.2 

FOR THE GENERAL EQUATION Y=MX+B, M= 1.27648 AND 

B=-. 482017 

THE COEFFICIENT OF DETERMINATION IS . 999636 

So, from the four selected data points, the general equation of 
a line that best describes the data has been determined. From this 
line and equation, the data can be extrapolated out through values 
of 12 for X. 




Fig. 3-1. To be used in solving a linear regression problem. 
76 



X 


Y 


1 


.79 


2 


2.07 


3 


3.35 


8 


9.73 


9 


11.01 


10 


12.28 


11 


13.56 


12 


14.84 
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LOGRITHMIC CURVE FIT 

This next program is very similar to linear regression, except 
that, instead of being a straight-line type of relationship, it is a 
function of the natural log (e) of X. Fitting to a logrithmic curve 
uses the general equation Y = MlnX + B. The logrithmic curve is a 
typical form in that one axis of the curve is compressed in relation 
to the other. The following examples will clearly show this 
relationship. The equations used to solve for the logrithmic curve 
fit are: 



Zy In x (2 lnx2y) 

S = 



2 (In x) 2 - -1 (2 In x) 2 



n 



T=-L(2y- S2 lnx) 
n J 



R = 



(2 y lnx - — (2 lnxSlny)) 2 
(S(lnx) 2 - 1 ©lnx) 2 ) (2y 2 -l(Sy) 2 ) 



Logrithmic Curve Fit Program 



5 'LOG CURVE FIT 

10 CLS 

20 PRINT" ENTER THE GROUPS OF DATA 

POINTS (X»Y> 
30 INPUT N*(1)»N*<2> 
40 IF N$(l)="/" GOTO 130 
50 Ni=VAL<N*(l)>:N2=VAL<N*(2)> 
60 A=A+N2*L0G(N1> 
70 B=B+LOG(Nl) 
75 I=I+N2 

80 C=C+(L0G(N1))*2 
85 F=F+N2*2 
?0 G=G+LOG(Nl) 
100 D=B*2 
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110 X=X+1 

120 GOTO 30 

130 E=(A-B*I/X)/<C-D/X) 

140 H=(I-E*G)/X 

150 PRINT" IN THE GENERAL FORM OF Y 

=M(LN<X))+B» M="E"AND X="H"." 
160 J=< <A-B*I/X>*2>/< <C-D/X)*(F-< ( 

I*2)/X))) 
165 PRINT"THE COEFFICIENT OF DETER 

MINATION IS"J 
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Sample Logrithmic Curve Fit Problems 

Determine the equation of the line shown in Fig. 3-1. As can 
he seen, the curve can not he stretched to produce a straight line 
and the data points seem to compress as the values for X get 
bigger. Therefore, this is a good candidate for logrithmic curve fit 
program. Four data points were selected. 



X 


Y 


2 


.5 


4 


1.5 


6.3 


2.25 


9.5 


2.75 



If the coefficient of determination is not high enough, additional 

points can be selected and the program run again. 

ENTER THE GROUPS OF DATA POINTS (X, Y) 

?2,.5 

?4,1.5 

? 6.3,2.25 

? 9.5,2.75 

?/,/ 

IN THE GENERAL FORM OF Y+M(LN(X))+B f M= 1.46487 

ANDB=-. 510028. 

THE COEFFICIENT OF DETERMINATION IS .997925. 

Again, with four point, you have been able to determine the 
equation of the line that produces the best fit through the data 
points. More points could have been added, though it is a matter of 
diminishing returns as to whether the additional time is worthwhile 
for the small increase in accuracy. For example, adding the 
additional points 1.45,0 and 7.5,2.5 changes M to 1.48962, X to 
- .541302, and raises the coefficient of determination to . 998624, a 
very small change in accuracy rating. The actual equation of the 
line for this example was Y = 1 . 51nX - 0. 5. 
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EXPONENTIAL CURVE FIT 

The last curve fit program is for data that exhibits an 
exponential type of progression. You'll notice that the curve in Fig. 
3-3 is somewhat similar to that of Fig. 3-2 in that it seems to 
approach some limit. This is a case where it might be easiest to run 
the data set in both the exponential and logrithmic curve programs 
to determine the best fit. However, in some cases, inferences can 
be made because of the data range to help select the appropriate 
program to use. 

The equations to solve the exponential curve fit problem are: 



Sxlny- J- (Zx) (2/ny) 



£x 2 — L. (2x) 2 



/ Slny Sx 

T = EXPi— — -S 



n / 



(2xlny--l(2x21nY)) 2 

R= " 



lnY) 2 \ 
n / 



Exponential Curve Fit Program 



5 'EXP. CURVE FIT 

10 CLS 

20 PRINT"ENTER THE GROUPS OF DATA 

POINTS <X»Y) 
30 INPUT N*(1)»N*(2> 
40 IF N*<1) = V" GOTO 130 
50 Nl=VAL<N*(i>) :N2=VAL(N$(2) ) 
60 A=A+N1*L0G(N2) 
70 B=B+L0G(N2) 
75 I=I+N1 
80 C=C+N1*2 
85 D=D+N1 
90 F=DC2 
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100 G=G+(L0G(N2)*2> 

110 X=X+1 

120 GOTO 30 

130 E=<A-B*I/X)/(C-F/X) 

140 H=EXP<B/X-E*I/X) 

150 PRINT" iN THE GENERAL FORM OF Y 

=M*ECBXi M="E"AND X="H" ." 
160 J=< <A-B*I/X>»2)/< (C-F/X)*(G-( ( 

BC2)/X>>> 
165 PRINT"THE COEFFICIENT OF DETER 

MINATION IS"J 
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Sample Exponential Curve Fit Problem 

Using Fig. 3-3, solve for the equation that best describes the 
line in the data points. 



X 


Y 


1 


.9 


3 


.35 


5 


.13 


7 


.04 



ENTER THE GROUPS OF DATA POINTS (X, Y) 

?1,.9 

?3,.35 

?5,.13 

?7,.04 

?/,/ 

IN THE GENERAL FORM OF Y=M*E!BX, B=-. 516547 AND 

M= 1.58822. 

THE COEFFICIENT OF DETERMINATION IS .997247. 

Running the same data set in the logrithmic curve fit program 
produces: B= .881006 and M=-. 452093 and a coefficient of 
determination of .991973. Since the two coefficients of determina- 
tion are so close, it might be wise to use additional points to find the 
best fit. However, using more data points is not without its prob- 
lems. Since there are more points, the amount of error tends to 
reduce the coefficients of determination. 

For example, consider the following data set: 



X 


Y 


1 


.9 


3 


.35 


5 


.13 


7 


.04 


2 


.51 


6 


.75 


10 


.015 



When applied to the exponential curve fit program, B now equals 
-.421978, M = 1.48592, and the coefficient of determination is 
.71439. You'll notice that even though the changes in B and M are 
relatively small, the coefficient of determination has decreased 
significantly. With the increased number of data points, any error is 
now more significant. For comparison, when using the logrithmic 
curve fit program, B = .32053, M = .817324, and the coefficient of 
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determination is .536147. So, increasing the number of data points 
has lowered the confidence levels overall for B and M, yet they 
have provided greater assurance that the correct curve fit program 
is being used. 
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FACTORIAL 

As mentioned earlier in this chapter, a factorial, N, is the 
simple progression N • (N- 1) • (N-2) . . . (2) • (1). Integers must 
be used. 

Factorial Program 

5 'FACTORIAL 

10 CLS 

20 INPUT"ENTER AN INTEGER" 5N 

30 H=l 

40 FOR 1=1 TO N 

50 M=M*I 

60 NEXT I 

70 PRINT"THE FACTORIAL OF"N"IS"M 
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Sample Factorial Problems 

ENTER AN INTEGER 

?9 

THE FACTORIAL OF 9 IS 362880 

ENTER AN INTEGER 

?24 

THE FACTORIAL OF 24 IS 6.20448E+23 
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CHI-SQUARE EVALUATION 

The last program in this chapter deals with chi-square tests or 
evaluations. This evaluation can be used to determine two pieces 
of information: "goodness of fit," and where two variables are 
related or independent. Since this is not a book devoted to 
statistics, but one that shows how many different types of 
programs can be computerized for everyday use, I am not going to 
explain the complete procedure for using the chi-square evalua- 
tion. However, this program will solve for chi-square, which can 
then be used to determine goodness of fit and any relationship 
between variables. 

The chi-square evaluation calculates the relationship between 
the observed frequency and the expected frequency in the 
distribution of samples. The formula is: 

(O.-E.) 2 



Chi-Square Evaluation Program 



5 ' 


'CHI-SQUARE EVAL. 


10 


CLS 


20 


PR I NT "ENTER THE OBSERVED FREQUE 




NCY, EXPECTED FREQUENCY INFORMAT 




ION" 


25 


INPUT 0$»E* 


30 


IF 0$="/" GOTO70 


40 


0=VAL ( 0$ ) : E=VAL ( E* ) 


50 


KI=KI+( <0~E)*2)/E 


60 


GOTO 25 


70 


PR I NT "THE CHI-SQUARE EVALUATION 




IS"KI 
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Chi-Square Sample Problems 

Determine Chi-Square for the following information. 



served 


Expected 


20 


30 


60 


90 


no 


120 


50 


30 


60 


30 



ENTER THE OBSERVED FREQUENCY, EXPECTED FRE- 
QUENCY INFORMATION 
? 20,30 
? 60,90 
? 110,120 
? 50,30 
? 60,30 

THE CHI-SQUARE EVALUATION IS 57.5 

As I said, this information will now be used to determine 
whether the assumed hypothesis is true or false. 
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Chapter 4 
Electronics Programs 




This next chapter is devoted to programs encountered in many 
facets of electrical engineering. Admittedly, there are an infinite 
number of equations used in the field of electrical engineering, but 
these programs represent a fair cross section. There are also 
programs which might not ordinarily be at one's fingertips for 
instant use. In this chapter I will start going into detailed 
explanations of the programs themselves. So far in this book, the 
programs have generally not been sufficiently complicated to 
warrant a closer examination. However, in this chapter, a number 
of programs contain smaller programs which can be used either 
singly or in combinations to solve other problems. 
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ANTENNA PATTERN 

This first program of the chapter will be used to determine the 
radiation pattern of vertical antennas. In many cases, the radiation 
pattern can be found by checking any of the many reference texts 
available which show patterns for "common" antenna configura- 
tions. The problem arises when certain factors dictate antenna 
positioning that does not fall into a common configuration. At such a 
time this program proves invaluable. 

To compute the radiation pattern from a group of vertical 
radiators requires that you know the antenna spacing, phase 
difference, individual antenna currents, and positional location 
between the antennas. Figure 4-1 shows a diagram that will be 
used to explain the designations used in the program. In this 
diagram there are three antennas, with one designated as the 
reference antenna. The other two are located both in spacing and 
distance from this reference antenna. In addition to the physical 
positions, the second two antennas are related to the reference by 
their individual antenna currents and phase. The four designators 
for the antennas are: 

Antenna Current C 

Phase relationship to the reference P 

Spacing from the reference S 

Angular Displacement from the reference AD 

The last quantity, angular displacement, is used to orient the 
pattern in either a relative or true plane. To explain, if you are only 
concerned with the pattern relative to the antennas, use zero as the 
angular displacement for Antenna 1 and 60 degrees as the angular 
displacement for Antenna 2 (see Fig. 4-1). If you want to determine 
the pattern base on True North, use the actual true bearing from 
the reference to the selected antenna. Lines 50 through 150 of the 
Antenna Pattern program are used to input this data. The first input 
statement, line 20, sets the number of points at which the data will 
be plotted. In effect, it determines the number of times the pattern 
information will be summed. Line 50 inputs the total number of 
verticals in the antenna system. Note that this value, E, is 
decreased by 1 to take into account the reference element. The 
complete equation to determine the antenna pattern is: 

N-l 
A(0)= j + 2 A 1 ej* | 
n=l 
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Lines 160 through 270 of the program do the actual summing to 
determine the relative amplitude of the pattern for each plotting 
point. Note that in lines 200 and 210 the formula for resolving a 
polar form of notation to rectangular coordinates is used. The 
respective X and Y quantities are then summed in lines 220 and 230 
each time through the loop. Once the required number of loops 
have been completed, as determined by the number of elements, 
line 260 converts the X and Y information back to a polar notation. 
From here, line 270 returns back to line 170 to repeat this entire 
process for the next plotting point. 

Lines 275-430 show the display portion of the program. In this 
section, the information is formatted into three columns for 
presentation. Lines 271 through 300 display the information in a 
relative amplitude format. More useful than the relative amplitude 
format is one that presents the information after being normalized 
to the largest value. This program will present either normalized 
data or normalized data expressed in dB. Lines 310 and 320 search 
the information file to determine the largest value for use in the 
normalized data. 

Antenna Pattern Program 

5 7 ANTENNA PATTERN 
10 CLS 

20 INPUT "ENTER NUMBER OF DEGREES 
BETWEEN PLOTTING POINTS" ;D 



0° 



ANTENNA 1 




ANTENNA 2 



90° 



REFERENCE ANTENNA 



Fig, 4-1. Reference antenna with two other variably located antennas. 
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30 DIMN(360/D) 
35 PRINT 

45 PR I NT "THE REFERENCE ELEMENT IS 
DESIGNATED AS ELEMENT 0" 

49 PRINT 

50 INPUT "ENTER THE TOTAL NUMBER OF 
ELEMENTS" ;E 

55 PRINT 

60 E=E-1 

70 FOR 1=1 TO E 

80 PRINT"FOR ELEMENT" I" ENTER THE A 

NTENNA CURRENT. " 
90 INPUT C(I) 

100 PRINT"ENTER THE PHASE DIFFEREN 
CE BETWEEN ELEMENT" I" AND THE R 
EFERENCE " CHR* < 1 3 ) " ELEMENT . " 

101 PR I NT "USE + FOR LAGGING AND - 
FOR LEADING. " 

110 INPUTP(I) 

120 PRINT"ENTER THE SPACING IN DEG 
REES <ONE WAVELENGTH = 360 DEG 
REES > " CHR* (13)" BETWEEN ELEMENT 
'T'AND THE REFERENCE" 

130 INPUTSU) 

140 PRINT"ENTER THE ANGULAR DISPLA 
CEMENT IN DEGREES BETWEEN THE 
REFERENCE "CHR$< 13) "ELEMENT AND 
ELEMENT" I 

150 INPUT AD(I) 

160 NEXT I 

170 FOR J=0 TO 360STEP D 

180 FOR 1=1 TO E 

185 M=C0S((J-AD(I))/57.29578) 

190 G=S<I)*M-P<I) 

200 X=C<I)*C0S(G/57.29578) 

210 Y=C<I)*SIN(G/57. 29578) 

220 XS=XS+X 

230 YS=YS+Y 

240 NEXT I 

250 XS=XS+1 

255 F=J/D 

260 N(F)=S6»R(XSC2+YSC2) 
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261 XS=0 

262 YS=0 

270 NEXT J 

271 CLS 

275 PRINTTAB(9)"MA6. " T AB < 29 ) " MAG. " 

TAB (50)" MAG. " 
280 FOR 0=0 TO 110 STEP D 
295 PRINTQ+DTAB(5)"DEG"TAB(8)N< (Q+ 

D ) /D ) TAB ( 20 ) Q+D+ 1 20TAB ( 25 ) " DEC 

" TAB ( 28 ) N ( ( Q+D+ 1 20 ) /D ) TAB ( 4 1 ) 6 

+D+240TAB ( 46 ) " DEG " TAB ( 49 ) N < ( 0+ 

D+240)/D) 
300 NEXT Q 

310 FOR R=0 TO 360/ D 
320 IF N<R)>NG THEN NG=N(R):NEXT R 

ELSE NEXT R 
340 INPUT "DO YOU WANT NORMALIZED D 

ATA?" ;ND* 

345 IF ND$="Y"GOTO 347 

346 END 

347 CLS 

348 PRINTTAB(9)"MAG. " TAB ( 29) "MAG. " 
TAB (50) "MAG." 

350 FOR 0=0 TO 110 STEP D 

360 PRINTQ+DTAB(5)"DEG"TAB(8)N< (0+ 
D ) /D ) /NGTAB ( 20 ) Q+D+ 1 20TAB ( 25 ) " 
DEG"TAB(28)N( (Q+D+120) /D> /NGTA 
B ( 4 1 ) O+D+240TAB ( 46 ) " DEG" TAB ( 49 
)N( (Q+D+240)/D)/NG 

370 NEXT Q 

380 INPUT "DO YOU WANT NORMALIZED D 
ATA EXPRESSED IN DB?"?DB* 

390 IF DB*="Y"THENGOTO 400 

395 END 

400 CLS 

410 PRINTTAB ( 9) " DB" TAB ( 29 ) " DB" TAB ( 
50)"DB" 

415 FOR 0=0 TO 110 STEP D 

420 PRINTQ+DTAB ( 5 ) " DEG " TAB ( 8 ) 20* ( L 
OG(N( <O+D)/D)/NG)/LOG<10) )TAB( 
2 1 ) 0.+DH- 1 20TAB ( 27 ) " DEG " TAB ( 30 ) 2 
0* ( LOG ( N ( < O+D+120 ) /D ) /NG ) /LOG ( 
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10 ) ) TAB ( 43 ) Q+D+240TAB ( 47 ) " DEG" 
TAB (51) 20* ( LOG ( N ( ( Q+D+240 ) /D ) / 
NG)/LOG(10) ) 
430 NEXT Q 

Where 

A(0) = pattern amplitude as a function of 
A 1 = relative amplitude of the nth element 
N = total number of elements 
!f» n = j3d n cos(0 n --0)-a 
a = phase of the nth element relative to the 

reference element in degrees (+ for lagging, 
-for leading phase) 
/3d n = electrical distance (in degrees) from the 
driven element to the nth element 
= spatial angle between elements. 

Sample Antenna Pattern Problem 

For this example, use Fig. 4-1. The spacing between the 
Reference Antenna and Antenna 1 is 0.5 wavelengths (180 
degrees) and between the Reference and Antenna 2 is 1 
wavelength. Antenna 2's angular displacement is 60 degrees. The 
current in Antenna 1 is 1 and in Antenna 2 is 2. The phase between 
the reference and Antenna 1 is 90 leading and between the 
reference and Antenna 2 is 180 degrees leading. Determine the 
normalized amplitude in dB every 30 degrees. 

ENTER NUMBER OF DEGREES BETWEEN PLOTTING 
POINTS? 45 

THE REFERENCE ELEMENT IS DESIGNATED AS ELE- 
MENT 

ENTER THE TOTAL NUMBER OF ELEMENTS? 3 
FOR ELEMENT 1 ENTER THE ANTENNA CURRENT. 
?1 

ENTER THE PHASE DIFFERENCE BETWEEN ELEMENT 1 
AND THE REFERENCE ELEMENT. Use + FOR LAGGING 
AND- FOR LEADING. 
?_90 

ENTER THE SPACING IN DEGREES (ONE WAVELENGTH = 
360 DEGREES) BETWEEN ELEMENT 1 AND THE REFER- 
ENCE 
?180 



96 



ENTER THE ANGULAR DISPLACEMENT IN DEGREES BE- 
TWEEN THE REFERENCE ELEMENT AND ELEMENT 1 
?0 

FOR ELEMENT 2 ENTER THE ANTENNA CURRENT. 
?2 

ENTER THE PHASE DIFFERENCE BETWEEN ELEMENT 2 
AND THE REFERENCE ELEMENT. USE + FOR LAGGING 
AND- FOR LEADING. 
?-180 

ENTER THE SPACING IN DEGREES (ONE WAVELENGTH = 
360 DEGREES) 

BETWEEN ELEMENT2 AND THE REFERENCE 
?360 

ENTER THE ANGULAR DISPLACEMENT IN DEGREES BE- 
TWEEN THE REFERENCE ELEMENT AND ELEMENT 2 
'60 



MAG. 
30 DEG. 940224 
60DEG2 
90 DEG 2. 51373 
120 DEG 4 



MAG. 
150 DEG 1.08759 
180 DEG 3. 16228 
210 DEG 2.40561 
240DEG1.12352E-06 



MAG. 
270 DEG. 593411 
300 DEG 2 
330 DEG 1.67844 
360 DEG 3. 16228 



DO YOU WANT NORMALIZED DATA?? Y 



MAG. 
30 DEG. 235056 
60 DEG. 5 
90 DEG. 628432 
120 DEG 1 



MAG. 
150 DEG. 271897 
180 DEG. 790569 
210 DEG. 601402 
240DEG2.8088E-07 



MAG. 
270 DEG. 148353 
300 DEG. 5 
330 DEG. 419609 
360 DEG. 79057 



DO YOU WANT NORMALIZED DATA EXPRESSED IN DB?? Y 



MAG. 
30DEG-12.5766 
60DEG-6.0206 
90 DEG-4. 03484 
120 DEG 



MAG. 
150DEG-11.3119 

180 DEG-2. 0412 
210DEG-4.4167 
240DEG-131.03 



MAG. 
270 DEG-16. 5741 
300DEG-6.0206 
330 DEG-7. 5431 
360DEG-2.04119 
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GAMMA MATCH 

The next two programs are also antenna-related, being 
designed to provide matching information between an antenna and 
a transmission system. The first, the Gamma Match Program, will 
determine the necessary dimensions for a Gamma Match on the 
driven element of a Yagi antenna. However, the Yagi antenna is not 
the only applicable type of antenna. Vertical antennas can also be 
shunt fed using the same type of matching arrangement. 

As with the last program on antenna patterns, this program 
uses many smaller operations to obtain the final answer, computing 
the impedance of two parallel conductors, polar to rectangular 
conversions, and Smith Chart operations. Again, the complete 
program will be explained in small sections so that you can 
understand the operations and use any portions for your own 
programs. 

First, see Fig. 4-2. This diagram illustrates the gamma match 
and shows the designations used in the program. Each is relatively 
self-explanatory. Note that the spacing between the driven 
element and gamma rod is the spacing between the inner edges of 
the two pieces. Each of the spacings is entered in inches, the 
impedances in ohms, and the frequency in megahertz. 

Lines 90 through 110 of the Gamma Match Program compute 
the impedance step-up ratio between the gamma rod and driven 
element by the following formula: 



r=|l + - 



2S « 
lo g"df 

log — 



where S = center-to-center spacing 

d z = diameter of driven element 

dj = diameter of rod 
This ratio must be computed to compensate for the unequal 
diameters of the element and gamma rod which act as a section of 
transmission line. 

Since these two pieces do act as a section of transmission line 
that has been short circuited at one end, they will add some 
additional reactance that must be taken into account during the 
calculations. Line 120 computes the impedance of this section of 
transmission line. The formula for this equation is: 
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FREQUENCY-F 
IMPEDANCE RA.XA 



f 



GAMMA ' :SL - 



DE CAPACITOR 4^ 



Fig. 4-2. The gamma match and designations used in the program. 

45 2 -d 2 2 -d 1 2 
Z Q = 60 cosh" 1 2dd ohms 

The next step is to determine the increased impedance of the 
driven element over its center-point impedance caused by its being 
fed off center (tap point of gamma rod). Line 130 does double duty 
by not only calculating this increased impedance, but also multi- 
plies this product by the previously determined step-up ratio. At 
this point, RI and XI represent the real and imaginary components 
of the inital input impedance. Lines 140 and 150 calculate the 
impedance after the RI and XI values have been added to the 
reactance caused by the short-circuited section of the transmission 
line. These new values, RI(2) and XI(2), represent the real and 
imaginary components which will be rotated around the Smith 
Chart. The amount of rotation is determined by the length of the 
gamma rod vs. the frequency of operation. Lines 160, 170, and 180 
perform the actual rotation calculations. The formula for this 
rotation is: 



/rii-hxiiV 

"\RI2+jXI2J* 



Z=| ^T^n »R 



where RI1 = R 

XL1 = RX K +X t 

K L 

RI2 = R -X.3L 

L K 

XI2 = R i X K 

The division performed in line 180 is a complex type of division. 
That is, it is performed in a polar notation by dividing the mag- 
nitudes and subtracting the phase angles. 
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The next group of lines starts by converting the M3 and A3 
values back to the rectangular type of notation in line 190. Next, 
lines 210 and 220 convert the series equivalent of the impedance to 
a parallel equivalent. This is done so that the reactance caused by 
the short-circuited section of transmission line can be directly 
added to the parallel equivalent. Once this reactance is added, line 
230, the parallel equivalent, is converted back to the series form in 
lines 240 and 260. And finally, the real part of the impedance is 
multiplied by the system impedance to produce the final real and 
imaginary parts. 

In line 270, the final real part of the impedance is compared to 
the transmission system impedance. If the difference is less than 1 
percent, the value of the gamma capacitor is calculated and 
displayed. Otherwise, the program will change the value of T, the 
distance to the tap point, and go through the entire series of 
calculations again. Initially, the change to T is one inch. However, 
once the difference is less than 5 ohms, but more than 1 percent, T 
is either incremented or decremented by 0.25 inch until a 
calculation shows a difference of less than 1 percent. 

The final group of lines, lines 350 to the end, are used to 
change any of the initial parameters for recalculations. 
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Sample Gamma Match Program 

Determine the correct tap point and gamma capacitor value for 
an antenna and Gamma Match system with the following charac- 
teristics: 

Driven Element Diameter 1 inch 

Gamma Rod Diameter .375 inch 

Inside Spacing 4 inches 

Transmission System Impedance 50 ohms 

Real Part of Antenna Impedance 35 ohms 
Imaginary Part of Antenna Impedance -35 ohms 

Frequency 16 Megahertz 

Distance to Tap Point 24 inches 

(Each time the program goes through a set of calculations that does 
not produce a difference of less than 1 percent, the tap point 
distance goes back through its calculations until the right tap point 
is found.) 

THIS PROGRAM IS DESIGNED TO DETERMINE THE 
FINAL TAP POINT AND VALUE OF THE SERIES CAPACITOR 
IN A GAMMA MATCH. THE PROGRAM WILL TAKE THE 
INITIAL DATA AND GO THROUGH MULTIPLE ITERATIONS 
UNTIL A SUITABLE MATCH HAS BEEN DETERMINED. THE 
FINAL TAP POINT WILL BE WITHIN ONE QUARTER OF AN 
INCH. THIS WILL PROVIDE A SUFFICIENTLY ACCURATE 
STARTING POINT FORANACTUAL ANTENNA. 
ENTER DIAMETER OF ELEMENT? 1 
ENTER DIAMETER OF GAMMA ROD? .375 
ENTER SPACING BETWEEN ELEMENT AND ROD? 4 
ENTER TRANSMISSION SYSTEM IMPEDANCE? 50 
ENTER RESISTIVE PORTION OF ANTENNA IMPEDANCE? 
35 

ENTER REACTIVE PORTION OF ANTENNA IMPEDANCE? 
-35 

ENTER FREQUENCYIN MEGAHERTZ? 16 
ENTER DISTANCE IN INCHES FROM CENTER OF DRIVEN 
ELEMENT TO CENTER OF GAMMA ROD SHORTING 
STRAP? 24 

24 INCHES 21. 5639 OHMS 

25 INCHES 24.0441 OHMS 

26 INCHES 26.7077 OHMS 

27 INCHES 29. 5588 OHMS 
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28 INCHES 32.6002 OHMS 

29 INCHES 35.8341 OHMS 

30 INCHES 39. 2608 OHMS 

31 INCHES 42.8794 OHMS 

32 INCHES 46.6872 OHMS 
32.25INCHES 47.6684 OHMS 
32. 5 INCHES 48.6609 OHMS 

THE FINAL RESISTANCE IS 49.6648 (A DIFFERENCE OF 

LESS THAN 1 PERCENT). THE NEW TAP POINT IS AT 32.75 

INCHES. THE GAMMA CAPACITOR IS 93.4315 PICOFARADS 

AGAIN? (1=YES, 0=NO) 

?1 

ENTER THE PARAMETER YOU WANT TO CHANGE 

ELEMENT DIAMETER 1 

ROD DIAMETER 2 

SPACING 3 

SYSTEM IMPEDANCE 4 

ANTENNA RESISTANCE 5 

ANTENNA REACTANCE 6 

FREQUENCY 7 

?7 

NEW FREQUENCY? 

?21 

24 INCHES 44.7599 OHMS 

THE FINAL RESISTANCE IS 49. 9176 (A DIFFERENCE OF 

LESS THAN 1 PERCENT). THE NEW TAP POINT IS 25 

INCHES. THE GAMMA 

CAPACITORIS 71.0716 PICOFARADS 

AGAIN? (1=YES, 0=NO) 

?0 

READY 

As can be seen in this example, the initial set of information 
had the tap point on the driven element too close to the center of 
the element. The program went through the calculations until a 
suitable tap point was found. Each time through the calculations, 
the tap point and impedance were displayed to show the progres- 
sion. Finally, the program reached a point where the increments 
were changed to 0.25 inch until the final point was reached. Once 
the correct point was found, the frequency was changed to show 
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how the same physical conditions could provide a suitable match at 
a different frequency. As it turns out, at the new frequency of 21 
megahertz, only one tap point increment was required. Note the 
difference between the values for the gamma capacitor. 
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OMEGA MATCH 

The Omega Match is also a method of matching antennas to 
transmission lines. In addition to the series capacitor and gamma 
rod of the Gamma Match, another capacitor is added in parallel 
across the open end of the Omega rod (see Fig. 4-8). Though a 
seemingly uncalled for addition to the Gamma Match, the added 
capacitor allows a person to set the tap point at one spot 
permanently. Therefore, instead of trying to move the tap point to 
find the correct matching point, the tap is set at a single spot and 
the two capacitors are adjusted for a suitable match. Some may 
deem this an unnecessary complication with little to be gained. 
However, depending upon the size and installation of the antenna, 
it may be much easier to settle for the additional capacitor rather 
than to try to change the gamma rod tap point. This is especially 
true when the tap point is out greater than three feet from the 
center of the element. Plus, it is relatively easy to attach a small 
motor to the capacitors to provide complete, remote tuning 
capabilities. 

The program for the Omega Match is very similar to the 
previous Gamma Match program. In fact, it is essentially identical 
up through line 200. At this point, a new factor is added, C. The 
Omega capacitor, since it is in parallel with the short-circuited 
transmission line added directly to the reactance presented by this 
line. The value of C is either incremented or decremented until a 
suitable match is found. 

Electrically, the Omega Match functions essentially as the 
Gamma Match. For the Omega Match to be effective, though, the 
tap point must intentionally be set closer to the center of the driven 
element than would normally be required for a Gamma Match; i. e. , 
if the normal tap point for a Gamma Match was at 25 inches, the 
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Fig. 4-3. The additional capacities which distinguish the omega match 
from the gamma match. 
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Omega Match would require that the point be set less than 25 
inches. The closer the tap point is to the correct point for a Gamma 
Match, the smaller will be the value required for the Omega 
capacitor. In effect, by setting the tap point closer to the center, 
you are adding a large amount of inductive reactance due to the 
short-circuited transmission line effect. By placing the Omega 
capacitor across the open end of this transmission line you will be 
able to "tune out" the effects of the transmission line reactance. 
This essentially allows you to match the antenna via two 
capacitors. 

The formulas used in the Omega Match program are the same 
as those used for the Gamma Match and won't be repeated. 



108 



^ iN 



X 
g o 

co l— 

go <I 

© x 



U 

E 



jl 
CO 

LU 

X 

cl 

ll "■ 

X 

CO O 

IU !— 

D < 

_J 21 



X 

*-■ 






U 
SI 



_J o 



2: 

cc 

a:: 
a. 



CL 
<t 

CO 
UJ 

cr 
in 

CO 



H <I 



y 

a: 



a: < 
Q 

LU - H 

- 9 uj 

= O X 

H CL* LU 

xL —J 

UJ < UJ 

21 21 

LU 21 X 

-J <C UJ 

UJ IP UJ 

LL Ll H 

O C UJ 

iXi 

CL* CL 

111 UJ CD 

t_« i "^ 

LlJ LU s — i 

21 21 O 

<C <I < 

w w CL 

p Cl iXf 



CL <C 
•* X 

LU = 
U UJ 
X o 
<c x 
p <r 

LU P 

r^ a. LU 
- X cl 

= j-i X 

LU HH 

x x < 

<£ X X 

Q LU X 

UJ H LU 

9z 2: H 

i- <t X 

U_ 
X O Ll 
UJ o 

I- X 
CO O X 

> I—! O 

CO H •"• 

x o cl 

O CL o 

w CL 
CO LU 
CO > LU 

t— ! J— I 1> 

CO CO h- 
X 1-1 o 

<r co <c 

CC UJ LU 

h- ir a: 



LU 

> 

ir 
a 



a: 

LU 

H 
X 
UJ 

O 



M O 

h or 
tr ll 
lu 
x co 

< LU 
CD X 

lu a 

X X 



H 



> 

O UJ 

x cj 

LU X 

3 <r 

■3 t- 

LU CO 

cr m 

Ll P 



CL 



CO 

cd 



0:0:0:0:0:0:0:0: 

LU 



LU LU UJ 111 LU LU LU 

hh'rhl-Hhh 
!X X X! X! X X! X! 
LU LU LU UJ LU LU UJ 



*<Cj~ J^ Ei i - n ^ i^-* * Lwhh* h ii n w Lhhl, U— fin mil X^ttm 

CD hhhbDDDDDDD 

lu xxxclclclclclclclcl 

xca*-.w>-<zxzzxzzz 

OJCtUCCHHMHMHM 

r. O CL CL a. 

d3C9QQ£3£3C3C3 

i-H in n <t o *+ r\i ro <t* m o r- co 



h co 

K + 
O iN 

X \ 
en a: 

Q 

O + 
H IN 

UJ 

a 

h 

CD 



H 
X 

UJ 
X 
LU 

_j in 

Ui CD 



I > 

X * 

* > 

x ^ 

a: -2? 

»2» CO 

CO + 

+ > 

x ^ 

— CD 

CD O 

O _I 

-J II 

li ~* 

i-i >- 

x » 

~ LU 

cr q 

p * 

* CO 

CO * 

* «- 



.*-«. iN 
■N +■ 

ir q 



.-. CM 

«N <- 

^ lU 

UJ Q 

O -' 

-w + 

i ^ 

~ CO 

CO * 

* CO 

CO * 



11 > 

X 

IS 



Q 
! 

Ui 

n 
<t 

n 

«N 

ro 

m 
* 

Ll 

o co 



■r-1 

-J 

M 

■w 

rr 

i3» 
0) 

+ 

_J 

M 



O 



■>0 
II 

M 



a: 
P 

UJ 

p 

C-4 



x LU 



iN 



0-4 
-*- 

a: 
p 



n 
<*■ 
ro 

iN 

n 

<t . 

si in 

i * 

UJ Ll 

n * 

x 
< 

* 
-J 



iN 

•*- 
LU 
P «N 

— n 
i - 

Ll 



^ i CO 

j^«. ^-s. t^i 



CO 

-3- 



^ x 

* li 

M ^ 

X T-l 

I! ~ 



- If 

il H 

N JH M 

SQSQ 

t-h in n <t 



109 



- LU 

<$ co 

Q -J 

I OJ 

Ul © 

re n s oo <:ijjq 

<f CM < < H h- X X < 

n m 2: < o h h- o: 

r-j o £i 'N z U9 «t 

hccoo ~ h <r <t n " m n 

■ in in o o cd re <t = co • o 

in o th * ^ ^ oc m - uj en o 

* CM CM -J O •■ OUJOUJI- -J LU w 

ll - ■ m cm ^ ll x h s: z x cl 

* h- I s - * r*3 i h~ O O *-« Ll O = 
h- in IP ^ LL IP OH l£3 o o z ■- 
•w * * ^ * <f II O » 2: CL h-i — . 

z -^ -~- *> ^<i r^ o cl « = <r lu = ll 

< ^ ^v f-..( Q- o h CO > ■ CL UhX 

h th im re t>4 j isi z q: co cr < z ■"• * 

* w ^. tn re uj • lu h <r q o h in = ll 

M j-, (-; <j- m~| ^h $ x O tXi <£ Ll. CC {— * s 

a: cr qc n oo <r * od un h lu co cr a: lu < <i ■ 

u N -v ri m I s - <5 in <i sr lu <r h lu ll ^ co 

~ .*-. — Q j>^.^([>^ M ETULlDX Ll CO <f Q 

iN th iN . »N Q <t 04 Z A- O LU O Oh- «-« »-» *«-* < 

•^ -w w * * . ^4 . m * LL o 21 a O P ■ LT 

m hh i-i re n # » n co co cr huz (— re<: 

xxx < m <t n m * 2: lu 1— cl <: < <; z * ll 

«. w ^ ^ #<* *mill X O = LL _ — ^ N 'w- 

^zz z ■- ^ m ^ e: o I-* h» z: -^ x = w *-* o 

^ l_ |_ M CO Z -w IX. II = ^Oh LlCL\h 

<$■<:< co cr »-i + ce ~ ■»* lu hcoiion cr <i cl 

^ ii I! * \ CO ^ ^ h LL CO < ►-» X H 3 «* CL LU = 

1 h cm re CO * ~ H —- cr _J & = < t~i ■<* 

LU < <C E x 't <!■ X 03 h UJ COiTLLCOZ CO H = O 

re ■■ •■ ^ '-'is ^xz M o * »-s i-H w co ** 

Kt u -[ m 11 z 11 1 z » h co *- <i <c ui *-* = 

re ■ > cr h cl cu f- ~ or <i h h h h ^ lu x = co 

CM <a- <a- X <C X re <C O O. CM Z O <f Z < O H ^ *-» 

re ^ .-s. .« J an <j~ i )N sa J--I <T ^H (™j ^ (^ 

« n cm ^ h * re ii re *- o o cl ■ o co < ^ cc: 

m U u o -s" o cm m in ^ h cl <c re cl z h • ^ o 

* ^. ^ in < cm re <c <r *-* o u * o co ■~- & h 
ll ^ — re ■• ro . «• h* co up cl i>\ a. w MhctH 

* th cm m in tfi m in h x <<t v <H co z x o 

{-^■■^ <j- »<t**«Si^:© h-CDVh~<t LULUO< 

^, .^ i-t ,%| |^ ^ N U. u » J S LU M3 H a: O r 0- 

Z X X < h ^Hif:^*Nt-! LUEILULU3 tr tC <Z 

<. — ^ i S^<N8h »N m 11 N I O ^ I CL -J LU Q O 

K 4- + •»-! « H ♦- «^* <IiiNQh ^hE <!CLh" 

* CM CM <*OlXi*Zh-Xv n = 2: = = O Z = w < 

w ^ ^ H re bb dc <t < x co » ^ \'<a h < to h u ^c: ^ 

x ^ ^ rj < q: t <r h + o ^ n " z ^ z ll n <r uj 

>^ ,-. ^ <^ ^. x iN ^ \ <N U h l O h h-i LU ^ HZ 

1 rn in »» d*j H ^ co -h u ^ ^ ii h q: . cr cr x uj v <c o 

j ^ v ,n o co ixr o 1 cl m co cc o cl x o cl h x * a 

N HHZoxxo^Q:E£r^^ oh H-cr ui 

iiircr*v*-«''w.^-i-wj|^co««sH^Si> = 151 

^ ^ ^ ^ rj q: ^ ^ a. x ^ _j Qj ^t H < o v a: h» z: h- 

'N ^ ^ I Z Ct h E x h h n < "I- U I <I O h iN Z = 51 

w n n H |j ii 11 ii ii ii *w ii o- a. in m -^ < 

w th cm re cc co <r cl h- m co u. ll ii ll <r <r ll - n cr u? . 

OCZEIErirEEKxEiiaMUHlOHl CL^c CO 

SQSQQSQSiSQSQ^ CO C> Q 

m^MDi>QH cm re <f Ui <3 nQ -0 n0 n 

^ ^ ^h ^1 ^, CM CM CM CM CA CM CM CM CM CM '.M 



110 





s 




cr <: 




LU 




LU X 




LD 




Q M H 




Z 




3 h- 




< 




H l-CLD 




X 




3 0£ DZH 




o 




H Q QL m Z 
Z Z « t»t 




C' 




H f~ (—!=■■ 


■3* 


1" 




« 3 •■ o- = Ll f- 




= = 


= = = = = = 


S CL -r-i = LU O' 


K 


h *h c^ rn <t in <i n oo 


c- Z UJ c- O LU H 1- 


3 


z 




£T w LU Z O 3 3 


CL 


< 




uj « ho<za.o. 


Z 


3 


LU 


H * 3 Z H <£ Z Z 


i„j 




a LU 


LUe-CL<£tf!f-HHf-i 


■ B 


3 IT 


lu z a 


I H Z P h o ■■ « 


= 


o yj 


a <r z 


<LU»-sLU(Xl<= = 


•—». 


> f- 


z h <r 


t-i h- «« a. lu ui »• o* 


O 1 


LxJ 


<: co h 


quj= z a: a: > h 


z 


££ z tr 


O m o 


z in h-i a z 


ii 


Hi < LU 


lu en <: 


!- < ID <I<Zh 


s 


I — H j — 


CL LU LU 


Z i-h Z Z Z Z Hi O 




LU Q ill 


z ctr ir > h» 


IxiQHliJZZDa- 


». 


z z 


HH O Z 


Z O H LU LU & 


en 


< !-■ < is <<tzh-i 


LUQ<C0HKLUQ. 


LU 


(TZh 


Z Z Z Z LU O 


J O i > Z Z £T < 


> 


<t 111 P 


HUjZZDCL 


Luct:cnco<:<:LLH- 


II 


a. z 


O H LU LU O 




^-i 


u a 


<. co h- h lu a. 


|iM*fr» #**■•** HPiTt* II iiimV **»•■> **j**J* -;^ mumV 


•w 


Q U J O 


a. > z z tr < 


LUiUlULULULULULU 




zxujcrixiiXi<:<:LLH 






lu h- = = 


= = = = = = 


= = = = = = = = 


t>- 






f— .(— u™{_j_j— . |«-[— 


ZZQIQQQSQQQQ 


Xv> Ji— £■ Ml J™ <C" >bw J5™ 


HH [JLJ LU -r-! T-i 

<J J U v v — -^ ^ -^ ^ -w 


t-1 t-I *r-( TT-i rH -»H 


hhhhmhhh 

irartrtrtrarcrot: 

CLCLCLCLIX.CLQLCL 


U?h-ZQJCuCQCQCGCQCQC£j 


< 


uj<:<:<c<<:<:<3:<:z 


in 


= 


Q - | | — | | J — j — J | 


^-i c-4 re <r- in o h- co co 


H 


H HHHI— f~ h-HHHH 


1! II 11 II i! II fl II 


Z 


o*zzzzzzzzzr> 


zzzzzzzzo 


h-i 


i — i t — i i— i 


H H H M H M U- 


L_ 


tru.o:Q:DCQ:a:Q:a:£KQ:zLj.LLLLLLLLLLLLLi-b 


CL 


HQ.O.Q.0LQ.CLQ.IICLH 


HH »— iHHh tl— 1 H-i (— 11X3 



^^iNn<fin%oNoochQ^iNn<}-in<iNcoi>isi 
intnininininminmm>o^NONO-<3<!NO<isO<3h- 



111 



Omega Match Example 

Determine the values for the Omega and Gamma capacitors 
for an antenna with the following characteristics: 

Driven Element Diameter 1 inch 

Omega Rod Diameter . 25 inch 

Spacing 4 inches 

System Impedance 75 ohms 

Resistive Part of Antenna Impedance 30 ohms 

Reactive Part of Antenna Impedance - 30 ohms 

Frequency 25 Megahertz 

Tap Point Distance 24 inches 

THIS PROGRAM WILL DETERMINE THE VALUES FOR THE 
SHUNT AND SERIES CAPACITORS IN AN OMEGA MATCH. 

ENTER DIAMETER OF ELEMENT? 1 

ENTER DIAMETER OF GAMMA ROD? .25 

ENTER SPACING BETWEEN ELEMENT AND ROD? 4 

ENTER TRANSMISSION SYSTEM IMPEDANCE? 75 

ENTER RESISTIVE PORTION OF ANTENNA IMPEDANCE? 

30 

ENTER REACTIVE PORTION OF ANTENNA IMPEDANCE? 

-30 

ENTER FREQUENCY IN MEGAHERTZ? 25 

ENTER DISTANCE IN INCHES FROM CENTER OF DRIVEN 

ELEMENT TO SHORTING STRAP? 24 

83.3972 OHMS 

80. 3289 OHMS 

77.4027 OHMS 

74.6119 OHMS 

THE TAP POINT IS TOO FAR OUT FOR AN OMEGA MATCH. 
TRY THE COMPUTATIONS AGAIN WITH ANOTHER TAP 
POINT. 
ENTER THE PARAMETER YOU WANT TO CHANGE 



ELEMENT DIAMETER 


1 


ROD DIAMETER 


2 


SPACING 


3 


SYSTEM IMPEDANCE 


4 


ANTENNA RESISTANCE 


5 


ANTENNA REACTANCE 


6 


FREQUENCY 


7 


TAP POINT 


8 
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?8 

NEW TAP POINT? 

?22 

67.5678 OHMS 
70.0705 OHMS 
72.6968 OHMS 
75.4535 OHMS 

THE FINAL RESISTANCE IS 75.4535 (A DIFFERENCE OF 

LESS THAN 1 PERCENT). THE TAP POINT IS AT 22 INCHES. 

THE GAMMA CAPACITOR 

IS 53.4325 PICOFARADS. THE OMEGA CAPACITOR IS 4 

PICOFARADS. 

AGAIN? (1=YES,0=NO) 

?1 

ENTER THE PARAMETER YOU WANT TO CHANGE 
ELEMENT DIAMETER 1 
ROD DIAMETER 2 

SPACING 3 

SYSTEM IMPEDANCE 4 
ANTENNA RESISTANCE 5 
ANTENNA REACTANCE 6 
FREQUENCY 7 

TAP POINT 8 

?8 

NEW TAP POINT? 
?23 

75.3624 OHMS 
THE FINAL RESISTANCE IS 75.3624 (A DIFFERENCE OF 
LESS THAN 1 PERCENT). THE TAP POINT IS AT 23 INCHES. 
THE GAMMA CAPACITOR 

IS 54.0296 PICOFARADS. THE OMEGA CAPACITOR IS 1 
PICOFARADS. 
AGAIN? (1=YES, 0=NO) 
?0 
READY 

This example was selected to illustrate several features of 
this Omega Match program. In the first run through the data, the 
program found that the tap point was too far out from the center of 
the driven element. This means that the reactive element seen at 
the open end of the short-circuited transmission line was capacitive 
rather than inductive. This would have required a variable 
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inductance to act as the compensating element. Next, when the tap 
point was changed to 22 inches, a suitable match could be found. 
You'll note that the value for the Omega or parallel capacitor was 
extremely small. This indicated that the tap point was close to that 
required for a normal Gamma Match. The last run through the data 
was with the tap point at 24 inches. This produced a value for the 
Omega capacitor of 1 picofarad, illustrating the fact that as the tap 
point gets closer to that required for a Gamma Match, the value of 
the Omega capacitor gets increasingly smaller. Conversely, as the 
tap is moved in, the Omega capacitor gets larger. 
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T-PAD 

The next three programs deal with three types of attenuators 
or pads: a T-pad, a Pi-pad, and a minimum loss pad. The T and Pi 
are named for their characteristic shapes. Each pad will provide a 
known amount of attenuation after determining the values by using 
the appropriate program. 

Before going into the specifics of the T-pad, I do want to 
illustrate one useful property of attenuators that is oftentimes 
overlooked. In addition to providing known amounts of attenuation 
in a match system, pads can be used to provide impedance 
matching. Consider Fig. 4-4: this diagram shows a VSWR 
measuring device connected to some component that exhibits 
some amount of return loss (VSWR). Return loss is a measure of 
VSWR; the higher the return loss, the lower the VSWR. A return 
loss of zero dB is an infinite VSWR. A 1.4: 1 VSWR corresponds to 
areturn loss of approximately 15 dB. 

In measuring return loss, a small amount of power is applied to 
the device under test. Some of this power is reflected back to the 
VSWR measuring device as a function of the VSWR or impedance 
match. The return loss is the difference, in dB, between the 
applied and reflected power. 

Now, consider Fig. 4-5. In this case the applied power is 
mostly absorbed by the termination. The difference between 
applied and reflected power is very large, producing a high return 
loss. Next look at Fig. 4-6. Now, the 50-ohm load has been 
replaced with an open circuit. This presents an infinite mismatch, 
causing all the applied power to be reflected back, producing a 
return loss of zero dB. Finally, in Fig. 4-7, a 10-dB attenuator has 
been placed between the VSWR measuring device and the open. 
The applied power, when passing between the VSWR measuring 
device and the open, is attenuated 10 dB. At the open, the applied 
power, now down 10 dB, is reflected back. In the return path back 
to the VSWR measuring device, the reflected power is also 
attenuated by 10 dB. Therefore, the difference between the 
applied and reflected power is 20 dB, twice the value of the pad. 
Therefore, by simply inserting a 10-dB pad before the open 
improves the return loss from to 20 dB. This corresponds to a 
decrease in VSWR from infinitely high to about 1.2: 1. 

This technique can be applied throughout circuitry to help 
maintain system and circuit impedance levels. Normally, in today's 
50-ohm system, 50-ohm pads are used. However, nothing pre- 
vents one from designing pads for different impedance levels. 
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Fig, 4-4. A VSWR measuring device connected to a component which 
measures return loss. 

Another good example is the output of a doubly balanced mixer. 
For best intermodulation performance, the output of the mixer 
should see 50-ohms for all frequencies. A pad could be inserted 
between the mixer and the next device to ensure a good 50-ohm 
termination at all frequencies. Nothing is free, however. This will 
only work if the system has enough gain to make up for the loss in 
the attenuator. 

Basically, this is a very simple program. The values for the 
resistors are computed from the following common-day formulas: 



K1 " z yo 05A -MJ 



R2 = 



R11Z 

10 05A - 



where Z = impedance 

A = attenuation value 
Once the theoretical values have been determined, the only 
problem is to obtain these values with practical components. The 
second part of the program lets you input the real values that will be 
used and then computes the actual attenuation. 



VSWR 

MEASURING 

DEVICE 



APPLIED POWER 

"El^^ 

<^-4 50-OHM 

/ LOAD 
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Fig, 4-5. How a VSWR measuring device with a 50 ohm load can absorb 
applied power. 
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Fig. 4-6. A VSWR measuring device with the 50 ohm load replaced by an 
open circuit 

T-Pad Program 



5 'T PAD CALCULATIONS 

10 CLS 

20 PRINT" ENTER DESIRED ATTENUATION 

VALUE" 
30 INPUT A 

40 PRINT-ENTER SYSTEM IMPEDANCE" 
50 INPUT I 
60 R1 = I*< <10H.05*A)-1)/<10* <.05*A 

) + l)) 

70 R2=<R1 + I)/<10M.05*A>-1) 

80 PRINT n Rl = ,, Ri ,, OHMS% " R2=" R2"0HMS 
ii 

90 PRINT"ENTER THE ACTUAL VALUE TO 

BE USED FOR Rl" 
100 INPUT Rl 



VSWR 

MEASURING 
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50-OHM 
SYSTEM 
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POWER 
~p 
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PAD 
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Fig. 4-7. A VSWR measuring device with a 10 dB attenuator placed be- 
tween the VSWR and the open circuit. 
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110 PRINTENTER THE ACTUAL VALUE T 
BE USED FOR R2" 

120 INPUT R2 

130 Z=(R2*(I+R1) )/(Rl+R2+I)+Rl 

140 PR I NT "WITH ONE BIDE ACTUALLY T 
ERMINATED IN"I"OHMS» THE OPPOS 
ITE"CHR*<13)"SIDE OF THE PAD W 
ILL APPEAR AS "Z" OHMS." 

150 E=R2/(I+R1+R2) 

160 DB=20*(LOG(1/E)/LOG<10> ) 

165 PRINT 

170 PRINT" WITH Rl ="Rl"OHMS AND R2 
="R2"0HMS» THE ACTUAL ATTENUA 
TI0N"CHR$(13)"IS"DB"DB. " 
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50-OHM 
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Fig. 4-8. The component designations for a 50 ohm pad with 17 dB of 
attenuation. 



T-Pad Examples 

You need a 50-ohm pad with 17 dB of attenuation. What are the 
required values? (See Fig. 4-8 for the component designations.) 
ENTER DESIRED ATTENUATION VALUE 
?17 

ENTER SYSTEM IMPEDANCE 
?50 

Rl = 37.6229 OHMS R2 = 14.413 OHMS 

ENTER THE ACTUAL VALUE TO BE USED FOR Rl 
?39 

ENTER THE ACTUAL VALUE TO BE USED FOR R2 
?10 

WITH ONE SIDE ACTUALLY TERMINATED IN 50 OHMS, 
THE OPPOSITE SIDE OF THE PAD WILL APPEAR AS 47.9899 
OHMS. 

WITH Rl = 39 OHMS AND R2 = 10 OHMS, THE ACTUAL 
ATTENUATION IS 19. 9127 DB. 

This example readily illustrates the usefulness of this pro- 
gram. Since the required resistors are not of normal value, 
practical components have been substituted to build the T-pad. In 
this example, using real-world components have changed the 
attenuation almost 3 dB. It would be better to use a 15-ohm resistor 
in place of the 10-ohm resistor. This would decrease the actual 
attenuation to 16.8 dB vs 19.9. The impedance would also change, 
actually getting closer to 50 ohms. 



119 



PI-PAD 

The next program covers the Pi-pad. This program is very 
similar to the previous T-pad program, providing the computed 
values and then allowing for real-world components. The best 
approach would be to determine the desired attenuation value and 
then to run both programs. From this dual information you can then 
select the best pad according to the components on hand. 

The equations used for determination of the resistor values in 
the Pi-pad are: 



1 10 05A + A 

R1= \10 05A - 1/ R2 = 



ZR1(10 05A -1) 



Z + Rl 
Sample Pi-pad Program 

5 'PI PAD 

10 CLS 

20 PR I NT "ENTER DESIRED ATTENUATION 

VALUE 
30 INPUT A 

40 PRINT"ENTER SYSTEM IMPEDANCE" 
50 INPUT Z 
60 R1=Z*(10* <.05*A)+1)/<10M.05*A> 

-1) 
70 R2=Z*R1*<18M.05*A>-1>/<Z+Ri> 
80 PRINT" THE REQUIRED RESISTOR FOR 
Rl IS"Rl"OHMS» WHILE R2 SHOULD 
BE "R2" OHMS 
90 PR I NT "ENTER THE ACTUAL VALUE TO 

BE USED FOR Rl" 
100 INPUT RN1 
110 PRINT"ENTER THE ACTUAL VALUE T 

BE USED FOR R2" 
120 INPUT RN2 
130 RI=(Z#<RN1+RN2)+<RN1*RN2> )/(Z+ 

RN1) 
140 RT=R1*RI/(RI+R1) 
150 PRINT" WITH ONE SIDE ACTUALLY T 

ERMINATED IN" Z" OHMS ? THE OPPOS 

SIDE WILL APPEAR AS "RT" OHMS 
160 E=Z*RN1/(RI*Z+RI*RN1> 
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170 EF=-20*<LOG<E)/LOG(10> ) 
180 PRINT"THE ACTUAL ATTENUATION W 
ILL BE"EF"DB. 
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Pi -PAD 
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Fig.4-9. Resistor designations required to determine the value of resistors 
needed for a 12 dB pad on a 600 ohm system. 

Sample Pi-pad Program 

What value resistors are required to make a 12-dB pad for a 
600 ohm system? See Fig. 4-9 for the resistor designations. 
ENTER DESIRED ATTENUATION VALUE 
?12 

ENTER SYSTEM IMPEDANCE 
?600 

Rl = 1118.97 R2 = 1002.54 

ENTER THE ACTUAL VALUE TO BE USED FOR Rl 
? 1200 

ENTER THE ACTUAL VALUE TO BE USED FOR R2 
?1000 

WITH ONE SIDE ACTUALLY TERMINATED IN 600 OHMS, 
THE OPPOSITE SIDE OF THE PAD WILL APPEAR AS 612.6 
OHMS. 

WITH Rl = 1200 OHMS AND R2 = 1000 OHMS, THE ACTUAL 
ATTENUATION IS 12 DB. 

You can see from this example that as the system impedance 
goes higher, the pad becomes more tolerant to resistance values. 
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Fig. 4-10. Values for 2 resistors and their loss through the pad. 

MINIMUM-LOSS PAD 

The last attenuator program is a special case known as a 
minimum- loss pad. The pad is characterized by the ability to match 
impedances, whereas the T- and Pi-pad are able to swamp out 
impedance mismatches. However, the minimum- loss pad is also 
characterized by a high impedance loss through the pad. If the high 
loss can be tolerated, this is one excellent way of matching 
between two circuits. The program will provide the values for the 
two resistors (see Fig. 4-10 for the component designations) and 
the loss through the pad. The only other stipulation is that the pad 
is not bi-directional; it can only be installed in the circuit in one 
way. 

The equations to determine the component values for the 
minimum loss pad are: 

Rl = Z x (1 - Z/L^ R2 = Z 2 / (1 - (Z/Z,))* 
Minimum Loss Pad Program 

5 'MINIMUM LOSS PAD 

10 CLS 

20 PRINT "ENTER THE IMPEDANCES TO 

BE MATCHED (HIGH? LOW) 
30 INPUT Z<I)»Z(2) 
40 R<l)=ZQ)*SQR<i-Z<2)/Z<l>) 
50 R(2)=Z<2)/SQR(1»Z(2)/Z(1>) 
60 PRINT"Rl = l, R(i) ,, OHMS u 
70 PRIN^'Rl^'R^T'OHMS" 
80 L=LOG ( < Z ( 1 > ~R ( 1 ) ) / Z ( 1 ) ) /LOG (10) 
90 PRINT"THE LOSS IS ^B^L-DB" 
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Fig, 4-11. Matching a 10 ohm output to a 50 ohm input. 

Example Minimum Loss Program 

As seen in Fig. 4-11, you wish to match the 10-ohm output of 
the first stage to the 50-ohm input of the second stage which is an rf 
amplifier. The amplifier has 15-dB gain. What are the values of the 
matching resistors and overall gain from the input of the first stage 
to the output of the second stage? 

ENTER THE IMPEDANCES TO BE MATCHED (HIGH, LOW) 
? 50,10 

Rl = 44. 7214 OHMS 
R2 = l 1.1803 OHMS 
THE LOSS IS -19. 529 DB. 
Overall gain, therefore, is-4.5dB. 
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PARALLEL/SERIES CONVERSIONS 

As mentioned in the section of the Gamma Match program, 
parallel/series conversions are used during impedance transfor- 
mations. There are two methods of presenting impedance data: 
series format and parallel format. In both cases, they consist of a 
real, resistive component in either series or parallel, with an 
imaginary, reactive component. Though the values are different, 
the overall impedance is still the same since this is just two 
different ways of presenting the same data. Generally, most 
information is presented in series format. However, in some 
cases, as in particular when paralleling components, it is easier to 
work with the data in a parallel format. In this way values can be 
added directly. 

One area where this parallel presentation of impedance 
information is prevelant is the output of rf power transistors. In 
many cases manufacturers will present the transistor's input 
characteristics in series form and the output characteristics in a 
parallel form. 

The equations for series/parallel conversion are: 
R 

n 




R=R(1 + (X/R) 2 ) 

R 

X = —£ — 

p X/R 

s s 

Note that when using the plus and minus signs for the reactances 
they will carry through the equations to give the correct sign for the 
answer. 

Parallel/Series Program 

5 MTOS AND S TO P CONVERSION 
10 CLS 

20 PRINT"ENTER •'"S"" IF CONVERSION 
IS PARALLEL TO SERIES OR ""p"" 
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IF" CHR$< 13) "CONVERSION IS SERI 
ES TO PARALLEL" 
30 INPUT C$ 
40 IF C*="P"GOTO 130 
50 PRINT" ENTER THE PARALLEL RESIST 

ANCE> RP (OHMS)" 
60 INPUT RP 
70 PRINT "ENTER THE PARALLEL REACTA 

NCE, XP (OHMS)" 
80 INPUT XP 
90 RS=RP/(1+(RP/XP)*2> 
100 XS=RS*RP/XP 
110 PRINT"THE SERIES RESISTIVE PAR 

T IS" RS" OHMS AND THE SERIES" CH 

R$( 13) "REACTIVE IS" XS" OHMS 
120 END 
130 PRINT"ENTER THE SERIES RESISTA 

NCE ? RS (OHMS)" 
140 INPUT RS 
150 PRINT"ENTER THE SERIES REACTAN 

CE» XS (OHMS)" 
160 INPUT XS 
170 RP=RS*(1 + (XS/RSH2) 
180 XP=RP/(XS/RS) 
190 PRINT"THE PARALLEL RESISTIVE P 

ART IS" RP" OHMS AND THE PARALLE 

L"CHR$( 13) "REACTIVE PART IS"XP 

" OHMS 
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Parallel/Series Examples 

The optimum source impedance for best noise figure of a 
microwave transistor is 12.5 +J0.5 ohms. What is this same 
impedance expressed as a parallel equivalent? 
ENTER S IF CONVERSION IS PARALLEL TO SERIES OR P 
IF CONVERSION IS SERIES TO PARALLEL 
?P 

ENTER THE SERIES RESISTANCE, RS (OHMS) 
?12.5 

ENTER THE SERIES REACTANCE, XS (OHMS) 
?.5 

THE PARALLEL RESISTIVE PART IS 12.52 OHMS AND THE 
PARALLEL 

REACTIVE PART IS 313 OHMS 
Express 80-J90 ohms as a series equivalent. 
ENTER S IF CONVERSION IS PARALLEL TO SERIES OR P 
IF CONVERSION IS SERIES TO PARALLEL 
?S 

ENTER THE PARALLEL RESISTANCE, RP (OHMS) 
?80 

ENTER THE PARALLEL REACTANCE, XP (OHMS) 
?-90 

THE SERIES RESISTIVE PART IS 44.6897 OHMS AND THE 
SERIES 
REACTIVE PART IS -39.7241 OHMS 
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MATCHING NETWORKS 

The four programs to follow will solve for the required 
reactive values in the following impedance matching networks. 
Since the answers are given in reactive values, they can be scaled 
to any frequency by using the appropriate reactance formula and 
frequency. The first program is for network A which is shown in 
Fig. 4-12. 

This network is only used when the device that is to be 
matched has a series real part less than the load impedance. If you 
were to run this program several times you would see that as the 
real part approaches the load impedance, the reactance of CI 
increases toward infinity. This network, though, can be used for 
RF power transistor design since their input and output impe- 
dances are generally very small. 

The equations to solve network A are: 

XL = QR + X c 

XC2=AR L 

xci = (B/AXB/Q) = B 
(B/AMB/Q) Q-A 



whereA= y^) _, 



B = R(1 + Q 2 ) 

Network A Program 

5 'MATCHING NETWORK A 

10 CLS 

20 PRINT" SELECT A VALUE FOR 0" 

30 INPUT Q 



NETWORK A 



XC 



DEViCETOBE 
MATCHED 




Fig. 4-12. Network A, for which Program 4-8 is designed. 
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40 PR I NT "ENTER THE LOAD IMPEDANCE" 

50 INPUT RL 

60 PRINT "WITH THE IMPEDANCE TO BE 
MATCHED IN A SERIES FORMAT, ENT 
ER"CHR*(13)"THE RESISTIVE (OHMS 
) AND REACTIVE ELEMENTS (OHMS). 

II 

70 INPUT R»XC 

80 M=SOR( <R*(1+Q*2)/RL)-1) 

90 N=R*(1+Q*2) 

100 XL(1)=Q#R+XC 

110 XC(1)=N/(Q-M) 

120 XC(2)=M*RL 

130 PRINT"XL="XL(l)"OHMS" 

140 PRINT"XCl="XC(l)"OHMS" 

150 PRINT"XC2="XC(2)"0HMS" 



129 



Network A Example 

Using the diagram shown in Fig. 4-13, determine the 
reactance values for the matching components. 
SELECT A VALUE FOR Q 
?3 

ENTER THE LOAD IMPEDANCE 
?75 

WITH THE IMPEDANCE TO BE MATCHED IN A SERIES 
FORMAT, ENTER THE RESISTIVE (OHMS) AND REACTIVE 
ELEMENTS (OHMS). 
P20.-30 
XL = 30 OHMS 
XC1 = 117. 027 OHMS 
XC2 = 96. 8246 OHMS 
READY 

You'll notice that, in this example, the reactance of the inductor 
matches the reactance of the output capacitance. In this way, the 
reactive element is cancelled and the two capacitors transform the 
real part of the output impedance up to the 75-ohm load impedance. 

To illustrate the problem encountered with a series real part 
close the load impedance, run the following information: Q = 3, ZL 
= 75, Real Part = 74, and Imaginary Part = -30 
SELECT A VALUE FOR Q 
?3 

ENTER THE LOAD IMPEDANCE 
?75 

WITH THE IMPEDANCE TO BE MATCHED IN A SERIES 
FORMAT, ENTER THE RESISTIVE (OHMS) AND REACTIVE 
ELEMENTS (OHMS). 
?74,-30 



-30 



NETWORK Q =3 

XL 96,82 fl 



30a 



117.03, 



DEVICE TO BE 
MATCHED 



o- 



,XC1 



\XC2 " L RL 

7sa 



ffi 



Fig. 4-13. Network Q-3, to be used in Network A problem. 
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XL = 192 OHMS 

XC1 = 33177. 5 OHMS 

XC2 = 223.327 OHMS 

When applied to the standard capacitive reactance formula, the 

33000+ ohms create a very small capacitor value. 
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NetworkB Example 

Network B, as shown in Fig. 4-14 is the Pi network which has 
been widely used to match the output of a vacuum-tube transmitter 
to a transmission line. 

The main limitation of this network deals with the selection of 
Q and the matching range. When a low Q is selected, the value of 
the parallel real part of the output impedance must also be low. 
This will be illustrated in the following examples. Notice that the 
program calls for the data to be input in a series equivalent format. 
If you wish to use the series format, delete line 71. Otherwise, 
leave line 71 in and input the data in a parallel format. The type of 
format you use will depend generally upon the device being 
matched. For example, some power transmitters present their 
output impedance data in a parallel equivalent, while others use a 
series equivalent format. Power tubes, on the other hand, use a 
parallel format. 

The equations to solve for the required component values are: 

XC1 = R/Q 



XC2 — RL 



R/RL 



'(Q 2 +1)~(R/RL) 



yt _QR + (RRL /XC2) 

xi WTi 



Network B Program 

5 ' MATCHING NETWORK B 

10 CLS 

20 PRINT-SELECT A VALUE FOR Q" 



NETWORK B 



XC 



ft? 




rr? 



rn 



DEVICETOBE 
MATCHED 



Fig. 4-14. Network B, for which Program 4-9 is designed. 
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30 INPUT Q 

40 PRINT" ENTER THE LOAD IMPEDANCE" 

50 INPUT RL 

60 PRINT" WITH THE IMPEDANCE TO BE 
MATCHED IN A SERIES FORMAT » ENT 
ER"CHR*(13)"THE RESISTIVE (OHMS 
) AND REACTIVE ELEMENTS (OHMS). 

II 

70 INPUT R»XC 

71 RP=R:XP=XC:GOTO 100 
80 RP=R*(1+(XC/R)*2) 
90 XP=RP/(XC/R> 

100 XC(l)=RP/0. 

110 XC(2)=RL*S0.R( (RP/RL)/( (QA2+1)- 

(RP/RL) ) ) 
120 XL=(G>*RP+(RP*RL/XC(2) ) )/(GH2+l 

) 
130 PRINT"XCl="XC(l)"OHMS" 
140 PRINT" XC2="XC(2) "OHMS" 
150 PRINT" XL=" XL" OHMS" 
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Sample Network B Problems 

Match the output of the vacuum tube transmitter shown in Fig. 
4-15 using a network with a Q of 5 and a load impedance of 50 ohms. 
(Since the diagram shows the output in a parallel format, leave line 
71 in the program.) 
SELECT A VALUE FOR Q 
?5 

ENTER THE LOAD IMPEDANCE 
?50 

WITH THE IMPEDANCE TO BE MATCHED IN A SERIES 
FORMAT, ENTER THE RESISTIVE (OHMS) AND REACTIVE 
ELEMENTS (OHMS). 
? 2000, 10 

?FC ERROR IN 110 

This example shows that with the high real part of the impedance to 
be matched and with the low Q of the network, there is no solution 
to the problem. To achieve a solution, the circuit Q will have to be 
raised. 

SELECT A VALUE FOR Q 
?15 

ENTER THE LOAD IMPEDANCE 
?50 

WITH THE IMPEDANCE TO BE MATCHED IN A SERIES 
FORMAT, ENTER THE RESISTIVE (OHMS) AND REACTIVE 
ELEMENTS (OHMS). 
? 2000, 10 

XC1 = 133.333 OHMS 
XC2 = 23. 1869 OHMS 
XL = 151.826 OHMS 



NETWORKQ = 5 



2000 



10 



rrr 



XL 
_rYYY*> ■ o 



m 

XC1 



DEVICETOBE 
MATCHED 



rn 

XC2 
NO SOLUTION 




Fig. 4-15. Network Q-5, to be used in Network B problem. 
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NETWORK Q = 2 
XL 



2871 

n 



DEVICE TO BE 
MATCHED 



ft? 

XC1 
12.5ft 



-75ft 



ft? 

XC2 
20.04ft 



Fig. 4-16. Network Q-2, to be used in Network B problem. 

In actual practice, the value for CI = CI - Output C Depending 
upon the frequency and the matching components, CI can actually 
be eliminated and the output capacitance of the tube will suffice for 
the matching capacitor. This, generally, happens only at frequen- 
cies greater than 40 to 50 megahertz. 

Match the output impedance of the circuit shown in Fig. 4-16 
to 50 ohms using a network with a Q of 2. 
SELECT A VALUE FOR Q 
?2 

ENTER THE LOAD IMPEDANCE 
?75 

WITH THE IMPEDANCE TO BE MATCHING IN A SERIES 
FORMAT, ENTER THE RESISTIVE (OHMS) AND REACTIVE 
ELEMENTS (OHMS). 
?25,5 

XC1 = 12.50HMS 
XC2 = 20. 0446 OHMS 
XL = 28. 7083 OHMS 

Network C 

Network C can be implemented by either of the two methods 
shown in Fig. 4-17. 

This program, though, is designed to solve for the matching 
network when the data is presented in a series equivalent format. 
In either case, the real part of the impedance to be matched must be 
less than the load impedance. Practically, this configuration yields 
the best component values when the real part is of a low value. 

The equations to solve network C are: 
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NETWORK C 
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XL2 
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/77 




DEVICE TO BE 
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Fig. 4-17. Two methods by which Network C can be implemented. 

XL1=XC 
XC1 = QR 



XC2 



= RL S / 



Rl 



XL2 = XC1 + 
Parallel Format 



RL-R 



/RRL\ 
\XC2| 



XC1 = QR1 



XC2 



=rlJ 



R 



RL-R 



XL2=XC1+(||^\ +XC 
Series Format 

network C Program 

5 'MATCHING NETWORK C 

10 CLS 

20 PRINT"SELECT A VALUE FOR Q" 
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30 INPUT 

40 PRINT" ENTER THE LOAD IMPEDANCE" 

50 INPUT RL 

60 PRINT" WITH THE IMPEDANCE TO BE 
MATCHED IN A SERIES FORMAT* ENT 
ER"CHR*(13)"THE RESISTIVE (OHMS 
) AND REACTIVE ELEMENTS (OHMS). 

I) 

70 INPUT R>XC 

100 XC(1)=R#0 

110 XC(2)=RL*SQR(R/(RL-R> ) 

120 XL(1)=XC(1)+<R*RL/XC<2) )+XC 

130 PRINT"XL="XL(1)"0HMS" 

140 PRINT" XC1="XC(1) "OHMS" 

150 PRINT" XC2="XC(2) "OHMS" 



137 



Sample Network C Programs 

Match the circuit shown in Fig. 4-18 using Network C with a Q 
of 5. 

SELECT A VALUE FOR Q 
?3 

ENTER THE LOAD IMPEDANCE 
?50 

WITH THE IMPEDANCE TO BE MATCHED IN A SERIES 
FORMAT, ENTER THE RESISTIVE (OHMS) AND REACTIVE 
ELEMENTS (OHMS). 
? 60,300 

?FC ERROR IN 110 
READY 

This illustrates the point that the load impedance must be higher 
than the series real part. Now, try solving the network by raising 
the load impedance to 75 ohms 
SELECT A VALUE FOR Q 
?3 

ENTER THE LOAD IMPEDANCE 
?75 

WITH THE IMPEDANCE TO BE MATCHED IN A SERIES 
FORMAT, ENTER THE RESISTIVE (OHMS) AND REACTIVE 
ELEMENTS (OHMS). 
? 60,300 

XL = 510 OHMS 
XC1 = 180 OHMS 
XC2 = 150 OHMS 
READY 



10 



If 



DEVICE TOBE 
MATCHED 



NETWORK Q = 3 



XC1 

-\(r 

180H 



XL2 
51011 



-r 



XC2 _ 
1500'X 




Fig. 4-18. Network Q-3, to be used in Network C problem. 
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xc 



NETWORK D 



XL1 



XL2 



XC1 



m 



> 



DEVICE TO BE 
MATCHED 



Fig. 4-19. Network D. 

Network D 

The last program covers Network D, the classical T type of 
network. An illustration of Network D is given in Fig. 4-19. The 
advantage of this network is that it will match load impedances that 
are both less than and greater than the series real part of the 
impedance to be matched. 

Equations to use this network are: 

XL1 = RQ + XC 
XL2 = RLB 



XC1: 



Q + B 



whereA = R(l + Q 2 ) 



B 



=vfp 



Network D Program 

5 'MATCHING NETWORK D 

10 CLS 

20 PRINT "SELECT 

30 INPUT 

40 PRINT 

50 INPUT 

60 PRINT 



A VALUE FOR 0." 



ENTER THE LOAD IMPEDANCE" 

RL 

WITH THE IMPEDANCE TO BE 
MATCHED IN A SERIES FORMAT, ENT 
ER"CHR$<13>"THE RESISTIVE (OHMS 
) AND REACTIVE ELEMENTS (OHMS). 



139 



70 INPUT R*XC 

80 M=R*< 1+0* 2) 

90 N=SQR(M/RL-1) 

100 XL<1)=R*Q+XC 

110 XL<2)=RL*N 

120 XC=M/<Q+N) 

130 PRINT" XL1="XL<1) 'OHMS" 

140 PR I NT " XL2= " XL < 2 ) " OHMS " 

150 PRINT" XC=" XC'OHMB" 



140 



Sample Network D Problems 

Match the circuit shown in Fig. 4-20A by using network D. 
The load impedance is 50 ohms and the circuit Q is 10. 
SELECT A VALUE FOR Q 
?10 

ENTER THE LOAD IMPEDANCE 
?50 

WITH THE IMPEDANCE TO BE MATCHED IN A SERIES 
FORMAT, ENTER THE RESISTIVE (OHMS) AND REACTIVE 
(OHMS) ELEMENTS. 
?5,10 

XLl = 60OHMS 
XL2 = 150.831 OHMS 
XC = 38. 7966 OHMS 

Match the circuit shown in Fig. 4-20B to a load impedance of 
50 ohms using network D with a circuit Q of 10. 
SELECT A VALUE FOR Q 
?10 
ENTER THE LOAD IMPEDANCE 



NETWORK Q = 10 



10 



DEVICE TO BE 
MATCHED 




© 



NETWORK Q = 10 



10 



75 



if 



760£1 



613.39a 

;xci 



DEVICE TO BE 
MATCHED 



ft! 

340. 180 




Fig. 4-20. Network Q-10 f to be used in Network D problem. 
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?50 

WITH THE IMPEDANCE TO BE MATCHED IN A SERIES 

FORMAT, ENTER THE RESISTIVE (OHMS) AND REACTIVE 

ELEMENTS (OHMS). 

?75,10 

XL1 = 760 OHMS 

XL2 = 613.392 OHMS 

XC = 340. 177 OHMS 

This illustrates that, regardless of the value of the series real part, 

the T network is able to perform the necessary matching. 
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CAPACIT1VE-DIVIDER NETWORK 

The last network to be discussed is the capacitive-divider 
transformer. This network uses a resonant circuit to match 
between a low and high impedance. Whereas the minimum loss pad 
had considerable loss, this circuit exhibits only minimal loss. In 
addition, since a tuned circuit is used, it shows a band pass 
characteristic and is therefore useful in minimizing harmonic 
energy. Since specific frequencies are used for the center fre- 
quency and bandwidth, the solution provided by this program is a 
listing of real component values rather than reactances. 

To solve the capacitive-divider transformer manually, use the 
following equations and Fig. 4-21: 

q m Center Freq. 
Bandwidth 



CT = 



Q 



xlO 6 



2 ^f) 

C2/C1 = a/|^ _i 

CT= CI x C2 
C1 + C2 

Capacitive-Oivider Network Program 

5 'CAPACITIVE DIVIDER TRANSFORMER 

10 CLS 

20 PR I NT "ENTER LOAD IMPEDANCE <>50 

OHMS ) " 
30 INPUT RL 

40 PR I NT "ENTER CENTER FREQUENCY (M 
HZ)" 




Fig. 4-21. A C-D Transformer. 
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50 INPUT F 

60 PRINT" ENTER BANDWIDTH (MHZ)" 

70 INPUT B 

80 Q=F/B 

90 CT=Q*lE6/<2*3. 1416*F*< RL/2) ) 

100 M=S<3R(RL/50)-l 

110 C1=CT/(M/(M+1)> 

120 C2=M*C1 

130 L=i/<4*3. 1416t2*F*2*CT*lE~6) 

150 PRINT" Cl = " CI "PICOFARADS" 

160 PRINT" C2="C2" PICOFARADS" 

170 PRINT" l=" L"MICROHENRYS" 
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Example Capacitive-Divider Network Problem 

Using the capacitive-divider transformer network shown in 
Fig. 4-22, determine the required component values for a 
bandwidth of 250 kHz at a center frequency of 10 MHz. 
ENTER LOAD IMPEDANCE (>50 OHMS) 
?500 

ENTER CENTER FREQUENCY (MHZ) 
?10 

ENTER BANDWIDTH (MHZ) 
P.250 

CI = 3724. 15 PICOFARADS 
C2 = 8052.66 PICOFARADS 
L = .0994715 MICROHENRYS 
As proof, determine the resonant frequency of the tuned circuit by: 



f= 



1 



where C = series capacitance of CI and C2 

1 



f = 



V2tt .099xl0- 6 «2.55xl0- 9 



f = 



.159 



.159 



V2.52xl0- ls 1.59 xl0~ 8 



= 10.01 MHz 



C13724pF 
C28052pF 
L.099/iH 







=nC1 




\ o 


/%- 






k \ 


>RL500ft 


INPUT 
— 


50 S s 






i o 



Fig. 4-22. C-D transformer network to be used in example problem. 
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SERIES-MATCHING SECTION 

This next program uses a rather unusual matching technique 
known as the series-matching section. Because of the materials 
involved, it is more applicable to matching antennas rather than 
circuits. This technique uses a group of different impedance 
transmission lines matching different impedances. To understand, 
see Fig. 4-23; this diagram illustrates that two different impedance 
sections of transmission line are inserted between the system 
transmission line and the load to be matched. 

The program computes the lengths of the matching sections in 
electrical degrees. For practical implementation, the velocity 
factor constant will have to be applied to find the physical lengths. 
Using real transmission lines the different impedance values that 
can be matched are somewhat limited, however, there is enough 
range to match most real-world type of requirements. At higher 
frequencies where stripline and microstrip construction are preva- 
lent, the impedance of the transmission line can be tailored to fit 
the needs of the matching problem. 

The equations used to determine the length of the matching 
sections are: 



L2 = tan~ 




R(N~±)*(R- 1)2 



LI = tan" 1 2 ^A 



R + XNB-1 



where N = Zmatch 
Zline 



MAINLINE 


MATCHING SECTION 


L1 


_> 






L2 




3i ZL 


Z 
ANY LENGTH 


™i r~ 


Z 


? 


Z MATCH 



Fig, 4-23. A series-matching section illustrating two difference impedance 
sections of transmission line. 
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Y _XL 

o 



In some cases, the values of impedance for the matching section 
and main line section are too close together. This may cause the 
quantity under the radical to be negative, an invalid solution. The 
limits on the impedance of the transmission line can be calculated 
by the following formulas: 

Z ma tch >Z rn P V / SWRorZ <Z v / SWR 

match line T match line v 

FOR MATCHING TO OCCUR 

The program automatically calculates these limits and gives an 
error message when the values are too close. 

Series Section Program 



5 ' 


'SERIES MATCHING SECTION 




10 


CLS 




20 


PRINT"ENTER LOAD IMPEDANCE 


<R»X 


30 


) 

INPUT R»X 




40 


PRIIMT"ENTER THE IMPEDANCE i 
E MAIN LINE" 


OF TH 


50 


INPUT Z(l) 




60 


PRINT"ENTER THE IMPEDANCE ' 
E MATCHING SECTION" 


OF TH 


70 


INPUT Z(2) 




71 


S=SQR(R*2+Xf 2>/Z(l> 




72 


IF S<1 THEN S=l/S 




75 


GOTO 110 




80 


N=Z(2)/Z(1) 




90 


R=R/Z(1) 




100 X=X/Z(1) 




10S 


i GOTO 125 




110 IF Z(2)>Z(1)*SQR(B) OR ZC 


2)<Z( 




l)/SQR(S)GOTO 120 ELSE PRINT" I 




NSUFFICIENT RANGE FOR MATCHING I 




. CHOOSE A NEW VALUE FOR 


THE" 




CHR*< 13) "MATCHING SECTION 1 


11 
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115 
120 
125 
130 

140 
150 



160 
170 
180 
185 
190 



GOTO 20 

GOTO 80 

ON ERROR GOTO 190 

B=SO.R< ( <R-1H2+(X*2> >/<R#(N-l/ 

N>* 2-<R-l)*2--<X*2) ) ) 

A=< (N-R/N)*B+X)/(R+X*B*N-1> 

IF ATN(A)*57.29578<0 THEN L1=A 

TN < A) *57. 29578+180 ELBE L1=ATN 

<A)*57.29578 

L2=ATN(B)*57. 29578 

PR I NT "SECT I ON 1 

PRINT" SECTION 2 

END 

PRINT" INSUFFICIENT 

ATCHING. CHOOSE A 



OR"CHR*( 13) "THE 
ON. " sGOTO 60 



'LI" DEGREES" 
■L2" DEGREES" 

RANGE FOR M 
NEW VALUE F 
MATCHING SECTI 
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Series Matching Section Examples 

Match the antenna shown in Fig. 4-24 with a series matching 
section to a 50-ohm transmission line. The impedance of the 
matching section is 75 ohms. 
ENTER LOAD IMPEDANCE (R,X) 
? 250, 100 

ENTER THE IMPEDANCE OF THE MAIN LINE 
?50 

ENTER THE IMPEDANCE OF THE MATCHING SECTION 
?75 

INSUFFICIENT RANGE FOR MATCHING. CHOOSE A NEW 
VALUE FOR THE MATCHING SECTION 
ENTER LOAD IMPEDANCE (R,X) 
? 250, 100 

ENTER THE IMPEDANCE OF THE MATCHING SECTION 
?93 

INSUFFICIENT RANGE FOR MATCHING. CHOOSE A NEW 
VALUE FOR THE MATCHING SECTION. 
ENTER LOAD IMPEDANCE (R,X) 
? 250, 100 

ENTER THE IMPEDANCE OF THE MAIN LINE 
?50 

ENTER THE IMPEDANCE OF THE MATCHING SECTION 
?300 

SECTION 1 - 24.8774 DEGREES 
SECTION 2 - 20. 051 DEGREES 

This example shows that you would have to use a section of 
300-ohm balanced line to match the load. This might present some 
practical problems, interfacing between the 50-ohm unbalanced 
and 300-ohm balanced lines. However, it illustrates the wide range 
of antenna matching that can be achieved by using the series- 
matching section. 



MAINLINE 


MATCHING SECTION 

1 1 


^250,100 


30011 


Z=50O 


Z = 50fi 


L2 L1 
20,051° 24,8774° 



Fig. 4-24, Antenna to be used in solving the series-matching section 
problem. 
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MICR0STR1P 

The last program in this chapter is used to determine the 
width of microstrip transmission lines for specific values of 
impedance. As I mentioned in the previous section, the microstrip 
transmission lines can be tailored in impedance to fulfill the 
matching requirements. In many other cases, the impedance of the 
microstrip is also tailored to meet the matching requirements, i. e. , 
low noise microwave amplifiers or vhf/uhf power amplifiers. 

Figure 4-25 illustrates the relationships that are used in this 
program. The microstrip is that portion of the copper cladding that 
is not etched away. To compute the width of the microstrip, two 
different formulas are used. Once the width is determined, the 
correct answer is chosen according to the width-to-height ratio. 
Lines 60 through 90 perform the actual computations for the two 
formulas, while line 180 determines which answer of the two to 
use. 

The two formulas to calculate the width of the microstrip are: 



w/h < 2 w/h = 



8 e 



e 2A ~2 



w/h>2 w/h= .L(B-l-(n(2B-l)+ -f- 
7T Z E 

(ln(B-l) + .39 




CIRCUIT BOARD 
THICKNESS =T 



COPPER CLADDING 
1 OUNCE 



Fig, 4-25, The relationships used in the microstrip program. 
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Z / E +1 E-l -11 \ 

whereA=^V — + — (.23+-=- ) 

60 2 E+l E ) 



377 it 
B = 



2Z VE 

o r 



Microstrip Program 

5 'MICROSTRIP CALCULATIONS 

10 CLS 

20 PR I NT "ENTER DESIRED IMPEDANCE 

F MICROSTRIP" 
30 INPUT Z 
40 PR I NT "ENTER DIELECTRIC CONSTANT 

OF CIRCUIT BOARD MATERIAL" 
50 INPUT E 
55 GOTO 100 
60 A=(Z/60)*SQR< (E+l)/2)+( <E-1)/(E 

+1) >*<.ll/E+.23> 
70 B=592. 1902/<Z*SQR<E> ) 
80 WH(1)=8*2.7183*A/(2.7183*<2*A)- 

2) 
90 WH ( 2 ) = . 6366* ( < B- 1 ) -LOG ( 2*B-1 ) + ( 

<E-1)/<2*E) )*<L0G<B-l>+.39-.61/ 

E)) 
95 GOTO 140 
100 PRINT'ENTER THICKNESS OF CIRCU 

IT BOARD (MILS) 
110 INPUT T 
120 PR I NT "ENTER THICKNESS OF COPPE 

R (OUNCES) 
130 INPUT 
135 GOTO 60 
140 0=0*1.4 
150 H=T-2*0 
160 R=WH(1)*H 
170 S=WH(2)*H 

180 IF WH(1)<2 PRINT"THE MICROSTRI 
P IS"R"MILS WIDE" ELSE PRINT"T 
HE MICROSTRIP IS"S"MILS WIDE" 
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Microstrip Examples 

Determine the width of a 50-ohm line on a glass epoxy board 
that is 1/16-inch thick. The copper clad is 1 ounce. 
ENTER DESIRED IMPEDANCE OF MICROSTRIP 
?50 

ENTER DIELECTRIC CONSTANT OF CIRCUIT BOARD 
MATERIAL 
? 5 (Glass epoxy) 

ENTER THICKNESS OF CIRCUIT BOARD (MILS) 
?62 

ENTER THICKNESS OF COPPER (OUNCES) 
?1 

THE MICROSTRIP IS 102.721 MILS WIDE. 
Next, compute the width of the same microstrip except on an 
alumina substrate. 

ENTER DESIRED IMPEDANCE OF MICROSTRIP 
?50 

ENTER DIELECTRIC CONSTANT OF CIRCUIT BOARD 
MATERIAL 
?8 8 

ENTER THICKNESS OF CIRCUIT BOARD (MILS) 
?62 

ENTER THICKNESS OF COPPER (OUNCES) 
?1 
THE MICROSTRIP IS 63.8949 MILS WIDE. 

This illustrates that as the dielectric constant of the circuit 
board material is raised, the lines can be of smaller width. By using 
substrates with a higher dielectric constant, microwave circuitry 
can be reduced significantly in size as compared to lower frequency 
type of circuits. 
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Chapter 5 
Geometry 




This last chapter is not a return to the eighth grade. However, it is 
devoted to the person who uses his computer for the more exotic 
problems. The first group of programs in this chapter will deal with 
the hyperbolic functions which are used in a lot of engineering 
problems. In addition, hyperbolic functions are quite useful in 
connection with differential equations. 
HYPERBOLIC FUNCTIONS 

The basic identities of the hyperboic sine and cosine are 
related to the coordinates of the points X, Y on the basic hyperbola 
X 2 - Y 2 = 1. Sinh and Cosh, the hyperboic sine and cosine, 
respectively, are identified as follows: 

SinhX-M>(e x -e- x ) 
CoshX = V2(e x + e" x ) 

The remaining hyperbolic functions are defined in terms of sinh x 
and cosh x as follows: 

sinhX 



tanhX = 
cothX = 
sechX = 
cschX = 



coshX 

coshX 
sinhX 

1 

coshX 

1 
sinhX 
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The graphs of the hyperbolic functions are shown in Figs. 5-1 
through 5-6. The hyperbolic cosine is an even function, cosh(- x) = 
cosh x, and the hyperbolic sine is an odd function, sinh(-x) = 
- sinh x. Therefore, one curve is symmetrical about the y-axis and 
the other is symmetrical with respect to the origin. Here, the 
hyperbolic functions behave like ordinary trigonometric functions. 
The hyperbolic functions are so important that their numerical 
values have been calculated and tabulated just as have the ordinary 
trigonometric functions. However, it is very inconvenient to have 
to stop work and look up a value in a table. Thus, the reason for the 
following she programs. In this way they can be incorporated 
directly into a program. As a good example, the impedance of the 
transmission line in the Gamma and Omega programs was 
calculated using the hyperbolic cosine. 
Hyperbolic Sine Program 

5 'HYPERBOLIC SINE 

10 CLS 

20 PRINT H ENTER X FOR HYPERBOLIC SI 

NE (SINH) OF X" 
30 INPUT X 

40 A=<<EXP<l)*X>-<EXPU>*-X)>/2 
50 PRINT"THE HYPERBOLIC SINE OF"X B 

IS" A 
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I 1 H 



2000 T 



1600 



1200- 



800- 



400-- 



-H 



-10-8 »6/<-4 ~-2 



-400. 



-800-* 



-1200- 



-1600. 



-2000 



2 4 6 8 10 



Fig. 5-1. A graph illustrating the hyperbolic sine. 
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Hyperbolic Sine Examples 

ENTER X FOR THE HYPERBOLIC SINE (SINH) OF X 

?9 

THE HYPERBOLIC SINE OF 9 IS 4051.54 

ENTER X FOR THE HYPERBOLIC SINE (SINH) OF X 

?.6 

THEHYPERBOLICSINEOF.6 IS .636654 

Note from Fig. 5-1 that with small values of X, the sinh is 

correspondingly small, greater values of X produce infinitely 

larger numbers. 
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Hyperbolic Cosine Program 

5 'HYPERBOLIC COSINE 

10 CLS 

20 PRINT"ENTER X FOR HYPERBOLIC CO 

SINE (COSH) OF X" 
30 INPUT X 

40 A==< (EXPU )♦*) + ( EXP(1)*~X) )/2 
50 PRINT" THE HYPERBOLIC COSINE OF" 

X"IS"A 
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Fig. 5-2. A graph illustrating the hyperbolic cosine. 
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Hyperbolic Cosine Examples 

ENTER X FOR HYPERBOLIC COSINE (COSH) OF X 

?-3 

THEHYPERBOLICCOSINEOF-3 IS 10.0677 

ENTER X FOR HYPERBOLIC COSINE (COSH) OF X 

?0 

THE HYPERBOIC COSINE OF IS 1 

You'll note that this example shows that the minimum value 
the hyperbolic cosine can assume is 1. Figure 5-2 also illustrates 
this fact, showing that the curve is symmetrical about the Y axis. 
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Hyperbolic Tangent Program 

5 'HYPERBOLIC TANGENT 

10 CLS 

20 PRINT" ENTER X FOR HYPERBOLIC TA 

NGENT (TANH) OF X" 
30 INPUT X 
40 A=MEXP<1)* X)-(EXP<l>t-X> )/( (EX 

P(l)» ...X). + <EX.P..<.1>* -X>> 
50 PRINT"THE HYPERBOLIC TANGENT OF 

H X"IS"A 
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Fig. 5-3. A graph illustrating the hyperbolic tangent 

Hyperbolic Tangent Examples 

ENTERXFORHYPERBOLICTANGENT(TANH) OFX 
?2 

THE HYPERBOLIC TANGENT OF 2 IS .964028 
As you run a hyperbolic tangent you'll find that it quickly 
approaches, but never reaches, 1. The TRS-80 is limited in the 
number of digits it can display. Even using double precision will 
only allow you to compute a hyperbolic tangent of about 8.3. To use 
double precision, change lines A in lines 40 and 50 to A#. 
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Hyperbolic Cotangent Program 

5 'HYPERBOLIC COTANGENT 
10 CLS 
20 PRINT "ENTER X FOR HYPERBOLIC CO 

TANGENT (COTH) OF X" 
30 INPUT X 
40 A=((EXP(1)» X) + (EXP(l)»--X>>/< (EX 

P<mX)-<EXP<l)»-X) ) 
50 PRINT" THE HYPERBOLIC COTANGENT 

OF" X" IS" A 
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Hyperbolic Cotangent Examples 

ENTER X FOR HYPERBOLIC COTANGENT (COTH) OF X 

?IE-7 

THE HYPERBOLIC COTANGENT OF 1E-07 IS 5.59241E+06 

ENTER X FOR HYPERBOLIC COTANGENT (COTH) OF X 

?5 

THE HYPERBOLIC COTANGENT OF 5 IS 1. 00009 

From these examples and Fig. 5-4 you can see that as the 
hyperbolic cotangent approaches zero, its value becomes infinitely 
big, and it approaches 1 as X becomes large. 



2-- 



I 1 1 i I 1 i H 

-1.6-1.4-1.2-1.0 -.8 -.6 -.4 -2 




H 1 1 1 1 1 1- 1 

.2 A £ .8 IX) 12 14 1.6 



Fig. 5-4. A graph illustrating the hyperbolic cotangent. 
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Hyperbolic Secant Program 

5 'HYPERBOLIC SECANT 

10 CLS 

20 PRINT"ENTER X FOR HYPERBOLIC SE 

CANT (SECH) OF X" 
30 INPUT X 

40 A=2/( <EXPQUX) + <EXP(1>*-X> ) 
50 PRINT "THE HYPERBOLIC SECANT OF" 

X"IS"A 
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Fig. 5-6. A graph illustrating the hyperbolic secant 



Hyperbolic Secant Examples 

ENTERXFORHYPERBOLICSECANT(SECH) OFX 

?0 

THE HYPERBOLIC SECANT OF IS 1 

ENTER X FOR HYPERBOLIC SECANT (SECH) OF X 

?50 

THE HYPERBOLIC SECANT OF 50 IS 3.8575E-22 

The maximum value for the hyperbolic secant is one. As X 
gets increasingly big, either positive or negative, the hyperbolic 
secant approaches zero. 

Generally, the hyperbolic functions behave much like e x /2, 
£ x /2, unity, or zero as shown in the following table: 



X (Large and Positive) 
CoshX^sinhx 

~y 2 e x 
tanhX^cothX^l 
sechX^csch^ 

~2e~ x ~0 



X Negative, |X| Large 
CoshX^-sinhX 

tanhX^coth^-1 
sechX^-cschX 
^2e x -*0 
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Hyperbolic Cosecant Program 

5 'HYPERBOLIC COSECANT 

10 CLS 

20 PR INT" ENTER X FOR HYPERBOLIC CO 

SECANT (CSCH) OF X" 
30 IMPUT X 

40 A=2/( <EXP<1)»X)-<EXPU)*-X)) 
50 PRINT" THE HYPERBOLIC COSECANT 

F"X"IS"A 
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.8 1.0 



Fig. 5-5. A graph illustrating the hyperbolic cosecant. 

Hyperbolic Cosecant Examples 

ENTER X FOR HYPERBOLIC COSECANT (CSCH) OF X 

?2 

THE HYPERBOLIC COSECANT OF 2 IS .275721 

ENTERX FOR THE HYPERBOLIC COSECANT (CSCH) OFX 

?-2 

THE HYPERBOLIC COSECANT OF - 2 IS - .275721 

The hyperbolic cosecant is very similar to the hyperbolic 
cotangent, except that instead of being asymptotic to Y=l it is 
asymptotic to Y=0. Again, as the value of X gets smaller, the 
hyperbolic cosecant becomes infinitely big. 
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INVERSE HYPERBOLIC FUNCTIONS 

The inverse hyperbolic functions perform the inverse opera- 
tion of the hyperbolic functions. That is, if you compute X = sinh Y, 
then the inverse hyperbolic secant is Y = sinh™ 1 X. The general 
curves for the inverse hyperbolic functions are shown in Fig. 5-7 
through 5-12. 

To compute the inverse hyperbolic functions, use the follow- 
ing formulas: 

sinh^x^lnCx + Cx 2 *!)' 5 ) 

] + x 
tantr l x = V4(n( - v 



sech-x = l„(-7 + (-?-l)') 



cosh™ 1 x 2 = secbr * — 

A 



cothr 1 x 2 = tcnlr'T 
csdr 1 x 2 = sinh -1 — 

Inverse Hyperbolic Sine Program 

5 MNV. HYPERBOLIC SINE 

10 CLS 

20 PRINT" ENTER X FOR INVERSE HYPER 

BOLIC SINE OF X" 
30 INPUT X 

40 A=LOG<X+< (X»2>+1)».5> 
50 PRINT"THE INVERSE HYPERBOLIC SI 

NE OF" X" IS" A 
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Fig. 5-7. A graph illustrating the inverse hyperbolic sine. 



Inverse Hyperbolic Sine Examples 

ENTER X FOR INVERSE HYPERBOLIC SINE OF X 

?5 

THE INVERSE HYPERBOLIC SINE OF 5 IS 2.3124 

ENTERXFORINVERSEHYPERBOLICSINE OFX 

?-5 

THE INVERSE HYPERBOLIC SINE OF -5 IS -2.3124 

These examples show that the inverse hyperbolic sine is a 
single value; that is, as X varies, Y produces a different answer for 
each value of X. The curve does not reverse to produce the same 
value for X as for - X. 
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Inverse Hyperbolic Cosine Program 

5 * INV. HYPERBOLIC COSINE 



10 


CLS 


20 


PRINT" ENTER X FOR INVERSE HYi 




BOLIC COSINE OF X" 


30 


INPUT X 


40 


A=LOG(X+( <XC2)-1H .5) 


50 


PR I NT "THE INVERSE HYPERBOLIC 




SINE OF" X" IS" A 45B=~A 



CO 



170 
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Fig. 5-8. A graph illustrating the inverse hyperbolic cosine. 

Inverse Hyperbolic Cosine Examples 

ENTER X FOR INVERSE HYPERBOLIC COSINE OF X 

?.5 

?FC ERROR IN 40 

Whereas cosh 0=1, the inverse hyperbolic function is limited 

to the range of greater than one. 

ENTER X FOR INVERSE HYPERBOLIC COSINE OF X 

?20 

THE INVERSE HYPERBOLIC COSINE OF 20 IS 3.68825 AND 

-3.68825 

Note that in the solution for the inverse hyperbolic cosine, for each 

value of X there are two corresponding values for Y. 
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Inverse Hyperbolic Tangent Program 



5=. •) 


INV. HYPERBOLIC TANGENT 


10 


CLS 


20 


PR I NT "ENTER X FOR INVERSE HY 




BOLIC TANGENT OF X" 


30 


INPUT X 


40 


A=.5*LOG< <1+X)/<1-X) ) 


50 


PRINT"THE INVERSE HYPERBOLIC 




NGENT OF"X H IS"A 



TA 
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Fig. 5-9. A graph illustrating the inverse hyperbolic tangent. 

Inverse Hyperbolic Tangent Examples 

ENTER X FOR INVERSE HYPERBOLIC TANGENT OF X 

?.4 

THE INVERSE HYPERBOLIC TANGENT OF .4 IS .423649 

ENTER X FOR INVERSE HYPERBOLIC TANGENT OF X 

?.l 

THE INVERSE HYPERBOLIC TANGENT OF. 1 IS .100335 

You'll note from these two examples and the curve shown in 
Fig. 5-9 that the inverse hyperbolic tangent is almost a straight-line 
function. However, as X approaches 1, the function becomes 
asymptotic to at the X = 1 and X = - 1 lines. 
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Inverse Hyperbolic Cotangent Program 

5 ' INV. HYPERBOLIC COTANGENT 

10 CLS 

20 PRINT"ENTER X FOR INVERSE HYPER 

BOLIC COTANGENT OF X" 
30 INPUT X 

40 A=.5*LOG<(l+l/X>/<l-l/X>) 
50 PRINT'THE INVERSE HYPERBOLIC CO 

TANGENT OF" X" IS" A 
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Fig. 5-10. A graph illustrating the inverse hyperbolic cotangent. 

Inverse Hyperbolic Cotangent Examples 

ENTER X FOR INVERSE HYPERBOLIC COTANGENT OF X 

?.5 

?FC ERROR IN 40 

This shows that the function Y = coth" l X is only valid for X greater 

than 1 and X less than - 1. The region running from - 1 through to 

1 is invalid. 

ENTER X FOR INVERSE HYPERBOLIC COTANGENT OF X 

?-6 

THE INVERSE HYPERBOLIC COTANGENT OF -6 IS 

-.168236 
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Inverse Hyperbolic Secant Program 

5 MNV. HYPERBOLIC SECANT 

10 CLS 

20 PR I NT "ENTER X FOR INVERSE HYPER 

BOLIC SECANT OF X" 
30 INPUT X 
40 A=L0G<1/X-K<<1/X+2>)~1>*.5) 

45 B=-A 

50 PRINT"THE INVERSE HYPERBOLIC SE 
CANT OF" X" IS" A" AND" B 
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Fig. 5-11. A graph illustrating the inverse hyperbolic secant. 

Inverse Hyperbolic Secant Examples 

ENTER X FOR INVERSE HYPERBOLIC SECANT OF X 

?1 

THE INVERSE HYPERBOLIC SECANT OF 1 IS 

ENTER X FOR INVERSE HYPERBOLIC SECANT OF X 

?.6 

THE INVERSE HYPERBOLIC SECANT OF .6 IS 1.09861 AND 

-1.09861 

The inverse hyperbolic secant acts in the same manner as the 
inverse hyperbolic cosine. Each unique value of X produces two 
answers, Y and -Y. Small values of X eventually become 
asymptotic to the Y axis, both plus and minus (see Fig. 5-11). 
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Inverse Hyperbolic Cosecant Program 



5 ' 


' INV. HYPERBOLIC COSECANT 


10 


CLS 


20 


PRINT'*ENTER X FOR INVERSE HYPER 




BOLIC COSECANT OF X" 


30 


INPUT X 


40 


A=LOGU/X+<( <1/X*2)>+1>*.5) 


50 


PRINT" THE INVERSE HYPERBOLIC CO 




SECANT OF" X" IS" A 
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Fig. 5-12. A graph illustrating the inverse hyperbolic cosecant. 

Inverse Hyperbolic Cosecant Examples 

ENTER X FORINVERSE HYPERBOLIC COSECANT OF X 

?0 

?/0 ERROR IN 40 

ENTER X FOR INVERSE HYPERBOLIC COSECANT OF X 

?8 

THE INVERSE HYPERBOLIC COSECANT OF 8 IS . 124677 

The curve of the inverse hyperbolic cosecant is not unlike the 
curve of the inverse hyperbolic cotangent. However, instead of 
small values of X becoming asymptotic to X = 1, the curve is 
asymptotic to the Y axis. Large values of X eventually become 
asymptotic to the X axis (see Fig. 5-12). 
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COMPLEX FUNCTIONS 

The next twelve programs deal with the complex trigonomet- 
ric functions, or the sine, cosine, ... of complex numbers and the 
inverse trigonometric functions. 

The equations to solve each of the following programs are: 

sin (x 4- £y) = sin x cosh Y + 1 cos x sinhY 
cos (x + tY) = cos x cosh Y - 1 sin x sinhY 

* n , ,. v sin2x-htsinh2Y 
** (x + t y )= cos 2x + cosh2Y 

cotCx + ly)^ 8in ?- t8inh ^ 
cosh2x-cos2Y 

1 



esc (x + iy) = 
sec(x + iy) = 



sin(x + iY) 

1 
cos (x + tY) 



sin -1 (x + iy) = sin -1 B + isgn(y) 1 n(A + V A 2 -l) 
where 



A*=fcVCX + l) 2 + y 2 + ftV(c-l) 8 + y 2 



B = VzViX + lf + f-VzV (c- I) 2 + f 

sgn(Y) = lify^O 
-lify 



cos" l (x + iy) = cos- ' B + i sgn(d) In (A + VA 2 -1) 
where 



A = y2V(x + l) 2 + y 2 + V2V(c-l) 2 + y 2 



B = V2 V(x + l) z + y 2 ~ V 2 V(c-lf + f 
sgn(Y) = lify 2=0 
-lify 
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tan- 1 (x + iy) = V 2 (7r-tarr , l±I-tan- 1 ^- + j- 

ln / (l+y) 2 +x 2 \ 
111 \Q-y) 2 +x 2 / / 

cot" ' (x + iy) = -77- - tan - ' (x + iy) 

csc- ' (x + iy) = sin" ' ( x + £y)J 

sec- l (x + iy) = cos-' L^y J 
Complex Sine Program 

5 'COMPLEX SINE 

10 CLS 

20 PRINT"ENTER COMPLEX NUNBER ( REA 

L, IMAGINARY)" 
30 INPUT R»I 
40 AR=SIN<R)*<<EXPU>* I+EXP(l)t-I) 

/2) 
50 AI=C0S(R)*(<EXP<1H I-EXP(l)t-I) 

/2) 
60 PRINT"THE COMPLEX SINE IS"AR;AI 

II T II 
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Complex SineExamples 

ENTER COMPLEX NUMBER (REAL, IMAGINARY) 

?2 —3 

THE COMPLEX SINE IS 9.1545 4.16891 

ENTER COMPLEX NUMBER (REAL, IMAGINARY) 

?3,-4 

THE COMPLEX SINE IS 3.85374 27.0168 
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Complex Cosine Program 

5 'COMPLEX COSINE 

10 CLS 

20 PRINT" ENTER COMPLEX NUNPFR ( RFA 

Li IMAGINARY)" 
30 INPUT R,I 
40 AR=COS(R)*< <EXP(DM+EXP<D* -D 

/2) 

50 AI=-(SIN(R)*( (EXP(l)iI-EXP( DA- 
D/2)) 

60 PRINT" THE COMPLEX COSINE 13" AR; 
AI" I" 
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Complex Cosine Examples 

ENTER COMPLEX NUMBER (REAL, IMAGINARY) 

?2 —3 

THE COMPLEX COSINE IS 4. 18963 9. 10923 

ENTER COMPLEX NUMBER (REAL, IMAGINARY) 

?-.2,5 

THE COMPLEX COSINE IS 72. 7307 14. 7419 
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Complex Tangent Program 

5 ' COMPLEX TANGENT 

10 CLS 

20 PRINT-ENTER COMPLEX NUNBER < REA 

Li IMAGINARY)" 
30 INPUT R, I 
40 D=COS<2*R)-K <EXP<1>* <2*I)+EXP(1 

]♦ <2*-I) )/2) 
50 AR=SIN<2*R)/<D) 
60 AI=((EXP<l)t <2*I)+EXP<1)* <2*-I> 

)/2)/(D) 
70 PRINT"THE COMPLEX TANGENT IB"AR 

?AI"i" 
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Complex Tangent Examples 

ENTER COMPLEX NUMBER (REAL, IMAGINARY) 
?-2,— 3 

THE COMPLEX TANGENT IS 3.76403E-03 1.00325 
ENTER COMPLEX NUMBER (REAL, IMAGINARY) 

?10,-5 

THE COMPLEX TANGENT IS 8. 28923E-05 .999963 
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Complex Cotangent Program 

5 'COMPLEX COTANGENT 

10 CLS 

20 PRINT" ENTER COMPLEX NUNBER ( REA 

L> IMAGINARY)" 
30 INPUT R»I 
40 D=-COS<2*R) + < (EXP(l)t <2*I)+EXP< 

1H <2*-I))/2) 
50 AR=SIN<2*R)/<D> 
60 AI=-((EXP(1H <2*I)-EXP<1H <2*-I 

>>/2>/<D) 
70 PRINT" THE COMPLEX COTANGENT IS" 

AR5AI" i" 
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Complex Cotangent Examples 

ENTER COMPLEX NUMBER (REAL, IMAGINARY) 

?2 3 

THE COMPLEX COTANGENT IS 3. 73971E-03 -. 996758 

ENTER COMPLEX NUMBER (REAL, IMAGINARY) 

?~12,.3 

THE COMPLEX CONTANGENT IS 1.18954 - .836285 
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Complex Secant Program 

5 'COMPLEX SECANT 
10 CLS 
20 PRINT "ENTER COMPLEX NUMBER ( REA 

L» IMAGINARY)" 
30 INPUT R, I 
40 AR=COS<R>*< (EXP(l)t I+EXP(l)f-I) 

/2) 
50 AI=-<SIN<R)*(<EXP<1)4 I-EXPUH- 

I)/2) ) 
60 M=SG>R(ARC2+AI*2) 
70 T=ATN(AI/AR) 
80 FR=l/M*<COS(-T> ) 
90 FI=1/M*<SIN<-T) ) 

100 PRINT"THE COMPLEX SECANT IS"FR 
?FI 
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Complex Secant Examples 

ENTER COMPLEX NUMBER (REAL, IMAGINARY) 

?4,-9 

THE COMPLEX SECANT IS 1.61332E-04 - 1.86794E-04 

ENTER COMPLEX NUMBER (REAL, IMAGINARY) 

? -.009, 10 

THE COMPLEX SECANT IS 9. 07961 E-05 - 8. 172O1E-07 
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Complex Cosecant Program 

5 'COMPLEX COSECANT 

10 CLS 

20 PRINT-ENTER COMPLEX NUNBER ( REA 

L, IMAGINARY)" 
30 INPUT R,I 
40 AR=SIN(R)*( <EXPU)M+EXP<1)*-I> 

/2) 
50 AI=COS<R>#( (EXPU)t I-EXP(1)+-I) 

/2) 
60 M=SQR(ARC2+AJt2) 
70 T=ATN(AI/AR) - 
80 FR=l/M*(COS(~T) ) 
90 FI = 1/M*<SIIM<-T> ) 

100 PRINT "THE COMPLEX COSECANT IS" 
FR;FI 
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Complex Cosecant Examples 

ENTER COMPLEX NUMBER (REAL, IMAGINARY) 

?3,-.5 

THE COMPLEX COSECANT IS .545987-1. 77002 

?-.3, 10 

THE COMPLEX COSECANTIS 2.68332E-05 8.67444E-05 
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Complex Arc Sine Program 

5 'COMPLEX ARC SINE 

10 CLS 

20 PRINT" ENTER COMPLEX NUMBER ( REA 

Li IMAGINARY) 
30 INPUT R» I 
40 M=SQR< <R+l)*2+(It 2) ) 
50 N=SQR< <R-l)t2+Ut2) > 
60 Z(l )-.5*M+.5*N 
70 Z<2)=.5*-M+.5*N 
80 Z(3)=.5*M+.5*-N 
90 I i 4 ) - . 5 *-M+ . 5 * -N 
100 Z(5)=.5*M-.5*N 
110 Z(6)=.5*-M-.5*N 
120 Z(7)=.5*M-.5*-N 

130 Z<8)=.5* M-.5*-N 

140 FOR P=5 TO S 

145 IF Z<P)<0 GOTO 160 

146 IF Z<P)>1 GOTO 160 

1 50 AR < P > =ATN < Z ( P > /SQR ( -Z < P ) *Z < P ) + 

1) ) 
160 NEXT P 

170 IF I=>0 THEN S=l i STOPS=- 1 
1S0 FOR P=l TO 4 
185 IF ZiPKl GOTO 210 
190 Al(P)=<LOG<Z<P)+SQR<Z(P)t 2-1) ) 

>*S 
200 A2 ( P ) = ( LOG < Z ( P ) -SOR < Z < P > t 2- 1 > ) 

)*S 
210 NEXT P 
220 PRINT" THE COMPLEX ARC SINE IS 



230 


FOR P=5 TO 


8 


240 


IF AR(P)=0 


GOTO300 


250 


FOR 6=1 TO 


4 


260 


IF A1(G)=0 


GOTO290 


270 


PRINTAR(P)! 


.AUG) 


2Q0 


PRINTAR(P)! 


»A2(G) 


290 


NEXT G 




300 


NEXT P 
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Complex Arc Sine Examples 

ENTER COMPLEX NUMBER (REAL, IMAGINARY) 

?5, 8 

THE COMPLEX ARC SINE IS - 

.556083 2.93872 

.556083 -2.93872 

Note that all the inverse complex trigonometric functions are 

double-valued. 

ENTER COMPLEX NUMBER (REAL, IMAGINARY) 

?l,-2 

THE COMPLEX ARC SINE IS - 

.4247079 -1.52857 

.427079 1.52857 
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Complex Arc Cosine Program 

5 'COMPLEX ARC COSINE 

10 CLS 

20 PR I NT "ENTER COMPLEX NUMBER < REA 

Li IMAGINARY) 
30 INPUT R,I 
40 M=SQR< <R+l)t 2+<It2) ) 
50 N=SQR< <R-l) + 2+<It2)) 
60 Z(1)=.5*M+.5*N 
70 Z<2)=.5*-M+.5*N 
80 Z<3)=.5*M+.5*-N 
90 Z<4)=.5*-M+.5*-N 
100 Z(5)=.5*M-.5*N 
110 Z<6) = .5*-M--.5*N 
120 Z(7)=.5*M-.5*-N 
130 Z(8)=.5*~M~.5*~N 
140 FOR P=5 TO 8 

145 IF Z(P)<0 GOTO 160 

146 IF Z<P)>1 GOTO 160 

150 AR < P ) =-ATN < Z ( P ) /SQR < -Z ( P ) *Z ( P ) 
+1) )+1.5708 

160 NEXT P 

170 IF I=>0 THEN S=l ELBE S=-l 

180 FOR P=l TO 4 

185 IF Z(P)<1 GOTO 210 

1 90 A 1 ( P ) =- ( LOG < Z ( P ) +SO.R ( Z < P ) t 2- 1 ) 

) )*S 
200 A2 ( P ) =- ( LOG ( Z ( P ) -SO.R < Z < P >♦ 2-1 > 

) )*S 
210 NEXT P 
220 PRINT" THE COMPLEX ARC COSINE I 

S — " 
230 FOR P=5 TO 8 
240 IF AR(P)=0 GOTO300 
250 FOR G=l TO 4 
260 IF A1(G)=0 GOTO290 
270 PRINTAR(P),A1(G) 
280 PRINTAR(P),A2(G) 
290 NEXT G 
300 NEXT P 
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Complex Arc Cosine Examples 

ENTER COMPLEX NUMBER (REAL, IMAGINARY) 

?l,-2 

THE COMPLEX ARC COSINE IS - 

1.14372 1.52857 

1.14872 -1,52857 

ENTER COMPLEX NUMBER (REAL, IMAGINARY) 

?23,8 

THE COMPLEX ARC COSINE IS - 

.334995 -3.88541 

.334995 3.88556 
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Complex Arc Tangent Program 

5 'COMPLEX ARC COTANGENT 

10 CLS 

20 PRINT-ENTER COMPLEX NUMBER ( REA 

L, IMAGINARY) 
30 INPUT R,I 

40 AR=.5*<3.141S9-ATN( <1+I)/R)-ATN 
<<1-I)/R)) 

30 AI=LOG<<<i + I)C2+RC2)/< <l-I)+2+R 

♦ 2> >/4 
60 PRINT "THE COMPLEX ARC COTANGENT 

IS 
70 PRINT3. 14159/2-AR ? AI 
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Complex Arc Tangent Examples 

ENTER COMPLEX NUMBER (REAL, IMAGINARY) 

?5, 8 

THE COMPLEX ARC TANGENT IS - 

1.51422 .0898434 

1.51422 -.0898434 



198 



Complex Arc Cotangent Program 

5 'COMPLEX ARC COTANGENT 

10 CLS 

20 PRINT"ENTER COMPLEX NUMBER ( REA 

L» IMAGINARY) 
30 INPUT R,I 
40 AR=.5*<3. 14159- ATN( <1+I)/R)-ATN 

< (1-D/R) ) 
50 AI=LOG< <<i+I)C2+RC2>/(<l-I>*2+R 

♦2) )/4 
60 PRINT" THE COMPLEX ARC COTANGENT- 
IS --» 
70 PRINT3. 14159/2~AR,AI 
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Complex Arc Cotangent Examples 

ENTER COMPLEX NUMBER (REAL, IMAGINARY) 

?2 3 

THE COMPLEX ARC COTANGENT IS - 

.160875 -.229073 

.160875 .229073 

ENTER COMPLEX NUMBER (REAL, IMAGINARY) 

?U 

THE COMPLEX ARC COTANGENT IS - 

.553574 .40236 

.553574 -.40236 
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Complex Arc Secant Program 



5 ' 


COMPl EX ARC SECANT 






10 


CLS 






20 


PRINT" ENTER COMPLEX NUMBER ( 


REAL, 




IMAGINARY) 






30 


INPUT R> I 






35 ' 


GOSUB 310 






40 1 


M=SQR< <R+l>r.2+< 1 12) ) 






50 1 


M=SQR( (R-1)[2+(K2) ) 






60 


Z(1)=.5*M+.5*N 






70 ; 


Z<2)=.5*-M+.5*N 






80 


Z(3)=.5*M+.5*-N 






90 


Z(4> = .5*—M+.5*-N 






100 


Z(5) = .5*M--.5*N 






110 


Z(6)=.5*-M~.5*N 






120 


Z<7)=.5*M-.5*-N 






130 


Z(8)=.5*-M-.5*-N 






140 


FOR P=5 TO 8 






145 


IF Z(P)<0 GOTO 160 






146 


IF Z<P)>1 GOTO 160 






150 


AR(P)= ATN(Z<P)/SQR(~ 

<P)+1) )+1.5708 


Z(P)*Z 




160 


NEXT P 






170 


IF I=>0 THEN S-l ELSE 


S=-l 




180 


FOR P=l TO 4 






185 


IF ZIPX1 GOTO 210 






190 


AKP)=-(LOG(Z(P)+SQR< 


Z(P)*2- 


-1) ) )*S 


200 


A2 ( P ) =- ( LOG < Z < P ) -SQR ( Z ( P ) ♦ 2- 


-1) ))*S 


210 


NEXT P 






220 


PR I NT "THE COMPLEX ARC 


SECANT IS -- 


230 


FOR P=5 TO 8 






240 


IF AR<P>=0 GOT 0300 






250 


FOR 6=1 TO 4 






260 


IF A1(G)=0 GOTO290 






270 


PRINTAR(P) ? A1(G) 






280 


PRINTAR(P)»A2<G) 






290 


NEXT G 






300 


NEXT P 






305 


END 






310 


F=SQR<RC2+IL2) 







201 



320 


T=ATN( I/R) 


330 


F=l/F 


340 


T=-T 


350 


R=F*COS(T) 


360 


I=F*SIIM(T> 


370 


RETURN 
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Complex Arc Secant Examples 

ENTER COMPLEX NUMBER (REAL, IMAGINARY) 

?8,-5 

THE COMPLEX ARC SECANT IS - 

1.48093 -.0563788 

1.48093 .0563784 

ENTER COMPLEX NUMBER (REAL, IMAGINARY) 

?-.03 .03 

THE COMPLEX ARC SECANTIS - 

.785853 -3.85313 

.785853 3.85321 
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Complex Arc Cosecant 

5 'COMPLEX ARC COSECANT 

10 CLS 

20 PRINT" ENTER COMPLEX NUMBER < REA 

Li IMAGINARY) 
30 INPUT R» I 
35 GOSUB 310 
40 M=SQR ( ( R+ 1 ) ♦ 2+ ( I * 2 ) ) 
50 N=SO.R( <R-1)*2+<I*2> ) 
60 Z( 1 >=.5*M+.5*N 
70 Z<2)==.5*~M+.5*N 
B0 Z(3)=.5*M+.5*—N 
90 Z(4)=.5*-M+.5*-N 
100 Z<5)=.5*M-.5*N 
1 10 Z(6)=.5*-M-.5*N 
120 Z<7)=.5*M-.5*-N 
130 Z<8)=.5*-M-.5*-N 
140 FOR P=5 TO 8 

145 IF Z(P)<0 GOTO 160 

146 IF Z(P)>1 GOTO 160 

150 AR < P ) =ATN ( Z ( P ) /SQR < Z ( P ) *'Z ( P ) t 

1) ) 
160 NEXT P 

170 IF I = >0 THEN S=l El SE S=--l 

180 FOR P=l TO 4 

185 IF Z(PXi GOTO 210 

190 A 1 ( P ) = ( LOG ( Z ( P ) +SO R ( Z ( P ) * 2 - 1 ) ) 

)*S 
200 A2 ( P ) = ( LOG ( Z ( P ) -SQR ( Z ( P ) * 2- 1 ) > 

)*S 
210 NEXT P 

220 PRINT "THE COMF1 EX ARC COSECANT 

IS - — " 
230 FOR P=5 TO 8 
240 IF AR(P)=0 GOTO300 
250 FOR 6=1 TO 4 
260 IF A1(G)=0 GOTO290 
270 PRINTAR(P)iAKG) 
280 PRINTAR(P)iA2(G) 
290 NEXT G 
300 NEXT P 
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305 


END 


310 


F=SQR(RC2+IC2) 


320 


T=ATN(I/R) 


330 


F=l/F 


340 


T=-T 


350 


R=F*COS<T) 


360 


I=F*SIN<T> 


370 


RETURN 
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Complex Arc Cosecant Examples 

?5,8 

THE COMPLEX ARC COSECANT IS - 

.0559826 -.089909 

.0559826 .0899086 

ENTER COMPLEX NUMBER (REAL, IMAGINARY) 

?-2,3 

THE COMPLEX ARC COSECANTIS - 

.150386 .231335 

.150386 -.231335 
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COMPLEX HYPERBOLIC FUNCTIONS 

The last group of trigonometric function programs are the 
complex hyperbolic functions. Programs to solve the complex 
hyperbolic functions are: 

sinh (x + iY) = - i sint (x + iY) 
cosh (x + iY) = cost (x + iy) 
tanh(x + iY) = -£tan£(x + iY) 
coth(x + £Y) = £coti(x + iY) 
csch (x + iy)=l esc £ (x + iY) 
sech (x + iy) = sec £ (x + iY) 

Complex Hyperbolic Sine Program 

5 ' C M P L F X H Y P E R B L 1 C SINE 

10 CLS 

20 PRINT'E'NTER COMPLEX NUNBER < REA 

Li IMAGINARY) " -^~ 

30 INPUT I»R 
40 AR=SIN(R)*< (EXP<1)U+EXP<1 )♦--! ) 

/2) 
50 AI=COS<R)#( (EXP(1)*I-EXP(1)*-I ) 

/2) 
60 PR I NT "THE COMPLEX HYPERBOLIC SI 

NE IS"AI5AR H I" 
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Complex Hyperbolic Sine Examples 

ENTER COMPLEX NUMBER (REAL, IMAGINARY) 

?3,3 

THE COMPLEX HYPERBOLIC SINE IS - 9. 91762 1.42075 

ENTER COMPLEX NUMBER (REAL, IMAGINARY) 

?-3 -3 

THE COMPLEX HYPERBOLIC SINE IS 9. 91762 - 1.42075 
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Complex Hyperbolic Cosine Program 

5 'COMPLEX HYPERBOLIC COSINE 

10 CLS 

20 PRINT" ENTER COMPLEX NUNBER < REA 

L» IMAGINARY)" 
30 INPUT I,R 
40 AR=COS < R ) * ( ( EX P < 1 ) ♦ I +EX P ( 1 ) ♦ - 1 ) 

/2) 
50 AI=<SIN(R>*< (EXP (1)* I-EXP(l)* -I 

)/2) ) 
60 PR I NT "THE COMPLEX HYPERBOLIC CO 

SINE IS"AR5AI" I" 
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Complex Hyperbolic Cosine Examples 

ENTER COMPLEX NUMBER (REAL, IMAGINARY) 

3, 2 

THE COMPLEX HYPERBOLIC COSINE IS -4. 18963 9. 10923 

ENTER COMPLEX NUMBER (REAL, IMAGINARY) 

2, 3 

THE COMPLEX HYPERBOLIC COSINE IS -3.72455 .511823 

Complex Hyperbolic Tangent Programs 

5 'COMPLEX HYPERBOLIC TANGENT 

10 CLS 

20 PRINT" ENTER COMPLEX NUMBER ( REA 

L» IMAGINARY)" 
30 INPUT I»R 
40 D=COS(2*R) + ((EXP(l)i* (2*1 )+EXP( 1 

)t (2*- I ) )/2) 
50 AR=SIN(2*R>/(B> 
60 AI = ( <EXP(1 >*<2*I)+EXP<1 )♦ (2*-I) 

)/2)/(D) 
70 PRINT" THE COMPLEX HYPERBOLIC TA 

NGENT IS"AI;AR" I" 
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Complex Hyperbolic Tangent Examples 

ENTER COMPLEX NUMBER (REAL, IMAGINARY) 

?1,2 

THE COMPLEX HYPERBOLIC TANGENT IS 1.21027 

-.243458 

ENTER COMPLEX NUMBER (REAL, IMAGINARY) 

?5,-3 

THE COMPLEX HYPERBOLIC TANGENT IS .999913 

2.53687E-05 
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Complex Hyperbolic Cotangent Program 

5 'COMPLEX HYPERBOLIC COTANGENT 

10 CLS 

20 PR I NT "ENTER COMPLEX NUNBER < RE A 

L» IMAGINARY)" 
30 INPUT I*R 
40 D=C0S<2*R>-(<EXPUH <2*I)+EXP<1 

)♦ <2*-I> >/2) 
50 AR=SIN<2*R)/(D> 
60 AI=-( (EXP( 1 )♦ (2*1 )-EXP( 1 )♦ <2*-I 

>>/2>/<D> 
70 PR I NT "THE COMPLEX HYPERBOLIC CO 

TANGENT IS"AI?AR" i" 
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Complex Hyperbolic Cotangent Examples 

ENTER COMPLEX NUMBER (REAL, IMAGINARY) 

?3,4 

THE COMPLEX HYPERBOLIC COTANGENT IS .999267 

-4.90118E-03 

ENTER COMPLEX NUMBER (REAL, IMAGINARY) 

?1,10 

THE COMPLEX HYPERBOLIC COTANGENT IS 1.08132 

-.272187 
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Complex Hyperbolic Secant Program 

5 'COMPLEX HYPERBOLIC SECANT 

10 CLS 

20 PRINT "ENTER COMPLEX NUMBER < REA 

L» IMAGINARY)" 
30 INPUT I»R 
40 AR=COS(R)*< (EXP( 1)* I+EXPU>»-I> 

/2) 
50 A1=-(SIN(R)*( (EXP(l)U-EXP(l)*- 

I )/2)) 
60 M=SG1R(AR*2+AI*2) 
70 T=ATN(AI/AR) 
80 FR=-l/M*<COS<-T> ) 
90 FI=1/M*(SIN(-T) ) 

100 PR I NT "THE COMPLEX HYPERBOLIC S 
ECANT IS"FR;FI 
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Complex Hyperbolic Secant Examples 

ENTER COMPLEX NUMBER (REAL, IMAGINARY) 
?1, 2 

THE COMPLEX HYPERBOLIC SECANT IS - .413149 - . 687527 
ENTER COMPLEX NUMBER (REAL, IMAGINARY) 
?2,1 

THE COMPLEX HYPERBOLIC SECANT IS - . 151176 . 226974 
This completes the programs and examples for the 
trigonometric functions. Though each is short, the basic equations 
and examples have been given to aid in understanding. Additional- 
ly, the computational part of the programs can be incorporated into 
larger, more complex programs. 
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Complex Hyperbolic Cosecant Program 

5 'COMPLEX HYPERBOLIC COSECANT 

10 CLS 

20 PRINT" ENTER COMPLEX NUMBER < REA 

L» IMAGINARY)" 
30 INPUT I»R 
A0 AR=SIN(R)*< <EXP<1)* I+EXP( 1>*-I) 

/2) 
50 AI=C0S<R)*(<EXP(1HI-EXP(1)*~I> 

/2> 
60 M=SQR(AR*2+AI*2) 
70 T=ATN(AI/AR) 
80 FR=-l/M*(COS<-T>> 
90 FI=-1/M*(SIN(-T) ) 
100 PRINT "THE COMPLEX HYPERBOLIC C 
OSECANT IS"FI;FR 
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Complex Hyperbolic Cosecant Examples 

ENTER COMPLEX NUMBER (REAL, IMAGINARY) 

?4 5 

THE COMPLEX HYPERBOLIC COSECANT IS -.0103816 

-.0351186 

ENTER COMPLEX NUMBER (REAL, IMAGINARY) 

?5, 4 

THE COMPLEX HYPERBOLIC COSECANT IS 8.80791E-03 

-.0101989 
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TRIANGLES 

The next four programs in this chapter deal with solving 
triangle problems. They are especially useful in surveying work. 
Solving for the unknown sides and angles of oblique triangles 
depends upon using the Law of Sines and Law of Cosines which are: 



Law of Sines- 



1 _ b - c 
Sin A sinB sinC 



Law of Cosines - c 2 = a 2 + b 2 - 2 ab cos C 

Two Sides, Non-included Angle Program 

5 '2 SIDES? NON- INCLUDED ANGLE 

10 CLS 

20 PRINT" ENTER THE TWO SIDES AND T 

HE NON- INCLUDED ANGLE" 
25 PRINT" ENTER FIRST SIDE* SECOND 

SIDE AND ITS ASSOCIATED ANGLE" 

30 INPUT BS>AS»AA 

31 QA=AA 

35 AA=AA*. 01745329 
40 Z=BS*SIN(AA)/AS 

50 BA=ATN(Z/SQR(~Z*Z+1))*57.2957B 
60 CA=180-QA-BA 
65 CA=CA*. 01745329 
70 CS=AS*SIN(CA)/SIN(AA) 
90 PRINT"THE THIRD SIDE IS"CS'S IT 
S ASSOCIATED ANGLE IS"CA*57.295 
78CHR$( 13)" DEGREES AND THE FIRS 
T SIDE'S ASSOCIATED ANGLE IS"BA 
CHR*< 13) "DEGREES." 
100 IF O.A<90 AND AS<BS GOTO 1 10 EL 

SE END 
110 PRINT" A SECOND SOLUTION EXISTS 

FOR THIS TRIANGLE." 
120 BA=180-BA 
130 CA=180-BA-AA*57. 29578 
140 CS=AS*SIN< CA*.01745329)/SIN(AA 
) 
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150 PRINT" FOR THE SECOND SOLUTION, 

THE THIRD SIDE IS"CS"» ITS"CH 

R*< 13) "ASSOCIATED ANGLE IS"CA" 

DEGREES AND THE FIRST SIDE'S"C 

HR$< 13) "ASSOCIATED ANGLE IS"BA 
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32 
KNOWN 



72.2948° 
36.3014 




25 
KNOWN 



Fig. 5-13. A triangle with two known sides and a non-included angle. 

Two Sides, Non-Included Angle Example 

Using Fig. 5- 13, solve for the unknown sides and angles. 
ENTER THE TWO SIDES AND THE NON-INCLUDED ANGLE 
ENTER FIRST SIDE, SECOND SIDE AND ITS ASSOCIATED 
ANGLE 
?35 25 41 

THE THIRD SIDE IS 36.3014, ITS ASSOCIATED ANGLE IS 
72.2948 

DEGREES AND THE FIRST SIDE'S ASSOCIATED ANGLE IS 
66.7052 
DEGREES. 

A SECOND SOLUTION EXISTS FOR THIS TRIANGLE. 
FOR THE SECOND SOLUTION, THE THIRD SIDE IS 16.5283, 
ITS ASSOCIATED ANGLE IS 25.7052 DEGREES AND THE 
FIRST SIDE'S 

ASSOCIATED ANGLE IS 113.295 

This example illustrates that it is possible to have two solutions for 
this triangle. When using degrees, they must be in a decimal 
format, not minutes and seconds. 
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Three Sides Program 

5 '3 SIDES 

10 CLS 

20 INPUT" ENTER THE THREE SIDES OF 
THE TRIANGLE"; A, B»C 

30 D=(A*2+B*2-C*2)/<2*A*B> 

40 CA=-ATN ( D/SQR ( -D*D+ 1 ) ) +1 . 5708 

50 Z=SIN(CA)/C 

60 CA=CA*57. 29578 

70 BA=ATN< <B*Z)/SQR(-B*2*Z*Z+1>) 

80 BA=BA*57. 29578 

90 AA=ATN( <A*Z)/SQR<-A»2*Z*Z+1> ) 

100 AA=AA#57. 29578 

110 PRINT" THE THREE CORRESPONDING 
ANGLES <A,B,C) ARE" CHR$< 13) AA" 
DEGREES" CHR* < 13) BA" DEGREES" CHR 
*< 13 >CA" DEGREES" 
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Fig. 5-14. A triangle with three sides, unknown angles. 

Three Sides Example 

In this example, the three sides are known, and you must 
determine the three included angles (use Fig. 5-14). 
ENTER THE THREE SIDES OF THE TRIANGLE? 
? 11, 14, 16 

THE THREE CORRESPONDING ANGLES (A, B, C) ARE 
42.3674 DEGREES 
59. 0558 DEGREES 
78. 5771 DEGREES 
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Two Angles and Opposite Side Program 

5 '2 SIDES, OPPOSITE ANGLE 

10 CLS 

20 INPUT" ENTER THE TWO ANGLES (AiB 

) AND THE OPPOSITE SIDE" SAA»BA» 

CS 
30 CA=180-AA-BA 
40 BS=CS*SIN(BA/57. 29578) /SIN (AA/5 

7.29578) 
50 CS=CS*SIN(CA/57. 29578) /SIN (AA/5 

7.29578) 
60 PR I NT "THE THIRD ANGLE IS"CA"DEG 

REES AND THE TWO UNKNOWN OPPOSI 

TE"CHR*( 13) "SIDES ARE"BS" AND" CS 
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12.5169 




12.5614 



KNOWN 



Fig. 5-15. A triangle with two known sides, opposite angles. 

Two Angles and Opposite Side Example 

Using Fig. 5-15, solve for the remaining angle and two sides. 
ENTER THE TWO ANGLES (A,B) AND THE OPPOSITE SIDE? 
? 23, 78, 5 

THE THIRD ANGLE IS 79 DEGREES AND THE TWO UN- 
KNOWN OPPOSITE SIDES ARE 12. 5169 AND 12. 5614. 
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Two Sides, Included Angle Program 

5 '2 SIDES, INCLUDED ANGLE 
10 CLS 

20 INPUT" ENTER THE TWO SIDES <A»B) 
AND THE INCLUDED ANGLE" ; AS, BS, 

CA 
30 CS=SQR < AS ♦ 2+BS *2- < 2*AS*B3*C0S ( C 

A/57.29578)) ) 
40 AA=ATN<AS*SIN<CA/57. 29578) /CS/S 

QR ( -AS*SIN ( CA/57 . 29578 ) / CS*AS*S 

IN(CA/57.29578)/CS+l) ) 
50 AA=AA*57. 29578 
60 BA=180-AA-CA 
70 PRINT"THE THIRD SIDE IS"CS"."CH 

R*<13)"THE CORRESPONDING ANGLES 
FOR THE TWO KNOWN SIDES ARE"AA 

CHR$( 13)" AND" BA" DEGREES." 
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Fig. 5-16. A triangle with two known sides, included angle. 

Two Sides, Included Angle Example 

Solve Fig. 5-16 for the unknown side and angles. 
ENTER THE TWO SIDES (A,B) AND THE INCLUDED ANG- 
LE? 
? 3,5,39 

THE THIRD SIDE IS 3. 26889. 

THE CORRESPONDING ANGLES FOR THE TWO KNOWN 
SIDES ARE 35.2787 AND 105. 721DEGREES. 
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GREAT CIRCLE CALCULATIONS 

The last program in this chapter will compute the distance in 
nautical miles and bearing (referenced to True North) between two 
points on the Earth. Information for this program must be entered 
as decimal degrees. Plus, Southern Latitudes and Eastern Lon- 
gitudes must be entered as negative numbers. 

The equations to manually calculate the distance and bearing 
are: 

A = (sin LAI • sin LAZ + COS LAI • COS LAZ • COS (L02-L01) 
D = 60 -Cos- 'A 

r ^- sin LA2 - (sin LAI cos D/6Q) 
sin(D/60)cosLAl 

H = cos~ 1 C 
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Great Circle Example 

Compute the great circle distance and bearing from San 
Francisco to Tokoyo. 

Latitude Longitude 

SanFrancisco 37.82 N 122.42 W 

Tokoyo 35.67N 139. 75 E 

WHEN ENTERING DATA, SOUTHERN LATITUDES AND 
EASTERN LONGITUDES MUST BE ENTERED AS NEGA- 
TIVE NUMBERS. 

ENTER LATITUDE AND LONGITUDE (LAT, LON) OF THE 
STARTING POINT 
? 37.82, 122.42 

ENTER LATITUDE AND LONGITUDE (LAT, LON) OF THE 
DESTINATION 
? 35.67, -139.75 

THE DISTANCE IS 4456.89 NAUTICAL MILES AT A BEARING 
OF 
303. 292 DEGREES. 

Again note that the Eastern Longitude of Tokoyo was entered 
as a negative number. 
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AUTO 

CLEAR 
CLEAR N 

CLOAD 



CLOAD? 



Appendix A 

Statements and 
Commands in BASIC 

turn on automatic line num- 
bering. Can be set for dif- 
ferent increments. 
set numerics to zero, 
strings to null, 
same as CLEAR but sets 
aside a certain amount of 
memory for strings (N). 
load program from cassette 
player. If no file-name is 
specified, the current pro- 
gram will be loaded, 
compares program on tape 
with program in memory. If 
programs are not equal 
"BAD" will be printed on 
video display. 
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CONT 
CSAVE 



DATA 
DELETE 

DIM 

EDIT 
END 

ERROR 



FOR— TO- 
FOR— TO-STEP- 



GOTO 



continue after BREAK or 
STOP in execution, 
saves current program on 
cassette tape, a file-name 
should be used so that 
program can be compared 
with what is in memory, 
hold DATA for accessby 
READ statement, 
delete program line or lines 
specified by INPUT 
statement. 

allocate storage for dimen- 
sional array with specified 
size per dimension, 
puts computer in EDIT 
mode for line specified, 
end execution of program, 
returns computer to 
command mode, 
if an error is encountered 
in the program, the com- 
puter prints the offending 
line number and the error 
code for that line, 
open a FOR-NEXT loop, 
open a FOR-NEXT loop 
with optional STEP. If no 
step is specified an incre- 
ment of one is used, 
branch to a specified line 
number in the program. 
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GOSUB 



IF-TH EN-ELSE 



INPUT 
INPUT #-1 
LET (var-exp) 



LIST 
NEXT var. 



NEW 



ON-GOTO 



ON-GOSUB 



branch to a sub-routine 
beginning with the speci- 
fied line number. 

test expression; if true, 
execute the statement and 
jump to the next program 
line. If it's false, the ELSE 
statement will be executed, 
await input from keyboard, 
input from cassette #1. 
assigns the value of an ex- 
pression to variable. The 
LET is optional with LEVEL 
II BASIC. 

list all program lines, or a 
specified range of lines, 
close the FOR - NEXT loop. 
The variable may be 
omitted, or a variable list 
may be used to close 
nested loops. 

deletes the entire program 
in memory and returns 
to command mode, 
branches to a specified 
line number after 
evaluating the expression. 
Otherwise go to the next 
statement. 

same as ON - GOTO but 
branches to the appropriate 
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PRINT exp. 



PRINT #-1 
PRINT 



PRINT USING 

RANDOM 

REM 

READ 

RESTORE 

RESUME 
RETURN 
RUN 



STOP 



GOSUB statement, 
outputs to the video display 
the value of exp, which 
can be a string expression, 
a numeric, a constant or a 
list of items, 
output to cassette #1 
begins printing at a 
specified area on the video 
display; a PRINT modifier. 
PRINT formatspecifier, out- 
put exp in a form 
specified by a string field, 
reseeds the RANDOM 
number generator, 
remark indicator. The rest 
of line is ignored, 
assigns values to specified 
variables starting with 
the current DATA element, 
resets the data pointer to 
the first item in the first 
DATA statement, 
returns from the error 
routine to the line specified, 
branch statement following 
last GOSUB executed, 
execute the program. If no 
line is specified, execution 
begins with the program's 
lowest line number, 
stop the program's exe- 
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SYSTEM 



TAB 

TROFF 
TRON 



cution and print a break 
message containing 
current line number. Pro- 
gram execution can resume 
with CONT. 

enter the monitor code for 
loading a machine lan- 
guage program from cas- 
sette tape player, 
print modifier; begins 
printing at specified TAB 
area of the video display, 
turn off the line-numbers 
trace. 

turn on the line numbers 
trace (very useful when 
debugging). 



String Functions 

ASC 

CHR$ 



FRE 
INKEY$ 



returns the ASCII code of 

first character in string 

argument. 

returns a one-character 

string defined by a code. 

If a contol function is 

specified by code, 

that function is acti-vated. 

returns amount of memory 

available for string storage. 

strobes the keyboard and 
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LEN 
LEFT$ 



MID$ 

MID$ (cont) 
RIGHTS 



STR$ 



STRINGS 



VAL 



returns a one- character 
string corresponding to 
the key pressed during 
strobe. 

returns the numeric 
length of string, 
returns the first n 
characters of a string 
(when n is a numeric 
expression). 

returns the substring of a 
string with length n and 
starting at position 
minastring(wherenand m 
are numeric expressions), 
returns the last n 
characters of a string 
(where n is a numeric 
expression), 
returns a string 
representation of an 
evaluated argument, 
returnsasequence of char- 
acters using the first char- 
acter. PRINTSTRINGS$ 
(50, "*"), for example, will 
print 50 asterisks on the 
video display, 
returns a numeric value 
corresponding to a 
numeric-valued string, such 
as PRINT VAL (A$+W$). 
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Appendix B 
Powers of Two 



t 



2' THREE BYTE BINARY 






1 


00000000 


00000000 


00000001 


1 


2 


00000000 


00000000 


00000010 


2 


4 


00000000 


00000000 


00000100 


3 


8 


00000000 


00000000 


00001000 


4 


16 


00000000 


00000000 


00010000 


5 


32 


00000000 


00000000 


00100000 


6 


64 


00000000 


00000000 


01000000 


7 


128 


00000000 


00000000 


10000000 


8 


256 


00000000 


00000001 


00000000 


9 


512 


00000000 


00000010 


00000000 


10 


1024 


00000000 


00000100 


00000000 


11 


2048 


00000000 


00001000 


00000000 


12 


4096 


00000000 


00010000 


00000000 


13 


8192 


00000000 


00100000 


00000000 


14 


16384 


00000000 


01000000 


00000000 


15 


32768 


00000000 


10000000 


00000000 


16 


65536 


00000001 


00000000 


00000000 


17 


131072 


00000010 


00000000 


00000000 


18 


262144 


00000100 


00000000 


00000000 


19 


524288 


00001000 


00000000 


00000000 


20 


1048576 


00010000 


00000000 


00000000 


21 


2097152 


00100000 


00000000 


00000000 


22 


4194304 


OlOOOOOO 


00000000 


00000000 


23 


8388608 


10000000 


00000000 


00000000 
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Appendix C 

Hexadecimal-Decimal 

Integer Conversion 



The following table provides for direct conversions 
between hexadecimal integers in the range 0-FFF and 
decimal integers in the range 0-4095. For conversion of larger 
integers, the table values may be added to the following 
figures: 

Hexadecimal Decimal Hexadecimal Decimal 




01000 


4 096 


20 000 


131 072 


02 000 


8 192 


30 000 


196 608 


03 000 


12 288 


40 000 


262 144 


04 000 


16 384 


50 000 


327 680 


05 000 


20 480 


60 000 


393 216 


06 000 


24 576 


70 000 


458 752 


07 000 


28 672 


80 000 


524 288 


08 000 


32 768 


90 000 


589 824 


09000 


36 864 


AOOOO 


655 360 


0A000 


40 960 


60 000 


720 8% 


OB 000 


45 056 


CO 000 


786 432 


ocooo 


49 152 


DO 000 


851 968 


0D000 


53 248 


EOOOO 


917 504 


OEOOO 


57 344 


FOOOO 


983 040 


OF 000 


61440 


100 000 


1 048 576 


10 000 


65 536 


200 000 


2 097 152 


HOOO 


69 632 


300 000 


3 145 728 


12 000 


73 728 


400 000 


4 194 304 


13 000 


77 824 


500 000 


5 242 880 


14 000 


8! 920 


600 000 


6 291 456 


15 000 


86 016 


700 000 


7 340 032 


16 000 


90112 


800 000 


8 383 608 


17 000 


94 208 


900 000 


9437 184 


18 000 


98 304 


AOOOOO 


10 485 760 


19000 


102 400 


BOO 000 


1 1 534 336 


1A000 


106 4% 


COO 000 


12 582 912 


IB 000 


110 592 


DOOOOO 


13 631 488 


ICOOO 


114 688 


£00 000 


14 680 064 


1O000 


118 784 


FOOOOO 


15 728 640 


IE 000 


122 880 


1000 000 


16 777 216 


IF 000 


126 976 


2000 000 


33 554 432 
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